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ABSTRACT 

Electronic  Properties  of  High- Temperature  Superconductors.  (August  1989) 
Brent  Armand  Richert,  B.S..  U.  S.  Air  Force  Academy; 

M.S.,  University  of  New  Mexico 
Chair  of  Advisory  Committee:  Dr.  Roland  E.  Allen 


We  have  developed  a  semiempirical  tight-binding  model  for  the  elec¬ 
tronic  energy  bands,  the  local  and  total  densities  of  states,  and  the 

v  v 

atomic  valences  in  the  high-temperature  superconductors  Lai.soSro.isCuO.}, 

Oo'- 

VBalCu30'7,  Bi'jSrjCyps,  BiVCaSrVCujO^1,  Tr2Ba2CuOs,  T^CaBa^Cu^oU 


Tl2Ca2Ba2Cu30'io,  TlCaVB^C^On,  BaPbfp.75Bio.25O3,  and  Ba^.^Kp.^BiO's 

A  single  tight-binding  model,  with  fully  transferable  parameters,  provides  a 
good  description  of  the  electronic  structures  of  all  these  materials. 

Calculations  of  the  changes  in  electronic  properties  associated  with  atomic 

- )  v ,  , ,  v.  v ,  «,  v  •  y  '  .  \a 

substitutions  in  YBa2Cu3'07,  Bi2CaSr2Cu20g,  and  Tl£CaBa$Cu2Otf  give  results 


3  ' 


in  agreement  with  expected  chemical  trends  and  consistent  with  observed 
changes  in  the  superconducting  properties.  For  example,  substitution  ofPb  for 

w-~v'  r  o. 

Bi  in  Bi.l/CaSr2Cu208  increases  the  concentratioruof  hole  carriers  within  the 
Cu02/! planes.  Similarly,  doping  with  Hg  or  Pb  in  TtA/CaBa2Cu20g  also  affects 

the  carrier  concentration,  with  Hg  creating  holes  and  Pb  destroying  them. 

v-  ,  K  v  .  s.  v 

Oxygen  vacancies  in  both  Lax  85Sro.i5CuO;4_y"',  and  \Ba2Cu307-y  act  as 
electron  donors.  This  is  consistent  with  the  observations  that  oxygen  vacancies 
degrade  the  superconductivity  and  metallic  conductivity  in  these  materials.  Lan¬ 
thanum  vacancies  in  La2-xCu04  donate  holes,  giving  the  same  electronic  effect- 


IV 


as  doping  with  divalent  metal  atoms  or  excess  oxygen  in  initially  stoichiometric 
La2Cu04- 

We  propose  a  specific  excitonic  mechanism  for  high-temperature  supercon¬ 
ductivity,  which  requires  insulating  metal-oxide  layers  adjacent  to  the  supercon¬ 
ducting  planes.  The  attractive  pairing  interaction  between  carriers  is  mediated 
by  two-dimensional  excitons  in  the  insulating  region.  No  specific  magnetic  prop¬ 
erties  in  the  metallic  region  are  needed  for  this  mechanism,  which  is  applicable 
to  all  currently-known  high- 71.  materials.  Reasonable  estimates  of  the  physical 
parameters  lead  to  a  pairing  interaction  that  is  sufficient  to  give  high- 71  super¬ 
conductivity.  The  theory  predicts  a  metal-oxide  (e.g.,  BaO)  exciton  that  may  be 
observable.  In  addition,  the  theory  predicts  that  high-temperature  superconduc¬ 
tivity  will  not  be  found  in  materials  such  as  Cai_xKxCu02  or  Cai_xVxCu02, 
since  these  materials  lack  the  insulating  metal-oxide  layers  necessary  to  support 
the  excitons. 
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CHAPTER  I 

INTRODUCTION 

The  discovery  of  high-temperature  superconductivity  by  Bednorz  and 
Muller1  has  been  followed  by  an  intense  search  for  new  materials,  and  an  equally 
intense  investigation  of  their  structural  and  electronic  properties.  Numerous 
systems  have  been  discovered,  with  a  current  record  superconducting  transition 
temperature.  Tc,  of  125  K.  The  discovery  of  high-temperature  superconductors 
without  copper,  principally  Bat_xKxBi03  with  Tc  ==  34  K,  has  placed  stringent 
constraints  on  the  development  of  a  comprehensive  theory  for  high-temperature 
superconductivity.  As  described  in  the  following  chapters,  we  have  addressed 
four  theoretical  problems:  (1)  the  basic  electronic  structure  of  high-Tc  super¬ 
conductors:  (2)  the  effects  of  atomic  substitutions:  (3)  the  effects  of  atomic 
vacancies;  and  (4)  the  mechanism  of  high-temperature  superconductivity. 

Chapter  II  provides  an  overview  of  the  literature  detailing  the  experimental 
discoveries  and  characterizations  of  high- temperature  superconductors.  We 
attempt  to  give  proper  credit  for  original  discoveries,  but  much  of  the  literature 
includes  concurrent  research  and  confirmatory  results  wrhich  are  also  important. 
Emphasis  is  placed  on  those  experiments  that  are  directly  related  to  our  own 
research,  but  results  that  are  needed  to  appreciate  the  complexity  of  these 
materials  are  also  included. 

Chapter  III  describes  the  development  of  a  semiempirical  tight-binding 
model  which  has  been  used  to  calculate  electronic  energy  bands,  local  and 
total  densities  of  states,  and  atomic  valences.  We  find  that  a  single  model, 

The  journal  model  used  is  Physical  Review  B. 
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with  fully  transferable  parameters,  provides  a  good  description  of  the  electronic 
structure  of  every  material  examined,  in  the  sense  that  results  in  the  local-density 
approximation  are  satisfactorily  reproduced. 

The  results  of  electronic  structure  calculations  based  on  this  model  are 
presented  in  Chapter  IV.  The  copper-oxide  superconductors  La185Sr0.i5CuO4, 
YBa2Cu3C>7,  Bi2Sr2Cu06,  Bi2CaSr2Cu20g,  Tl2Ba2Cu06,  Tl2CaBa2Cu208, 
Tl2Ca2Ba2Cu30in,  and  TlCa3Ba2Cu40ui  and  the  bismuth-oxide  superconduc¬ 
tors  BaPbo.75 Bi0.25O3  and  Bao.6Ko.4Bi03,  have  been  examined.  Although  these 
materials  have  very  different  crystal  and  magnetic  structures,  common  features 
in  the  electronic  structures  provide  strong  clues  regarding  the  mechanism  of 
superconductivity. 

We  present  in  Chapter  V  a  study  of  atomic  substitutions  in  YBa2Cu307, 
Bi2CaSr2Cu208,  and  Tl2CaBa2Cu208,  to  determine  the  modifications  in  the 
electronic  structure.  Substitutions  on  the  Ba,  Cu,  and  0  sites  were  considered 
for  the  case  of  YBa2Cu807.  In  the  remaining  materials,  we  studied  doping  on 
the  Bi  or  T1  site.  The  results  are  compared  with  the  experimentally-observed 
changes  in  the  superconductivity  of  these  materials. 

A  critically  important  factor  in  these  superconductors  is  the  oxygen  content. 
Calculations  of  the  electronic  effects  of  oxygen  vacancies  in  Laj  85Sro.15CuO4._y 
and  YBa2Cu3C>7_y  are  reported  in  Chapter  VI.  The  modification  of  the  density 
of  states  and  the  shift  in  the  Fermi  energy  were  calculated  for  0  <  y  <  1.0.  The 
effect  of  lanthanum  vacancies  in  undoped  La2CuC>4  is  also  presented. 

The  most  important  theoretical  question  is  the  specific  mechanism  of  high- 
Tc  superconductivity.  In  Chapter  VII  we  propose  a  mechanism  that  involves 
two-dimensional  excitons  in  the  insulating  metal-oxide  layers  (e.g.,  LaO,  BaO, 
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or  SrO)  adjacent  to  the  superconducting  layers  (e.g.,  CuC>2  or  Bi02).  This 
mechanism  requires  no  specific  magnetic  properties  for  the  conduction  region, 
and  is  applicable  to  all  currently-known  high-temperature  superconductors.  We 
have  calculated  matrix  elements  for  the  interaction  between  charge  carriers  and 
the  metal-oxide  excitons,  and  have  estimated  the  attractive  interaction  leading 
to  superconductivity. 
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CHAPTER  II 

REVIEW  OF  EXPERIMENTS 
Stoichiometry  and  Crystal  Structure 

The  discovery  of  superconductivity  in  the  La-Ba-Cu-0  system  by  Bed- 
norz  and  Miiller1  in  1986  initiated  an  explosion  of  research  activity  in  high- 
temperature  superconductivity.  Their  polycrystalline  samples,  of  nominal  com¬ 
position  BaxLa5_xCu50i5_y ,  indicated  an  onset  of  superconductivity  at  Tc  ^ 
30  K  in  resistivity  measurements.1  The  superconducting  phase  with  Tc  =  35  K 
was  identified  as  the  solid  solution  La2-xBaxCu04,  x  ==  0.15,  which  has  the 
tetragonal  K2NiF4  structure  at  room  temperature.2  This  oxygen-defect  per- 
ovskite  is  characterized  by  Cu02  planes  surrounded  by  layers  of  LaO,  giving  two 
formula  units  per  bet  cell.  The  transition  temperature  was  optimized  by  replac¬ 
ing  Ba  with  Sr,  giving  Tc  =  38  K  for  La185Sro.i5Cu04.3'4  Neutron  diffraction 
studies  of  this  phase  indicate  a  transition  to  an  orthorhombic  crystal  structure 
below  200  K,  associated  with  a  buckling  of  the  Cu02  planes  and  a  doubling 
of  the  unit  cell.5  The  corner-sharing,  tilted  Cu08  octahedra  have  an  in-plane 
copper-oxygen  bond  length  of  1.89  A  and  an  out-of- plane  bond  length  of  2.41  A.5 

Stoichiometry  has  been  found  to  play  a  critical  role  in  several  different 
respects.  The  superconducting  properties  are  highly  sensitive  to  the  oxygen 
content,  with  the  highest  Tc  resulting  from  full  oxygen  occupancy.3’4  Excess 
Sr  doping  with  x  >  0.15  apparently  leads  to  charge  compensation  through  the 
formation  of  oxygen  vacancies,  with  a  depression  of  Tc .3’6,7  Annealing  under 
high  oxygen  pressure  suppresses  the  formation  of  oxygen  vacancies,  and  allows 
superconductivity  at  Tr  ==  36  K  in  La2_xSrxCu04  up  to  x  =  0.24. 8  However, 
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the  transition  temperature  decreases  for  x  >  0.24,  with  no  superconductivity 
observed  beyond  x  —  0.32,  although  metallic  conductivity  persists.8 

The  undoped  material  La2Cu04  has  an  orthorhombic  structure  at  room 
temperature,9  with  a  possible  monoclinic  distortion  reported  below  10  K.10 
This  parent  compound  exhibits  superconductivity  of  a  filamentary  nature.  The 
small  volume  fraction  (much  less  than  1%)  of  superconducting  region  was  at¬ 
tributed  to  either  excess  oxygen  or  lanthanum  vacancies  within  the  material.11 
The  lanthanum-deficient  material  La2_xCu04  shows  a  sharper  resistive  tran¬ 
sition  and  a  greater  diamagnetic  susceptibility  than  the  nominal  La2C'u04 
superconductor,12  indicating  that  La  vacancies  can  indeed  contribute  to  the 
superconductivity. 

The  discovery  of  superconductivity  in  the  Y-Ba-Cu-0  system13  with  a  Tc 
of  93  K  surpassed  the  previous  theoretical  limits  on  Tc  and  broke  the  temperature 
barrier  of  77  K,  the  boiling  point  of  liquid  nitrogen.  The  superconducting  phase 
was  rapidly  isolated14-16  as  YBa2Cu307,  often  called  the  1:2:3  phase  to  refer 
to  the  metal  stoichiometry.  This  material  has  an  orthorhombic  crystal  structure 
(a  =  3.856  A,  b  =  3.870  A,  and  c  =  11.666  A)  characterized  by  double  layers 
of  Cu02,  separated  by  Y  ions  and  surrounded  by  BaO  layers,  stacked  with 
one-dimensional  chains  of  CuO. 17-21  The  chains  form  a  linear  coordination  of 
copper  and  oxygen  atoms  along  the  b  axis  of  the  crystal,  with  a  bond  length  of 
1.93  A.18  The  Cu02  layers  are  distorted  from  perfect  planes  by  displacements 
of  the  oxygen  atoms  along  the  c  axis  toward  the  Y  ion.17  The  plane  copper  sites 
are  five-coordinated  with  oxygen,  with  an  average  in-plane  interatomic  distance 
of  1.95  A  and  an  out-of-plane  distance  of  2.38  A.18  These  bond  lengths  are 
comparable  to  the  Cu-0  distances  in  Lai .ssSro.isCuC^.  The  small  difference  in 


6 


the  a  and  b  lattice  parameters  allows  the  formation  of  twin  boundaries  along  the 
1101  direction  in  the  orthorhombic  phase.22 

A  high-temperature  orthorhombic-to-tetragonal  phase  transition  near 
700  °C  in  YBa2Cu3C>7  is  found  to  be  an  order-disorder  transition  in  which 
the  oxygen  atoms  in  the  CuO  chains  are  disordered  onto  the  normally  vacant 
adjacent  sites.23  Oxygen  is  lost  from  YE^CusOr  under  ambient  atmosphere  at 
temperatures  over  400  °C.  and  the  orthorhombic-to-tetragonal  transition  occurs 
near  an  oxygen  content  of  6. 5. 23  The  decrease  in  oxygen  stoichiometry  intro¬ 
duces  vacancies  on  the  oxygen  chain  sites.24  Other  structural  effects  are  also 
associated  with  oxygen  vacancies,  such  as  an  increase  in  the  c  lattice  parameter 
and  an  increase  in  the  Cu-0  out-of-plane  interatomic  distance  with  decreasing 
oxygen  content.25 

Oxygen  vacancies  in  YBa2Cu307_y  have  a  dramatic  effect  on  the  super¬ 
conductivity.  The  oxygen  vacancies  lead  to  a  depression  of  Tc  from  above  90  K 
to  about  55  K  for  0.1  _  y  <  0.5.  and  a  loss  of  the  superconductivity  and 
even  metallic  conductivity  for  y  >  0.5. 26,27  The  plateau  in  Tc  near  55  K  for 
0-3  <  y  '£  0.5  is  associated  with  a  long-range  ordering  of  the  oxygen  vacancy 
sites  along  the  6  axis,  parallel  to  the  chains,  and  a  short  correlation  length 
along  the  a  axis.28,29  Speculation  on  the  importance  of  the  integrity  of  the  CuO 
chains  for  high-temperature  superconductivity  centered  on  such  issues  as  one¬ 
dimensional  superconductivity  along  the  chains,30  oxygen-defect  enhancement 
of  pairing  for  superconductivity,31  and  coupling  between  the  Cu02  planes  and 
CuO  chains.23,25  However,  further  experiments  have  indicated  that  the  one¬ 
dimensional  nature  of  the  CuO  chains  is  unimportant  for  high-temperature 
superconductivity.  Experiments  with  the  systems  YBa2-xLaxCu307_y  and 
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Y1_xCaxBa2Cu3  07-y  demonstrate  that  the  transition  temperature  may  be  var¬ 
ied  without  changing  the  oxygen  content.32,33  The  Tc  is  correlated  with  the 
carrier  density  within  the  CuC>2  planes,  while  charges  associated  with  the  CuO 
chains  are  localized  and  do  not  contribute  to  superconductivity.32  In  particular, 
the  material  Yo.85Cao.isBa2Cu306  is  superconducting  with  Tc  =  50  K,  but  has 
a  tetragonal  structure  and  complete  oxygen  vacancies  in  the  chain  region.34  The 
oxygen  in  the  chains  thus  serves  primarily  as  an  electron  reservoir  for  the  rest  of 
the  structure.33  In  the  same  manner  that  a  saturation  and  decrease  in  Tc  is  seen 
in  La2_xSrxCu04  for  x  >  0.24, 8  doping  with  Ca  in  Yi _xCaxBa2Cu307  leads  to 
a  depression  of  Tc  from  93  K  to  =  80  K  as  the  carrier  concentration  increases.34 

A  second  superconducting  phase  in  the  Y-Ba-Cu-0  system  was  first 
observed  as  a  defect  structure  in  YBa2Cu307  that  increased  the  c  axis  from 
11.7  A  to  13.6  A.36,37  This  phase,  with  Tc  —  81  K.  was  isolated  in  thin  film 
samples38,39  as  YBa2Cu40s,  with  an  orthorhombic  unit  cell  and  a  doubled  c 
axis  of  27.2  A.  The  structure  includes  two  CuO  chains,  linked  together  in  an 
edge-sharing,  square-planar  CUO4  coordination.40  The  remaining  structure  is 
unchanged  from  that  of  YBa2Cu307-  Intergrowths  of  single-chain  YBa2Cu307 
and  double-chain  YBa2Cu40g  can  also  form,  resulting  in  a  Y2Ba4Cu7015 
phase.41 

Bulk  samples  of  YBa2Cu4  0s  have  been  formed,41,42  along  with  a  series 
of  samples  with  rare  earths  replacing  Y.42  These  samples  all  exhibit  high- 
temperature  superconductivity  with  Tc  ranging  from  57  to  81  K.  The  oxygen 
stoichiometry  of  this  phase  is  relatively  constant  compared  to  the  \Ba2Cu3O7 
phase,  perhaps  because  of  the  additional  copper  coordination  of  the  oxygens 
in  the  double  chains.40-42  The  oxygen  content  is  thermally  stable  up  to  a 
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temperature  of  850  °C,  beyond  which  the  structure  decomposes.41  Thus,  it  is 
difficult  to  study  the  variation  of  Tc  with  oxygen  content  in  YBa2Cu408. 

The  next  series  of  discoveries  raised  the  record  transition  temperature 
for  the  copper-oxide  superconductors  to  over  100  K.  Michel  et  al.i3  reported 
that  partial  substitution  of  bismuth  on  the  La/Sr  site  in  La2_xBizSrx_zCu04 
increased  the  transition  temperature  Tc  from  38  K  to  42  K.  They  also  reported44 
a  new  family  of  superconductors  in  the  Bi-Sr-Cu-0  system  with  Tc  from  7  to 
22  K.  The  absence  of  rare  earth  elements  distinguished  this  system  from  the  other 
cuprate  superconductors.  The  superconducting  phase  was  tentatively  identified 
as  Bi2Sr2Cu207_y ,  with  a  substructure  of  Bi202  units  stacked  with  the  copper- 
oxide  perovskite  substructures.44  The  discovery  of  superconductivity  with  onset 
Tys  up  to  114  K  in  the  Bi-Ca-Sr-Cu-0  system40,46  extended  this  system  into 
true  high-temperature  superconductivity. 

The  phase  with  Tc  =  85  K  was  identified47,48  as  Bi2CaSr2Cu208.  Inter¬ 
growths  of  a  higher-Tc  phase  gave  magnetic  and  resistive  onset  temperatures 
above  110  K,  but  the  resistance  of  the  mixed-phase  sample  did  not  go  to  zero 
above  85  K.  The  85  K  phase  was  found  to  have  a  layered  structure.48-02  There 
are  two  Cu02  planes,  separated  by  a  Ca  ion  and  surrounded  by  SrO  layers.  This 
substructure  is  then  layered  with  two  BiO  planes,  giving  the  stacking  sequence 
Bi0-Sr0-Cu02-Ca-Cu02-Sr0-Bi0  along  the  c  axis.  A  solid  solution  for  Sr 
and  Ca  allows  mutual  substitution  on  these  sites,  with  Tc  optimized  for  the 
indicated  stoichiometry.47-50  The  crystal  structure  is  pseudotetragonal,49  with 
a  =  3.81  A  and  c  =  30.52  A,  but  with  a  long-range  incommensurate  modula¬ 
tion  of  approximately  5v  2  cells  27  A)  along  the  [110  direction.48-00  Instead 
of  the  twinning  observed  in  YBa2Cu307,  90°  twist  boundaries  are  observed  as 
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stacking  faults  perpendicular  to  the  c  axis.53,54  The  relatively  large  separation 
(=  3.3  A)  between  the  BiO  double  layers  gives  a  micaceous  nature  to  the  bulk 
material.48'49 

The  incommensurate  structural  modulation  is  associated  with  disorder 
within  the  BiO  layers,48  52,54  although  additional  disorder  is  observed  in 
the  SrO  and  Cu02  layers.35  Electron  microscope31,56  and  scanning  tunneling 
microscope37  images  detail  the  Bi  atom  superstructure,  giving  some  support  to 
a  structural  model  that  requires  a  vacant  Bi  site  every  nine  or  ten  unit  cells  to  fit 
the  modulation.33,3'  However,  a  model  for  excess  oxygen  atoms  within  the  BiO 
bilayer  can  account  for  the  superstructure38  and  the  excess  oxygen  stoichiome¬ 
try  (measured  by  thermogravimetric  analysis  as  8.25  per  formula  unit49).  This 
model  compares  B^CaS^C^Og  with  the  isostructural  material  Bi2Sr3pe209, 
which  has  a  commensurate  modulation.  This  isostructural  compound  accom¬ 
modates  an  extra  oxygen  every  ten  cells  to  stabilize  the  strained  BiO  bilayer. 
In  Bi2CaSr2Cu208,  one  extra  oxygen  for  every  9.52  bismuth  atoms  only  ap¬ 
proximately  fits  the  copper-oxygen  sublattice,  leading  to  the  incommensurate 
distortion.58 

Bi2CaSr2Cu2  08  is  relatively  insensitive  to  oxygen  loss  compared  to 
YBa2Cu307.  Heating  in  argon  leads  to  an  oxygen  loss  of  0.1  per  formula 
unit,  but  superconductivity  persists.49  Annealing  in  oxygen  leads  to  a  slight 
excess  in  the  oxygen  stoichiometry,49  perhaps  associated  with  the  incommensu¬ 
rate  modulation.58 

The  Tc  =  114  K  phase  proved  to  be  difficult  to  isolate.  Although  clear 
resistive  and  magnetic  onsets  of  superconductivity  are  observed  near  114  K, 
a  zero-resistance  state  above  100  K  could  not  initially  be  found.  The  finite 
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resistance  was  attributed  to  incomplete  grain  connectivity  for  the  higher-Tc 
phase,  which  grew  in  thin  shells  around  the  85  K  phase.59  Annealing  the  samples 
in  oxygen  at  885  °C  for  several  days  enhanced  the  Meissner  effect  of  the  higher- 
Tc  phase,39  and  produced  zero  resistance  at  107  K.  Addition  of  Pb  to  this 
material  also  allowed  zero  resistance  up  to  107  K. 50,60,61  Pb  substitutes  for  Bi 
and  stabilizes  the  higher- Tc  structure.50,61 ,62 

Bi2Ca2Sr2Cu30io  was  identified  as  the  114  K  superconductor.63,64  The 
structure  of  this  phase  is  the  same  as  for  B^CaS^C^Og,  but  with  an  additional 
CaCuC>2  slab,  giving  a  c  axis  of  37.1  A  and  a  stacking  sequence  of  BiO-SrO- 
Cu02-Ca-Cu02-Ca-Cu02-Sr0-Bi0.64  Stacking  faults  in  this  phase  occur  by 
intergrowths  of  the  85  K  phase,  which  have  double  rather  than  triple  CuC>2  lay¬ 
ers.  These  intergrowths  have  prevented  a  definitive  determination  of  the  atomic 
substructure  and  possible  incommensurate  modulation  in  Bi2Ca2Sr2Cu3Oio. 
However,  disorder  in  the  BiO  layers  is  expected  to  occur  in  this  phase.64 

A  third  phase  B^S^CuOg  has  a  Tc  of  only  6  K. 64-66  This  phase  is  similar 
to  that  proposed  by  Michel  et  a/.44  The  bet  structure  with  a  =  3.81  A  and 
c  =  24.61  A  h  as  a  single  CuC>2  plane  per  unit  cell,65  similar  to  La2-xSrxCu04- 
The  plane  is  sandwiched  by  SrO  layers  and  stacked  with  double  BiO  layers  to  give 
the  c  axis  ordering  Bi0-Sr0-Cu02-Sr0-Bi0.  The  CuOg  octahedra  are  tilted 
as  in  La2-xSrxCu04,  distorting  the  copper-oxide  region  from  a  perfect  plane.66 
A  small  orthorhombic  or  monoclinic  distortion  is  observed  along  with  90° 
twist  boundaries  about  0011.66  An  incommensurate  superstructure  is  found  in 
B^S^CuOe,  which  is  again  related  to  disorder  in  the  BiO  layers.64,66  Vacancies 
on  the  Sr  site  have  been  observed,63  which  may  account  for  some  variations  in 
Tc  reported  for  this  phase.64-66 
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Another  set  of  copper-oxide  materials  containing  T1  has  yielded  the  highest 
Tc' s.  Superconductivity'  in  the  Tl-Ba-Cu-0  system  was  seen  above  80  K.67  The 
phase  with  Tc  =  83  K  was  identified  as  T^I^CuOe,  and  is  structurally  very 
similar  to  B^S^CuOg,  although  the  CuC>2  layers  are  planar.68  The  addition 
of  Ca  to  this  system  increases  Tc  over  110  K. 68,69  The  structure  of  the  phase 
Tl2CaBa2Cu20g  (with  Tc  =  112  K)  is  similar  to  that  of  B^CaS^C^Og ,  but 
with  no  evidence  of  a  long-range  modulation.70  71  The  phase  Tl2Ca2Ba2Cu3O10 
has  Tc  =  125  K.  and  has  a  triple  layer  of  CuC>2  planes.71-73 

All  these  materials  are  body-centered  tetragonal  and  have  one,  two,  or 
three  layers  of  Cu02  planes  separated  by  Ca  ions,  sandwiched  by  BaO  layers, 
and  stacked  with  double  T10  layers.66,70  73  The  TIO  layers  are  separated 
by  ^  2.0  A,  significantly  smaller  than  the  distance  between  the  BiO  planes 
in  the  analogous  bismuth  cuprates,  so  that  the  T1  materials  are  more  plate¬ 
like  than  micaceous.66,70-'2  For  Tl2CaBa2Cu20g  and  Tl2Ca2Ba2Cu3Oio,  the 
copper-oxy'gen  planes  are  separated  by  ^  3.2  A,  similar  to  the  multiple-layer 
Bi  cuprates.'0- 13  A  solid  solution  for  T1  and  Ca  has  been  observed  for  these 
materials. 66,70,72  Although  there  is  displacement  of  the  T1  and  O  atoms  from 
their  ideal  positions,  no  long-range  modulation  occurs  in  these  phases.66,70-73 
Short-range  correlations  (less  than  5  A)  in  the  thallium  and  oxygen  disorder 
have  been  observed  in  neutron  diffraction  data.74,75 

Another  series  of  phases  in  the  Tl-Ca-Ba-Cu-0  system  has  only  a  single 
TIO  layer  per  unit  cell.  The  stoichiometries  TlCaBa2Cu207  (with  Tc  =  85  K) 
and  TlCa2Ba2Cu309  (with  Tc  =  110  K)  have  a  primitive  tetragonal  structure 
with  double  or  triple  Cu02  layers  stacked  with  a  TIO  monolayer.76,1 7  The  single- 
layer  Cu02  phase  TlBa2CuOs  does  not  exhibit  superconductivity.'6  The  phase 
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TlCa3Ba2Cu40n  was  discovered  to  contain  four  consecutive  CuC>2  layers78,79 
and  has  Tc  —  122  K. 

Under  certain  fabrication  conditions,  all  of  the  T1  phases  are  observed  to 
have  intergrowths  of  different  numbers  of  CuC>2  layers  and  either  double  or  single 
T10  layers,  resulting  in  a  change  in  the  transition  temperature  compared  with 
the  pure  phase.80  The  materials  with  T10  monolayers  generally  have  increased 
Tc's  with  an  increase  of  intergrowths,  while  the  T10  bilaver  systems  have  lower 
Tc' s  with  more  intergrowths. 

The  increase  in  Tc  with  additional  copper-oxide  layers  observed  for  the 
Bi  and  T1  cuprates  led  to  speculation  that  even  higher  TVs  could  be  found 
by  increasing  the  number  of  adjacent  CuC>2  layers.  However,  the  discovery  of 
Tl2Ca3Ba2Cu40i2  (with  Tc  =  119  K)  and  TlCa4Ba2Cu50i3  (with  Tc  =  110  K) 
shows  a  breakdown  of  the  simple  increase  in  Tc  with  the  number  of  copper-oxide 
planes.78,81  Instead,  a  saturation  (and  decrease)  in  Tc  is  seen  for  additional 
layers  of  Cu02-  Also,  the  material  Cao.86Sr0.i4Cu02  has  been  synthesized.82 
This  tetragonal  structure  is  the  limit  of  a  very  large  number  of  adjacent  CuC>2 
planes.  However,  the  material  is  semiconducting,  even  when  doped  with  alkali 
metals  such  as  Na  or  K,  or  with  trivalent  elements  such  as  Y  or  La,  so  that  Tc 
is  not  simply  dependent  on  adjacent  copper-oxide  planes.82 

A  series  of  high-temperature  superconductors  was  discovered  in  the  system 
Pb2.4Sr2Cu308,  where  .4  =  Y,  La,  Pr,  Nd,  Sm,  Eu,  Gd,  Dy,  Ho,  Tm,  Yb, 
or  Lu.83  Doping  with  Ca  or  excess  Sr  induces  superconductivity  in  these 
materials.  In  particular,  the  phase  Pb2 Yo.sCao.sS^CusOs  has  Tc  =  68  K. 
The  orthorhombic  structure  consists  of  double  CuC>2  planes  separated  by  the 
Ca  or  rare  earth  element,  as  in  YBa2Cu3C>7.  The  copper  perovskite  structure 
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is  sandwiched  by  SrO  layers,  then  stacked  with  a  PbO-Cu-PbO  structure.  The 
copper  layer  separating  the  PbO  planes  is  completely  deficient  in  oxygen,  but 
changes  in  the  oxygen  stoichiometry  are  accommodated  by  inserting  oxygen 
in  this  layer.83  However,  the  oxidation  of  this  copper  layer  is  not  necessary  in 
order  for  superconductivity  to  exist.  The  double  CuC>2  planes  are  separated 
by  ^  3.5  A,  somewhat  larger  than  the  3.2  A  typically  found  in  the  Bi  and  T1 
cuprates.  A  superconducting  phase  with  the  rare  earth  Tb  has  also  been  formed 
in  this  system.84 

The  copper-oxide  superconductors  containing  Bi,  Tl,  and  Pb  all  have 
layers  of  CuC>2  stacked  with  metal-oxide  layers.  There  are  no  one-dimensional 
CuO  chains  in  any  of  these  materials,  structurally  distinguishing  them  from 
YBa2Cu307.17  The  existance  of  high-temperature  superconductivity  without 
the  chains  indicates  that  this  one-dimensional  CuO  structure  is  not  necessary. 
However,  the  common  feature  of  Cu02  planes  within  these  layered  materials  is 
thought  to  play  a  major  role  in  the  specific  mechanism  of  superconductivity. 

The  newest  set  of  cuprate  superconductors  is  closely  related  to  the  La2CuO.} 
system.  The  rare  earth  cuprates  f^CuO-j  with  R  =  Pr,  Nd,  Sm.  Eu,  or  Gd  have  a 
structure  very  similar  to  La2Cu04,  but  with  the  out-of-plane  oxygens  rotated  45° 
to  face-centered  positions.  These  materials  are  not  metallic  or  superconducting, 
even  with  substitution  of  divalent  Sr  on  the  rare  earth  site.85  However,  partial 
substitution  of  Ce  (with  a  formal  valence  of  four)  in  i?2-xCexCu04_y  ( R  = 
Pr.  Nd,  or  Sm)  dopes  the  system  with  electrons  and  allows  superconductivity 
over  24  K.86  In  particular,  the  stoichiometry  Ndi.gsCeo.isCuC^  has  Tc  =  30  K, 
is  strictly  tetragonal,  and  has  undistorted  Cu02  planes.87  Annealing  under  a 
reducing  atmosphere  to  introduce  oxygen  vacancies  ( y  ==  0.07)  enhances  the 
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formation  of  a  bulk  superconductor  in  this  phase.86 

Prior  to  the  discovery  of  superconductivity  in  the  La-Ba-Cu-0  system, 
the  highest  transition  temperature  in  an  oxide  superconductor  was  13  K  for 
BaPbi_yBiy03  at  y  =  0.25. 88  This  material  was  studied  extensively  after  its 
discovery  in  1975,  particularly  since  its  transition  temperature  wras  unusually 
high  compared  with  other  oxide  materials  of  that  time  period.  The  supercon¬ 
ductivity  persists  in  BaPbi_yBiy03  from  9  to  13  K  for  the  solid  solution  range 
0.05  *_]}'_  0.30. 88  The  end  member  BaPbOs  is  metallic.89  while  BaBiC>3  is 
semiconducting.88  All  these  materials  have  the  basic  cubic  perovskite  structure, 
but  with  distortions  that  change  the  symmetry  of  the  unit  cell:90  BaPbOs  is 
orthorhombic,  BaPbo.75Bio.25O3  is  tetragonal,  and  BaBiOs  is  monoclinic.  The 
transition  from  semiconductor  to  superconductor  to  metal  with  Pb  doping  in 
BaPb^yBiyOs  somewhat  parallels  the  behavior  of  Sr-doped  La2_xSrxCu04.8 

Initial  attempts  to  increase  the  trar.aidon  temperature  of  BaPbj  _y  Biy03 
by  addition  of  alkali  metals  such  as  Na  or  K  were  unsuccessful,  but  sharper 
transitions  were  seen  for  these  materials.88,91  However,  high-temperature  super¬ 
conductivity  with  Tc  up  to  22  K  was  found  in  the  Bai_x(K.  Rb)xBi03  system.92 
The  single  phase  Ba0.6Ko.4Bi03  has  Tc  =  30  K, 93,94  and  with  annealing  can 
exhibit  Tc  =  34  K.95  This  material  is  strictly  cubic93,94,96  with  a  =  4.29  A. 
No  ordering  of  Ba  and  K  is  observed  in  the  superconducting  phase,96  and  the 
structure  remains  cubic  in  the  temperature  range  14  -  300  K.97 

The  similarity  in  structure  and  composition  of  BaPbo.75Bio.25O3  and 
Bao.6Ko.4Bi03  suggests  that  these  materials  should  share  the  same  properties  of 
superconductivity.  In  addition,  the  maximum  Tc  in  these  copper-free  materials 
is  at  least  comparable  to  that  of  the  copper  oxides  (particularly  La2 _xSrxCu04 ), 
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indicating  that  the  specific  properties  of  copper  and  its  d  electrons  are  not  re¬ 
quired  for  high-temperature  superconductivity. 

In  summary,  Table  2.1  lists  the  stoichiometry,  superconducting  transition 
temperature,  and  layered  structure  of  the  high-temperature  superconductors. 
The  diversity  of  these  systems  is  remarkable,  yet  many  similarities  are  seen  in 
the  physical  properties.  This  places  strong  constraints  on  theoretical  models 
that  must  describe  all  of  these  systems. 

Superconductivity  and  Transport  Properties 

The  characterization  of  the  high-temperature  superconductors  revealed 
many  similarities  among  the  superconducting  and  normal  state  properties  of 
these  materials,  along  with  several  distinct  differences.  Some  of  the  results 
are  consistent  with  the  standard  Bardeen,  Cooper,  and  Schrieffer  (BCS)  the¬ 
ory  of  superconductivity,98  but  with  notable  exceptions.  The  superconduct¬ 
ing  Tc  of  the  bulk  materials  wras  found  to  increase  slightly  with  pressure 
for  the  materials  La2_x(Ba,  Sr)xCu04,"  YBa2Cu3O7,10°  Bi2CaSr2Cu208,46 
and  Tl2CaBa2Cu208,69  while  BaPbi_yBiy03  has  a  negative  pressure  effect.91 
Hall  coefficient  measurements  in  the  normal  state  indicate  p-type  con¬ 
ductivity  (so  that  the  charge  carriers  are  holes)  in  La2_xSrxCu04, 101,102 
YBa2Cu307,103  Bi2CaSr2Cu208,104  and  Tl2Ca2Ba2Cu3Oio.105  The  Hall  “con¬ 
stant”  Rh  =  —  1/nec  actually  varies  with  temperature  (roughly  as  T_1)  in 
these  materials.101-105  The  bismuth-oxide  superconductors  BaPb!_yBiy03  and 
Bai_xKxBi03  also  have  hole  carriers.106,107  The  trend  of  an  initial  increase  in 
Tr  with  hole  concentration,  n,  is  generally  correlated  to  doping  with  lower  valent 
elements  (such  as  Sr2-  replacing  La3*  and  K1*  replacing  Ba2'1')  or  to  directly 
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TABLE  2.1.  Structures  of  high-temperature  superconductors. 


Stoichiometry 

Tc(  K) 

Metal  stacking  sequence 

Reference 

Laj  ssSro  15CUO4 

38 

L  a-  C  ii”  L  a 

1 

YBa2Cu307 

93 

Ba-Cu-Y-Cu-Ba-Cu 

13 

YBa2Cu408 

81 

Ba-Cu-Y-Cu-Ba-Cu-Cu 

38,39 

Bi2Sr2Cu06 

6 

Bi-Bi-Sr-Cu-Sr 

44,65 

Bi2CaSr2Cu208 

85 

Bi-Bi-Sr-Cu-Ca-Cu-Sr 

45 

TI2  Ba2Cu06 

83 

Tl-Tl-Ba-Cu-Ba 

67 

Tl2  CaBa2Cii208 

112 

Tl-Tl-Ba-Cu-Ca-Cu-Ba 

68 

Tl2Ca2Ba2Cu30io 

125 

Tl-Tl-Ba-Cu-Ca-Cu-Ca-Cu-Ba 

69 

Tl2Ca3Ba2Cu4  0i2 

119 

Tl-Tl-Ba-Cu-Ca-Cu-Ca-Cu-Ca-Cu-Ba 

78 

TlCaBa2Cu20/ 

85 

Tl-Ba-Cu-Ca-Cu-Ba 

76 

TlCa2  Ba2Cu30g 

110 

Tl-Ba-Cu-Ca-Cu-Ca-Cu-Ba 

76 

TICa3Ba2Cu40u 

122 

Tl-Ba-Cu-Ca-Cu-Ca-Cu-Ca-Cu-Ba 

78,79 

TlCa4Ba2CusOi3 

110 

Tl-Ba-Cu-Ca-Cu-Ca-Cu-Ca-Cu-Ca-Cu-Ba 

81 

Pb2^t  0  aCao  sSr2Cu308 

68 

Pb-Cu~Pb-Sr-Cu-(Y,Ca)-Cu-Sr 

83 

-Vdi  SsCeo.lS CUO4 

30 

Nd-Cu-Nd 

86,87 

BaPbo  75B10.25O3 

13 

Ba-(Pb,Bi) 

88 

Bao  6^o  4B1O3 

34 

(Ba,K)-Bi 

92-95 
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increasing  n  through  the  addition  of  more  charge-carrying  planes  (as  in  the  Bi 
and  T1  cuprates).  However,  the  materials  #2-xCexCu04  become  superconduct¬ 
ing  by  doping  Ce4^  on  the  trivalent  rare  earths  f?,  so  that  electron  carriers  are 
added,  as  confirmed  by  the  negative  Hall  coefficient86 

The  anisotropic  crystal  structure  of  the  copper-oxide  superconductors  is 
clearly  demonstrated  in  resistivity  measurements  for  in-plane  and  out-of-plane 
conduction.  The  resistivity  along  the  c  axis  (perpendicular  to  the  Cu02  planes) 
is  generally  several  orders  of  magnitude  larger  than  the  in-plane  resistivity  just 
above  Tc,  with  anistropy  values  of  103  in  La2-xSrxCu04,108  1  00  to  250  in 
YBa2Cu307, 109,110  and  10°  in  B^CaS^C^Og .1U  The  in-plane  resistivity  shows 
a  linear  decrease  with  temperature  down  to  Tc  for  those  samples  with  sharp 
superconducting  transitions,  while  the  resistivity  parallel  to  the  c  axis  indicates 
an  onset  of  localization  at  temperatures  just  above  Tc.  For  cubic  Ba!_xKxBi03, 
resistivity  data  show  hopping  conductivity  at  all  stoichiometries,107  in  contrast 
to  the  metallic  behavior  seen  in  the  cuprates. 

The  pairing  of  charge  carriers  in  the  superconducting  state  was  demon¬ 
strated  in  flux  quantization  measurements.112,113  Flux  vortices  trapped  within 
YBa2Cu307  are  quantized  to  the  value  hc/2e,  consistent  with  2e  pairs.  Voltage 
steps  observed  in  the  ac  Josephson  effect  in  La2_xSrxCu04  and  YBa2Cu307  are 
integral  multiples  of  h/2e  times  the  frequency  of  the  applied  microwave  radia¬ 
tion,  confirming  the  2e  pairing  in  these  materials.114-116  The  pair  wave  function 
appears  to  be  a  singlet  BCS-like  state,  with  no  nodes  observed  in  the  magnetic 
penetration  depth  as  determined  by  muon-spin-relaxation  measurements117  and 
dc  magnetization.118  The  observation  of  the  dc  Josephson  effect  in  tunnel¬ 
ing  between  the  singlet-paired  superconductors  Pb,  Sn,  and  A1  and  the  high- 
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temperature  superconductors  La2-xSrxCuC>4  and  YBa2Cu307  indicates  that 
the  pairs  must  be  of  the  same  symmetry  in  all  these  materials.114-118  Also,  dc 
Josephson  tunneling  has  been  observed  in  all  heterojunctions  of  the  supercon¬ 
ductors  Pb,  \Ba2Cu3O7,  Bi2CaSr2Cu20g,  and  Tl2Ca2Ba2Cu30io,  confirming 
that  all  these  materials  have  singlet  pairing.119 

The  characteristic  energy  gap  A  (at  T  =  0)  of  BCS  superconductivity  has 
proven  difficult  to  determine  in  the  high-temperature  superconductors.  The  BCS 
ratio  2  A  kgTc  =  3.53  represents  the  vveak-coupling  ratio  of  the  energy  gap  to  the 
transition  temperature.98  Initial  measurements  in  tunneling  and  infrared  reflec¬ 
tivity  experiments  gave  a  range  of  values  from  2  to  10  for  this  ratio,  at  least  par¬ 
tially  attributable  to  sample  inhomogeneity.120  Recent  tunneling121  and  infrared 
reflectivity122  measurements  on  YBa2Cu3C>7  yield  a  ratio  between  3  and  4,  con¬ 
sistent  with  weak  coupling.  Tunneling  measurements123  of  BaPbi_yBiy03  give 
a  ratio  of  3.5,  while  reflectivity  measurements124  give  the  ratio  3.2,  again  consis¬ 
tent  with  the  BCS  value.  However,  the  energy  gap  measurements  in  the  other 
materials  are  still  controversial,  with  tunneling,125'126  infrared  reflectivity,127 
and  high-resolution  photoemission128  experiments  yielding  a  ratio  2A/  kBTc  in 
the  strong-coupling  range  4.5  -  8  for  B^CaS^C^Os- 

Another  disputed  measurement  for  these  materials  is  the  change  in  the 
specific  heat,  AC,.,  at  the  transition  temperature,  A  small  change  of  ~  4  percent 
in  Cv  for  YBa2Cu307  has  recently  been  measured.129  Combined  with  a  free- 
electron  estimate  of  the  Sommerfeld  constant  7,  this  gives  a  ratio  <\C\,/~/Tc  near 
the  vveak-coupling  value  of  1.43.98,129 

The  coherence  length,  is  anisotropic  and  quite  short  in  these  materials. 
Estimates  of  £  are  made  from  measurements  of  the  upper  critical  field,  Hc 2,  in 
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these  type  II  superconductors  by  using  the  approximation  Hc 2  5=  hc/2eir£2  for 
one  flux  quantum  per  area  7r£2.  For  example,  in  BijCaS^C^Os  an  in-plane 
coherence  length  of  =  31  A  was  found,  while  along  the  c  axis,  £c  =  4  A.130 
The  coherence  length  in  the  c  direction  is  thus  shorter  than  the  intercell  spacing 
(:=  12  A)  of  the  double  Cu02  layers  in  this  phase.  Interpretation  of  these 
estimates  is  severely  limited  if  Hc 2  is  determined  by  resistivity  measurements, 
since  a  melting  of  the  flux  lattice  is  observed  well  below  Tc  at  all  fields  greater 
than  the  lower  critical  field,  Hci,  resulting  in  finite  flux-flow  resistivity.131,132 
Only  magnetization  measurements  of  Hc 2  may  be  considered  accurate. 

A  particular  property  of  BCS  superconductivity  is  the  isotope  effect,  which 
predicts  that  the  transition  temperature  is  proportional  to  M~a.  where  M  is 
the  mass  of  the  atom  and  a  =  0.5  for  single-element  superconductors.98  Sub¬ 
stitution  of  180  for  160  was  studied  to  determine  the  partial  aox  in  the  multi- 
component  high-tempcrature  superconductors.  Results  for  La2_xSrxCu04  gave 
aox  in  the  range  0.16  -  0.37,  smaller  than  the  BCS  value  but  significant.133,134 
Measurements  for  YBa2Cu3C>7  gave  a  much  smaller  value  aox  <  0.02. 130-138 
A  small  but  nonzero  isotope  shift  has  been  reported  for  180  substitution  in 
multiphase  Bi-Ca-Sr-Cu-0  materials.139  In  the  bismuth-oxide  superconductors 
BaPbj-yBiyOa  and  Ba!_xKxBi03,  the  180  isotope  effect  has  been  reported140 
as  aox  =  0.2  —  0.25,  but  a  value  as  large  as  0.4  for  Bai_xKxBi03  has  been 
suggested.141  Although  in  general  the  isotope  effect  results  from  electron-phonon 
interactions,  the  specific  value  of  a  can  vary  significantly  from  the  BCS  value 
even  in  well-characterized  low- temperature  superconductors,  so  that  small  (or 
large)  isotopic  shifts  in  Tc  do  not  necessarily  limit  the  specific  mechanism  of 
high- temperature  superconductivity. 
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Magnetic  Properties 

Magnetic  ordering  has  been  observed  in  the  non-superconducting  phases 
of  the  copper  oxides.  The  oxygen-deficient  material  La2Cu04_y  exhibits  three- 
dimensional  antiferromagnetic  ordering  on  the  Cu  sites,142-145  with  the  spins  of 
magnitude  ==  0.4  fig  aligned  in  the  copper  planes.  A  canting  of  the  spins  by  0.17° 
from  the  planes  gives  rise  to  a  weak  local  ferromagnetic  moment,  which  vanishes 
at  long  range  since  the  canting  is  opposite  in  alternate  copper  layers.146  The  Neel 
temperature.  Ty,  is  very  sensitive  to  y,  ranging  from  Ty  =  295  K  for  y  =  0.03 
to  r.v  ^  0  K  for  y  =  0.142  At  temperatures  well  above  Ty ,  a  two-dimensional 
quantum  spin-fluid  state  is  observed,  with  instantaneous  correlations  extending 
over  distances  of  200  A.147  No  superconductivity  has  been  observed  in  samples 
exhibiting  antiferromagnetism.  However,  local  magnetic  moments  have  been 
observed  in  superconducting  phases  of  La2_xSrxCu04  (0.07  <  .r  <  0.15)  at 
temperatures  below  3  K,  with  fields  at  least  an  order  of  magnitude  smaller 
than  those  associated  with  long-range  antiferromagnetic  order  in  the  insulating 
state.148,149 

The  parent  compound  YBa2Cu306  is  an  antiferromagnetic  insulator.  The 
Neel  temperature  of  YBa2Cu3  07_y  is  again  sensitive  to  the  oxygen  stoichiom¬ 
etry,  with  a  maximum  Ty  ^  500  K  for  y  =  1.0,  dropping  to  Ty  =  0  K  for 
y  ==  0.6.lo°  An  indication  of  magnetic  ordering  at  temperatures  below  10  K 
was  observed  in  superconducting  phases  of  YBa2Cu307_y ,  possibly  associated 
with  local  antiferromagnetic  order.151  Two  distinct  three-dimensional  antiferro¬ 
magnetic  orderings  are  found  in  YBa2Cu30g.  The  high-temperature  ordering 
is  associated  with  the  copper  plane  sites  only,  while  a  Ty  =  210  K  transition 
occurs  for  ordering  on  the  copper  plane  and  chain  sites.152  A  two-dimensional 
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spin-wave  excitation  is  seen  in  the  non-superconducting  material  above  the  Neel 
temperature.153 

The  bismuth  and  thallium  cuprates  do  not  have  non-superconducting 
parent  compounds  as  do  La2_xSrxCu04  and  YBa2Cu3(>7,  in  the  sense  that 
the  stoichiometric,  undoped  materials  are  superconducting.  However,  these 
materials  can  be  made  non-superconducting  by  various  atomic  substitutions. 
In  particular,  replacement  of  Ca  by  Y  gives  the  phase  Bi2  YSr2Cu20a,  which 
is  insulating. U4~l06  This  phase  exhibits  antiferromagnetic  order  on  the  copper 
sites  with  T \  210  K.107-158  Substitution  of  Y  for  Ca  to  form  the  insulating 

phase  TlYBa2Cu207  also  results  in  three-dimensional  antiferromagnetic  order 
with  r.v  •  350  K.109  Presumably  other  phases  in  the  Bi  and  T1  cuprate  systems 
will  exhibit  this  antiferromagnetic  behavior  with  similar  dopants;  however,  such 
measurements  have  not  yet  been  reported. 

The  rare-earth  cuprates  #2Cu04,  R  =  Pr,  Nd,  Sm,  Eu,  or  Gd.  are 
antiferromagnetic  insulators.  For  example,  the  Neel  temperature  of  Pr2Cu04 
is  270  K.8°  Doping  with  Ce  to  create  the  superconducting  phase  destroys  the 
antiferromagnetic  order.85 

No  magnetic  structure  has  been  seen  in  BaBi03,  BaPbi_yBiy03,  or 
Baj_xKxBi03.160  This  is  as  expected,  since  magnetism  is  not  generally  asso¬ 
ciated  with  valence  s  states  such  as  those  of  Pb  and  Bi.  The  set  of  bismuth 
oxide  high-temperature  superconductors  is  thus  unique  in  its  lack  of  magnetic 
ordering. 
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CHAPTER  III 

TIGHT-BINDING  MODEL 

In  this  chapter  we  present  the  details  of  the  tight-binding  model  used  to 
calculate  the  electronic  structure  of  the  high-temperature  superconductors.  The 
details  of  the  electronic  structure  provide  a  critically  important  foundation 
for  understanding  both  the  normal-state  and  superconducting  properties.  In 
Chapter  IV,  the  electronic  structure  results  of  this  model  are  compared  with 
other  calculations  and  with  experiment.  In  Chapter  V,  this  model  is  used  to 
examine  the  effect  of  atomic  substitutions,  while  vacancy  effects  on  the  electronic 
structure  are  presented  in  Chapter  VI. 

The  validity  of  any  model  for  the  electronic  structure  is  determined  bv  its 
agreement  with  experimental  results.  The  tight-binding  model  has  a  distinct 
advantage  in  this  respect  over  other  techniques,  since  one  can  obtain  a  better 
representation  of  the  electronic  structure  by  using  a  semiempirical  fit  to  the 
experimental  measurements.  Also,  the  parameters  of  the  model  have  a  strong 
chemical  meaning,  representing  effective  atomic  energies  and  bond  strengths. 
The  results  are  interpretable  as  physical  quantities  such  as  valences  and  local 
densities  of  states.  The  tight-binding  approach  does  have  less  predictive  power 
than  techniques  that  use  the  local-density  approximation  (LDA),  in  which  the 
electronic  states  are  calculated  self-consistently  from  first  principles,  but  with 
electron  correlations  and  exchange  effects  approximated  by  a  local  potential.161 
The  LDA  often  provides  a  good  description  of  ground-state  properties,  but  is 
less  accurate  in  representing  the  excited  states.  The  tight-binding  approach  is 
more  useful  for  many  applications  because  of  its  simplicity.162-166 
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Both  tight-binding  and  LDA  models  suffer  qualitative  failures  in  highly- 
correlated  systems.  These  problems  carry  over  to  high-temperature  super¬ 
conductors,  with  La2Cu04  and  YBa2Cu308  predicted  to  be  nonmagnetic 
metals  rather  than  antiferromagnetic  insulators. 142-145,150  Predictions  for  the 
highly-correlated  copper  d  states  may  be  in  disagreement  with  photoemission 
measurements,167  which  indicate  that  the  occupied  d  states  are  more  strongly 
bound  by  1  —  1.5  eV.  Also,  electron-energy-loss  spectroscopy  measurements168 
indicate  that  the  states  nearest  Ep  are  primarily  of  oxygen  p  character,  with  lit¬ 
tle  or  no  Cu  d  admixture.  These  failures  within  the  one-electron  approximation 
are  partially  attributable  to  the  neglect  of  many-body  correlation  effects.169  It 
may  ultimately  be  possible  to  treat  the  electronic  states  accurately  using,  for 
example,  the  variational  quantum  Monte  Carlo  method.1'0 

We  find  in  Chapter  IV'  that  a  single  tight-binding  model,  with  fully 
transferable  parameters,  provides  a  good  description  of  the  electronic  structure 
of  the  high-temperature  superconductors,  in  the  sense  that  the  LDA  results  are 
satisfactorily  reproduced.  The  fitted  tight-binding  parameters  provide  a  basis 
for  extrapolating  to  other  related  classes  of  materials.  Some  features  of  the  one- 
electron  description  appear  to  be  in  at  least  qualitative  agreement  with  photo¬ 
emission,  inverse  photoemission,  electron-energy-loss  spectroscopy,  and  positron 
annihilation  measurements  of  the  high-temperature  superconductors. 

Numerical  calculations  were  performed  in  double  precision  on  a  Digital 
Equipment  Corporation  V'AX  11  '782,  a  VAX  8650,  and  a  VAX  8800.  Matrix 
diagonalizations.  matrix  inversions,  and  numerical  integrations  were  performed 
using  software  from  the  International  Mathematical  and  Statistical  Libraries 
(IMSL).  Detailed  codes  are  available  from  the  author. 
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Electronic  Energy  Bands 

The  energy  bands  of  the  high-temperature  superconductors  are  calculated  in 
the  one-electron  approximation,  using  a  semiempirical  tight-binding  model  with 
valence  p  and  s  orbitals  included  for  the  0  atoms  and  d  and  s  orbitals  for  all  the 
metal  atoms.  Valence  p  states  were  included  for  calculations  with  Al,  Hg,  Tl, 
Pb,  and  Bi.  Atomic  core  levels  were  neglected  in  all  cases,  as  were  the  localized 
/  states.  The  effective  atomic  energy  parameters  es,  ep,  and  ej  were  fitted  to 
the  LDA  energy  bands  of  La2Cu04,  BaBiOs,  and  BaPbC>3,  beginning  with  the 
nominal  values  from  the  table  of  Harrison.165  Energy  parameters  not  given  in 
this  table  were  extrapolated  from  the  nearest  atom  for  which  a  value  is  given. 
The  parameter-  referenced  to  a  zero  of  energy  at  the  ionization  limit.  We 
found  it  nece-  oary  to  shift  the  Cu,  Pb,  and  Bi  energy  parameters  significantly 
from  ‘•heir  tabulated  values  (which  were  originally  applied  to  semiconductors 
and  insulators100,171 ).  We  also  shifted  the  energy  parameters  of  the  neighboring 
atoms  to  maintain  the  tabulated  difference  in  energies  of  the  atomic  states, 
preserving  the  physical  interpretation  of  these  parameters  as  effective  energies. 
The  resulting  atomic  parameters  are  listed  in  Table  3.1. 

The  tight-binding  Hamiltonian  matrix  has  elements 

HQt„3]{k)  =  ^  (aiH'3jJ)  exp[ik  ■  (fjj  -  r,)]  ,  (3.1) 

i 

where  a  and  3  label  orbitals,  and  the  sum  is  over  all  atoms  that  are  equivalent 
to  atom  j  in  the  unit  cell  and  that  neighbor  atom  i.  The  diagonal  elements  of 
the  Hamiltonian  are  the  atomic  energies  in  Table  3.1.  The  interatomic  matrix 
elements  <ai  H  3jJ)  are  determined  by  the  primitive  tight-binding  matrix 
elements  and  by  the  relative  coordinate  from  atom  i  to  atom  jJ.162A6°  The 


TABLE  3.1.  Tight-binding  parameters  for  high-temperature  superconductors. 


(eV) 

*p  (eV) 

td.  (eV) 

u  (A) 

0 

-29.0 

-14.0 

K 

-4.2 

-3.2 

1.2 

Ca 

-5.4 

-3.2 

1.2 

Cu“ 

-12.0 

-14.0 

0.95 

Sr 

-5.0 

-6.8 

1.6 

Y 

-5.5 

-6.8 

1.6 

Ba 

-4.5 

-6.6 

1.6 

La 

-4.9 

-6.6 

1.6 

T1 

-14.8 

-8.3 

-23.0 

1.0 

Pb6 

-18.0 

-9.4 

-29.0 

1.0 

Bi“ 

-21.2 

-10.5 

-35.0 

1.0 

“Parameters  fitted  to  La2Cu04 
6  Parameters  fitted  to  BaPbC>3. 
“Parameters  fitted  to  BaBiC^. 
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n2r3J2 

nidm  —  dfJ2  (l  =  S,p-t  m  =  <7,7t)  , 


where  me  is  the  mass  of  the  electron.  The  matrix  element  varies  with  the 


indicated  inverse  power  of  the  interatomic  separation,  d,  a  result  from  scaling 
in  the  linear  combination  of  atomic  orbitals  (LCAO)  theory.160  The  strength 
of  the  interaction  is  given  by  the  dimensionless  variable  i/wm,  where  l  refers 
to  the  angular  momentum  quantum  number  (5,  p,  or  d )  of  the  first  orbital, 
/'  is  the  angular  momentum  quantum  number  of  the  second  orbital,  and  m 
describes  the  total  angular  momentum  of  the  interaction  ( a  or  n  bonding). 
I'u'm  is  unaffected  under  exchange  of  l  and  /'  if  the  sum  of  the  parities  is  even, 
but  changes  sign  if  the  sum  of  the  parities  is  odd.162  The  strengths  rjS3cr,  etc., 
were  originally  fitted  from  free-electron  theory  to  the  energy  bands  of  silicon 
and  germanium.165171  We  have  adjusted  these  parameters  to  better  fit  the  LDA 
results  for  La2Cu04,  BaBi03,  and  BaPbC>3.  We  obtain  r/,,a  =  —1.1,  r}spi7  =  0.9, 
Hppo  —  l.o,  fjppw  =  0.6,  Psdcr  ■  —1.6,  Tjpd^  =  2.5,  and  T]pdi r  =  1.4.  A  larger 

value  qspcr  =  1.4  is  used  for  the  Bi  and  T1  cuprate  calculations,  but  the  earlier 
calculations  for  La2_xSrxCu04  and  YBa2Cu3C>7  were  insensitive  to  this  value. 
The  extent  of  the  d-orbital  interaction  is  characterized  by  the  length  r^.  The 
parameter  was  fitted  only  for  Cu. 

The  solution  of  the  matrix  form  of  the  time-independent  Schrodinger 


equation 


H(k)  v{k.n)  =  E(k,n)  v{k.n) 


(3.4) 
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gives  tight-binding  wave  functions  t'(k,n)  with  components  representing  the 
atomic  orbital  basis  of  the  model.  The  wave  functions  are  normalized  such  that 

Y  h'a(k,n) j2  =  1  (3.5) 

a 

for  the  orbital  components  a.  The  wave  functions  and  the  energy  eigenvalues 
are  both  functions  of  the  wave  vector  k.  The  energy  eigenvalues  are  plotted  as 
dispersion  relations  in  k  along  particular  symmetry  lines  of  the  Brillouin  zone 
for  each  structure  examined.1'2 

The  results  for  La2_xSrxCu04  and  YBa2Cu3C>7  did  not  require  any  second- 
neighbor  interactions  to  reproduce  the  LDA  energy  bands.  In  particular,  the 
oxygen-oxygen  interactions  were  optimally  zero.  However,  as  described  in 
Chapter  IV,  calculations  for  the  Bi  and  T1  cuprates  required  second-neighbor 
Bi-Bi  and  Tl-Tl  interactions  in  the  rocksalt-like  BiO  and  T10  layers  to  provide 
satisfactory  agreement  with  the  LDA  results.  In  those  cases,  the  same  scaling 
and  strengths  were  used  for  Vwm  as  for  the  nearest-neighbor  interactions,  but 
no  d  —  d  interactions  were  allowed. 

Virtual  Crystal  Approximation 

The  doping  with  Sr  in  La2_xSrxCu04,  with  Pb  in  BaPbj.yBiyOs,  and 
with  K  in  Baj_xKxBi03  requires  a  modification  of  the  energy  parameters  for  the 
atoms  in  the  solid  solution.  These  substitutions  are  treated  in  the  virtual  crystal 
approximation,173  in  which  the  random  distribution  of  dopants  is  assumed  to 
give  an  average  atom  on  a  particular  site.  For  example,  the  doping  La2-xSrx  is 
treated  by  replacing  ea(La)  with  (1  —  .r/2)ej(La)  —  (x/2)e4(Sr),  and  by  removing 
x  electrons  per  chemical  unit.  This  approximation  gives  a  good  description  of 
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dimetal  alloys  in  the  range  where  the  atomic  energies  are  similar,173  and  we 
expect  it  to  hold  in  the  context  of  the  high-temperature  superconductors.  In 
the  particular  case  of  La2_xSrxCu04,  the  energy  parameters  e,  and  for  La 
are  almost  identical  to  the  corresponding  parameters  for  Sr,  so  the  only  effect 
of  doping  with  Sr  is  to  add  hole  carriers  to  the  system. 


Density  of  States 


The  local  density  of  states  for  both  spins  is  calculated  from 


2 


p(E)  =  -  -  Tr  Im  Go{E) 


(3.6) 


where  Tr  indicates  a  trace  over  those  orbitals  associated  with  a  given  site.  The 
Green’s  function  is  given  in  the  spectral  representation  by 

<_'(£,  n.)  u^(k.  n ) 


C'o(E)  =  Y.  “■* 


k,n 


E  -  E(k.n)  -  iS 


(3.7) 


where  E(k.n)  and  v(k.n)  are  the  electronic  energy  eigenvalue  and  the  wave 
function  (in  the  tight-binding  representation)  for  the  nth  band  and  one  of 
the  .V  sample  wave  vectors  k  (with  weight  w ^)  within  the  irreducible  part 
of  the  Brillouin  zone.174  For  these  calculations,  .V  =  64  for  the  orthorhombic 
structures,  .V  =  24  for  the  tetragonal  and  bet  structures,  and  .V  =  20  for  the 
cubic  structures.  A  finite  value  S  —  0.2  eV  was  used,  smoothing  the  results  with 
a  Lorentzian  lineshape.  The  total  density  of  states  is  simply  the  sum  of  the  local 
densities  of  states  of  each  atom  in  the  unit  cell.  The  Fermi  energy  is  calculated 
by  integrating  the  total  density  of  states  up  to  the  number  of  valence  electrons 
for  each  material  (which  is  consistent  with  an  average  of  the  highest  occupied 
energy  state  across  the  sample  wave  vectors).  The  value  of  the  density  of  states 


29 


at  the  Fermi  energy,  p(Ep),  is  compared  with  those  from  other  calculations; 
however,  the  most  reliable  use  for  p(Ep)  in  this  work  is  as  a  measure  of  change 
for  various  substitutions  and  vacancy  effects,  as  detailed  in  Chapters  V  and  VI. 

Atomic  Valence 

The  formal  valence  of  an  atom  (also  known  as  its  oxidation  state)  is  a 
well-defined  quantity  that  provides  a  chemical  prediction  for  the  stoichiometry 
and  bonding  in  compounds.95  For  example,  the  formal  valences  La3r,  Cu2~, 
and  02_  correctly  describe  the  stoichiometry  in  LajCuC^.  In  the  context  of 
this  work,  however,  we  are  interested  in  the  real  charge  that  is  associated  with 
a  particular  atom  bonded  in  the  solid.  If  we  consider  the  solid  to  form  from 
originally  neutral  atoms,  then  charge  transfer  will  occur  between  the  atoms  as 
a  result  of  the  bonding.  The  valence  of  each  atom,  defined  to  be  the  number  of 
electrons  transferred  from  the  atom  when  it  is  bonded  in  the  solid,  is  calculated 
from 

E(k.n)<EF 

An,  =  m  -  2  ^  n),2  ,  (3.8) 

k,a,n 

where  n,  is  the  number  of  valence  electrons  on  the  free  atom  i,  v°(k,n)  is  the 
component  of  the  wave  function  corresponding  to  the  valence  orbital  a  on  this 
atom,  and  the  sum  [over  all  states  with  E{k,n)  <  Ep  is  a  weighted  average  over 
the  N  sample  wave  vectors  k.  The  An  obtained  is  the  real  charge  associated 
with  all  the  valence  basis  states  on  a  particular  atom.  This  technique  can  also 
determine  the  relative  occupancy  of  different  orbitals  on  the  same  atom,  such 
as  the  p j.,  py ,  and  p:  orbitals  of  oxygen. 
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CHAPTER  IV 

ELECTRONIC  STRUCTURE 

In  this  chapter  we  present  the  results  of  the  tight-binding  calculations  for 
the  electronic  structure  of  the  high-temperature  superconductors.  We  have 
calculated  the  electronic  energy  bands,  the  local  and  total  densities  of  states, 
and  the  atomic  valences  for  the  representative  materials  La^ssSro.^CuO.!, 
YBa2Cu307,  BUS^CuOg,  B^CaSrjC^Og,  Tl2Ba2Cu06,  Tl2CaBa2Cu20g, 
Tl2Ca2 Ba2C U3O10 ,  TlCa3Ba2Cu40ni  BaPbg .75Bio.2s03,  and  Bao.eKo.-tBiOs. 
These  results  are  compared  with  other  theoretical  work  and  with  experiments 
that  probe  the  electronic  properties  of  the  high-temperature  superconductors. 
A  summary  of  the  similarities  and  differences  in  these  materials  is  provided  at 
the  end  of  the  chapter. 

Lai.ggSro.isCuC^ 

The  bet  crystal  structure  of  La185Sro.i5Cu04  is  taken  from  experiment,9 
but  the  small  orthorhombic  distortion  and  the  doubling  of  the  unit  cell  are 
ignored  for  simplicity.  The  small  displacements  associated  with  the  tilting  of  the 
CuOg  octahedra  and  the  buckling  of  the  Cu02  planes,0  although  fundamentally 
significant,  would  have  an  imperceptible  effect  in  this  calculation.  The  lattice 
parameters  for  the  bet  cell  are  a  =  3.782  A  and  c  =  13.249  A;  we  use  the 
primitive  unit  cell  with  one  formula  unit  per  cell.  The  partial  substitution  of  Sr 
for  La  is  treated  in  the  virtual  crystal  approximation  as  described  in  Chapter  III. 

The  calculated  energy  bands  for  Laj.gsSro.isCuC^  are  shown  in  Fig.  4.1 
along  several  symmetry  lines  in  the  Brillouin  zone,1'2  with  the  zero  of  energy 
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shifted  to  Ep.  Comparison  with  the  LDA  calculations1 73-1 78  show's  very  similar 
structure.  The  principle  bands  within  10  eV  below  the  Fermi  energy  are  the 
seventeen  bands  contributed  by  Cu  d  and  0  p  states.  The  valence  bandwidth  at 
the  AU  midpoint  (along  the  [100]  direction)  is  7.0  eV,  compared  with  6.5  eV  in 
Ref.  175.  The  bandwidth  at  X  (along  the  110  direction),  the  largest  dispersion 
of  the  bands,  is  9.7  eV,  compared  with  9.2  eV  in  Ref.  178,  9.4  eV  in  Ref.  175, 
9.5  eV  in  Ref.  176.  and  10.3  eV  in  Ref.  177.  The  single  prominant  Cu(d)-0(p) 
antibonding  state  protruding  above  Ep  at  X  has  an  indirect  overlap  with  the 
lowest-lying  conduction  band  at  Z  of  0.9  eV,  compared  with  1.0  eV  in  Refs.  175 
and  176,  and  0.6  eV  in  Ref.  178.  The  relatively  flat  bands  near  5  eV  arise  from 
the  La/ Sr  d  states,  with  the  corresponding  s  bands  about  2  eV  higher  in  energy. 
The  O  s  bands  lie  approximately  18  eV  below  Ep.  The  bands  are  found  to 
have  almost  two-dimensional  character,173-178  with  little  variation  along  TZ  in 
the  kz  direction,  perpendicular  to  the  CuC>2  planes.  There  are  certain  details 
that  are  not  reproduced  in  our  calculation;  for  example,  the  highest  occupied 
spates  along  the  symmetry  line  AU  are  consistently  about  0.5  eV  nearer  Ep  in 
the  LDA  models,  and  some  degeneracies  in  the  non-bonding  oxygen  p  bands 
at  T  and  at  Z  are  broken  in  the  self-consistent  calculations.  However,  the 
level  of  agreement  between  the  tight-binding  bands  and  the  LDA  bands  is 
very  satisfactory,  especially  since  the  range  of  results  in  the  LDA  calculations 
(which  ideally  use  the  same  Hamiltonian)  includes  the  tight-binding  bands.  The 
parameters  of  the  tight-binding  model  are  thus  considered  properly  fitted  for 

Lai  .85Sro.i5Cu04 , 

The  local  densities  of  states  for  Lax.gsSro.isCuC^  are  shown  in  Fig.  4.2, 
with  0(1)  and  0(2)  respectively  labeling  the  in-plane  and  out-of-plane  oxygen 


cell) 
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sites.  Since  0(2)  forms  a  more  ionic  bond  to  Cu,  the  p  state  peak  for  this 
site  is  narrower  than  for  the  0(1)  site.  The  unoccupied  La  d  bands  peak 
near  5  eV  with  little  dispersion.  The  Cu  d  and  0  p  states,  which  have  the 
same  effective  atomic  energy  in  this  model,  are  almost  fully  occupied,  with 
some  Cu-O(l)  hybrid  states  dispersed  above  the  Fermi  energy,  corresponding 
to  the  antibonding  bands  of  Fig.  4.1.  The  total  density  of  states  at  the  Fermi 
energy,  p{Ep ).  is  1.9  states  eV  cell,  compared  with  the  LDA  values  ranging  from 
1.1  states,  eV  cell1'6  to  2.1  states  eV  cell.1'7  The  metallic  character  arises  from 
the  Cu  and  0(1)  sites,  while  the  La  and  0(2)  layers  are  insulating,  in  the  sense 
that  the  local  density  of  states  at  Ep  is  very  small  for  these  atoms. 

The  atomic  valences  for  La^sSro  15Cu04  are  listed  in  Table  4.1.  As  noted 
in  Chapter  III,  this  valence  does  not  equal  the  formal  valence  associated  with 
the  atom,  but  rather  the  real  charge  that  is  distributed  over  the  atomic  basis 
states  of  the  tight-binding  model.  The  valence  of  0(2)  is  significantly  larger  in 
magnitude  than  that  of  0(1),  reflecting  the  charge  transfer  associated  with  the 
ionic  bond  of  the  out-of-plane  0(2)  to  the  La  ion.  The  lower  charge  on  0(1) 
reflects  the  more  covalent  nature  of  the  in-plane  bond  with  Cu. 

Photoemission  measurements  of  La2-xSrxCu04  reveal  the  occupied  elec¬ 
tronic  states.179-181  A  low  density  of  states  at  Ep  is  observed.  The  Cu(d)-0(p) 
hybrid  bands  are  concentrated  within  8  eV  below  the  Fermi  energy,  but  the 
peak  in  the  emission  intensity  is  more  tightly  bound  by  1  —  2  eV  with  respect 
to  the  calculated  density  of  states.180,181  This  discrepancy  occurs  for  both  the 
tight-binding  and  the  LDA  models,1  75  -  1  78  and  may  be  related  to  correlation 
effects  in  the  Cu  d  bands.180 

Inverse  photoemission182-184  experiments  and  electron-energy-loss  spec- 
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TABLE  4.1.  Valences  An  for  Laj.gsSro.isCuC^  and  YBa2Cu307. 


La/Y 

Sr  Ba 

Cu(l)Q 

Cu(2) 

0(1) 

0(2) 

0(3) 

0(4) 

Lai  85Sro.i5Cu04 

2.12 

1.12 

1.08 

-1.14 

-1.51 

YBa2Cu307 

1.77 

1.44 

1.34 

0.98 

-1.12 

-1.09 

-1.11 

-1.21 

,!Cu  site  in  Lai.gsSro.isCuO^ 
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troscopy  (EELS)168,1 ' 9  measurements  probe  the  unoccupied  states  above  Ep 
in  La2_xSrxCu04.  The  vacant  La  d  bands  were  observed  at  5.8  eV  above  the 
Fermi  energy,183  consistent  with  the  tight-binding  result  (==  5  eV)  in  Fig.  4.2. 
The  LDA  calculations  placed  these  bands  much  lower,  within  2  eV  of  Ep.175 
The  La  /  bands  (neglected  in  the  tight-binding  model)  form  a  sharp  peak  near 
8.7  eV.183  The  EELS  data168  indicate  that  only  0  px  and  py  character  is  present 
at  Ep,  with  no  out-of-plane  0  pz  character  observed.  The  extent  of  Cu  d  char¬ 
acter  at  the  Fermi  energy  has  not  been  resolved,  but  the  tight-binding  and  LDA 
results  have  some  C'u(d)-0(p)  admixture  in  the  antibonding  bands  crossing  Ep. 

\Ba2Cu3O7 

The  orthorhombic  crystal  structure  used  for  YBa2Cu307  has  lattice  param¬ 
eters  a  =  3.856  A,  b  —  3.870  A,  and  c  =  11.666  A.18  For  this  calculation,  we 
assume  the  oxygens  are  completely  ordered  in  the  chains  along  the  b  axis.17  No 
adjustment  of  the  tight-binding  parameters  was  made  for  this  calculation,  so  the 
results  test  the  transferability  of  the  parameters  between  materials.  We  choose 
the  symmetry  points  of  the  orthorhombic  Brillouin  zone  as  shown  in  Fig.  4.3, 
with  the  path  TXMY  in  the  central  plane  of  the  Brillouin  zone,  and  ZX'M'Y' 
shifted  along  kz  to  the  zone  boundary.  The  symmetry  line  TY  is  parallel  to  the 
CuO  chains. 

The  energy  bands  for  YBa2Cu307  are  shown  in  Fig.  4.4.  The  complex 
bands  just  below  Ep  arise  from  the  36  states  in  the  Cu(d)-0(p)  manifold. 
The  bands  are  again  two-dimensional,  with  little  dispersion  along  TZ  and 
very  similar  results  along  TXMY  and  ZX'M'Y'.  Three  antibonding  Cu-0 
bands  protrude  above  Ep,  with  one  (associated  with  the  CuO  chain)  having 
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little  dispersion  along  MY,  perpendicular  to  the  CuO  chains.  The  other  two 
bands  are  associated  with  the  Cu02  planes,  and  have  dispersion  similar  to  the 
corresponding  band  in  La185Sr0.15CuO4.  Satisfactory  agreement  is  found  with 
the  LDA  calculations. 185-188  Tor  example,  the  valence  bandwidth  at  M  is  9.7  eV, 
compared  with  7.3  eV  in  Ref.  185  and  8.9  eV  in  Ref.  187.  The  CuO  chain  band 
peaks  1.2  eV  above  Ep  at  Y,  compared  with  1.4  eV  in  Ref.  185,  1.7  eV  in 
Ref.  187,  and  2.1  eV  in  Ref.  188.  For  this  material,  the  range  of  results  is  larger 
than  for  La2_xSrxCu0.j:  however,  the  overall  structural  features  are  in  good 
agreement  among  the  tight-binding  and  LDA  calculations,  indicating  that  the 
tight-binding  parameters  are  transferable  to  YBa2Cu307. 

The  local  densities  of  states  for  the  metal  atoms  in  YBa2Cu3  07  are  shown 
in  Fig.  4.5,  with  those  of  the  oxygen  sites  given  in  Fig.  4.6.  Cu(l)  and  0(1)  refer 
to  sites  in  the  CuO  chain:  Cu(2),  0(2),  and  0(3)  to  the  Cu02  plane  sites:  and 
0(4)  to  the  out-of-plane  oxygen  in  the  BaO  layer.  The  Y  d  and  Ba  d  bands  have 
narrow,  unoccupied  peaks  above  Ep.  0(4)  also  has  a  sharp  spectrum,  related 
to  the  ionic  bond  formed  with  Ba.  The  Cu(2)-0(2.3)  sites  show  broad  features 
just  below  the  Fermi  energy  characteristic  of  the  pda  covalent  bonds  formed  in 
the  Cu02  planes.  The  Cu(l)-0(1)  chain  sites  have  somewhat  narrower  peaks. 
The  total  density  of  states  for  YBa2Cus07  is  3.2  states/eV  cell,  similar  to  the 
LDA  results  of  3.0  states/eV  cell185  and  3.4  states/eV  cell.187  The  metallic 
contribution  to  p{Ep)  is  primarily  from  the  Cu02  plane  sites,  but  with  an 
additional  contribution  from  the  CuO  chains. 

The  valences  An  are  given  in  Table  4.1.  The  valence  for  the  plane  Cu(2) 
site  is  close  to  that  of  C'u  in  Laj  8sSro.i5Cu04,  but  Cu(l)  is  more  ionic.  0(4) 
has  the  largest  valence  of  the  oxygens,  reflecting  the  ionic  bond  formed  with  the 
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Fig.  4.5.  Local  densities  of  states  for  the  metal 
atoms  in  YBa„Cu  0  . 
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Ba  ions. 

The  occupied  electronic  states  of  YBa2Cu307  have  been  determined  in 
photoemission  experiments.16'  189-194  The  peak  of  the  Cu(d)-0(p)  hybrid  bands 
is  more  tightly  bound  by  1  -  1.5  eV  with  respect  to  the  calculated  density 
of  states,  again  attributed  to  localization  effects  in  the  Cu  d  bands.167  The 
oxygen  s  bands  appear  near  —20  eY, 189,191  consistent  with  the  local  density  of 
states  results  in  Fig.  4.6.  Angle-resolved  photoemission  measurements  indicate 
a  dispersion  of  ==  0.3  eY  along  the  TM  symmetry  line  for  the  uppermost 
antibonding  Cu-0  band.194  smaller  than  the  dispersion  predicted  by  the  tight- 
binding  and  LDA  models.  However,  the  angular  dependence  gives  support  to 
the  use  of  band  theory  to  describe  the  electronic  states. 

Inverse  photoemission193-197  and  EELS168,198-201  measurements  reveal  the 
unoccupied  states  in  YBa2Cu307.  The  Ba  d  and  Y  d  bands  were  observed 
at  6.6  eV  and  8.6  eY.19°  respectively,  somewhat  higher  than  the  tight-binding 
results  of  ^  5  eV  and  7  eY  in  Fig.  4.5.  The  Ba  /  bands  are  located  at 
13.4  eY,  well-removed  from  the  region  of  interest  near  the  Fermi  energy.193 
The  EELS  data168  198  show  0  p£  and  py  character  at  Ep  (again  with  no  0 
p.  character),  consistent  with  the  observation  of  oxygen  hole  states  above  Ep 
in  inverse  photoemission  data.197  The  amount  of  Cu  d  character  at  the  Fermi 
energy  has  not  been  determined. 

Another  probe  of  the  electronic  structure  is  positron-annihilation  spec¬ 
troscopy.  Interpretation  of  the  results  is  often  difficult  because  of  trapping  of 
the  positrons  at  vacancies  and  defects  rather  than  in  specific  electronic  bands. 
However,  angular  correlation  studies  of  La2_xSrxCu04  and  YBa2Cu307  indi¬ 
cate  that  the  core  atomic  states  can  be  accurately  studied,202  203  and  are  found 
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to  be  localized  as  in  band  theory. 

B^S^CuOg 

We  have  extended  the  original  tight-binding  model  to  the  superconductors 
containing  thallium,  lead,  and  bismuth.  In  addition  to  the  parameters  previously 
fitted  to  La2Cu04,  the  valence  s,  p ,  and  d  energies  for  Tl,  Pb,  and  Bi,  and  the 
Paper,  Pppcr-  and  r?ppT  interatomic  matrix  element  strengths  are  needed.  (In  the 
calculations  for  La! asSro.isCuC^  and  YBa2Cu3C>7,  no  pper  or  ppir  interactions 
occurred  between  nearest-neighbor  atoms,  and  the  results  were  insensitive  to 
the  value  r}sp(T  =  0.9  used  for  those  materials.)  To  obtain  the  new  parameters, 
we  fitted  the  LDA  energy  bands  of  Mattheiss  and  Hamann204  for  BaBiC>3  and 
BaPbC>3.  The  energy  parameters  for  Tl  were  extrapolated  from  those  of  Pb  and 
Bi.  The  resulting  atomic  parameters  are  listed  in  Table  3.1;  the  interatomic 
strengths  are  q3pcr  =  1.4,  qppi7  =  1.5,  and  qppK  =  -0.6,  as  listed  in  Chapter  III. 

In  order  to  find  better  agreement  between  the  tight-binding  and  LDA  bands 
for  Bi2Sr2Cu06  and  B^CaS^C^Og  (discussed  in  the  next  section),  we  have 
included  the  second-neighbor  Bi— Bi  interactions,  using  the  same  interatomic 
strengths  and  distance  scaling  as  for  nearest  neighbors,  but  neglecting  d-d 
interactions.  (To  be  precise,  interactions  between  the  large  Bi  atoms  are  included 
when  the  interatomic  separation  is  less  than  3.9  A.)  We  regard  these  Bi— Bi 
interactions  as  justified  by  the  large  covalent  radius  of  Bi  and  the  relatively  small 
lattice  parameter  a  =;  3.81  A  in  these  phases,49,65  even  though  only  nearest- 
neighbor  interactions  were  needed  for  Lai.gsSro.isCuC^  and  YBa2Cu3C>7.  We 
ignore  the  long-range  structural  modulation, 48-30,64-66  which  should  have  only 
a  minor  effect  on  the  electronic  properties.  The  displacements  of  the  atoms  in 
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the  BiO  layers  are  also  neglected. 48,33,66  Although  this  approximation  allows 
direct  comparison  between  the  tight-binding  results  and  the  LDA  calculations 
(which  also  used  the  idealized  atomic  positions),  there  may  be  significant  effects 
in  the  electronic  structure  associated  with  the  BiO  disorder,  particularly  for  the 
Bi  p  bands. 

We  use  the  bet  crystal  structure  and  idealized  atomic  positions63  of 
Bi2Sr2Cu06,  with  a  =  3.810  A  and  c  =  24.607  A.  The  energy  bands  are  shown 
in  Fig.  4.7  along  the  symmetry  lines  of  the  bet  cell.1'2  with  D  labeling  the  point 
at  the  Brillouin  zone  boundary  in  the  100  direction.  The  bands  show  two- 
dimensional  structure,  such  as  the  symmetry  of  the  bands  along  TX  and  ZX.  A 
single  Cu(d)-0(p)  antibonding  banu  crosses  Ep  along  TX,  peaking  1.5  eV  above 
the  Fermi  energy  at  X,  compared  with  1.6  eV  in  the  LDA  calculation.205  The  Bi 
p  states  form  dispersive  bands  up  to  =5:  3  eV  above  Ep.  For  this  idealized  struc¬ 
ture,  two  Bi  p  bands  dip  below  the  Fermi  energy  at  D  to  a  minimum  of  —0.7  eV. 
The  LDA  calculation  finds  the  Bi  p  bands  to  cross  below  Ep  to  -0.1  eV.203 
The  disorder  in  the  BiO  layers  may  decrease  the  dispersion  and  prevent  these 
occupied  electron  pockets  from  forming.66,203 

The  local  densities  of  states  for  the  metal  atoms  in  B^SraCuOg  are  shown 
in  Fig.  4.8,  with  those  for  oxygen  in  Fig.  4.9.  Oxygen  site  0(1)  is  in  the  CuOa 
plane,  0(2)  in  the  SrO  layer,  and  0(3)  in  the  BiO  double  layer.  The  Bi  p  bands 
give  a  small  contribution  to  p(Ep)  from  their  dispersion  to  Ep ,  but  the  Sr  d 
bands  are  unoccupied.  The  occupied  Bi  s  bands  peak  near  -10  eV,  while  the 
Bi  d  bands  lie  far  below  Ep  near  -24  eV.  Cu  and  0(1)  display  broad  features 
characteristic  of  the  in-plane  pda  bonds.  0(2)  has  relatively  narrow  p  bands, 
while  0(3)  shows  significant  mixing  with  the  Bi  s  and  p  bands.  This  material 
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TABLE  4.2.  Valences  An  for  Bi-Ca-Sr-Cu-0  superconductors. 


Bi 

Ca 

Sr 

Cu 

0(1) 

0(2) 

0(3) 

Bi2Sr2Cu06 

1.83 

1.15 

0.87 

-1.05 

-1.18 

-1.20 

B^CaS^C^Og 

1.80 

1.51 

1.09 

0.84 

-1.09 

-1.14 

-1.18 
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has  p(Ef)  =  1.6  states  eV  cell. 

Table  4.2  lists  the  valences  for  B^S^CuOe-  The  charges  associated  with 
Cu  and  0(1)  are  smaller  than  those  of  the  Cu02  plane  sites  of  Lai.ssSro.isCuC^ 
and  of  YBa2Cu3  07  in  Table  4.1.  It  is  conceh'able  that  the  lower  carrier  density 
associated  with  the  single  C11O2  layer  accounts  for  the  low  superconducting 
transition  temperature  of  6  K  in  this  phase,  although  the  single  Cu02  layer  of 
Lai.85Sro.i5Cu04  allows  superconductivity  up  to  38  K.  The  BiO  disorder  may 
also  contribute  to  the  low  Tc.66 

Bi2CaSr2Cu20s 

The  bet  crystal  structure  and  idealized  atomic  positions  are  used  for 
B^CaS^C^Os,49  again  neglecting  the  long-range  structural  modulation  and 
the  BiO  disorder.48 The  lattice  parameters  are  a  =  3.814  A  and  c  =  30.52  A 
for  the  bet  cell.49  The  energy  bands  for  B^CaS^C^Os  along  the  symmetry  lines 
of  the  Brillouin  zone  for  the  bet  crystal  structure172  are  shown  in  Fig.  4.10.  These 
bands  are  quite  similar  to  the  results  of  several  LDA  calculations.205-209  Strong 
two-dimensional  character  is  observed,  with  little  dispersion  in  the  direction 
perpendicular  to  the  copper-oxygen  planes.  The  dispersive  Cu(d)-0(p)  bands 
have  a  width  of  about  12  eV  at  X,  somewhat  larger  than  the  LDA  bandwidths 
of  8.7  eY20°  and  9.5  eV.206  Two  pda  antibonding  bands  cross  the  Fermi  energy 
from  below  to  a  peak  of  1.6  eV  at  X,  compared  with  1.7  eV  in  Ref.  205,  2.2  eV 
in  Refs.  207  and  209,  2.4  eV  in  Ref.  206,  and  2.8  eV  in  Ref.  208.  As  in  the 
LDA  calculations,  two  Bi  p  bands  disperse  down  to  the  Fermi  energy,  with  one 
band  crossing  below  Ef  to  -0.7  eV  at  D.  This  forms  occupied  electron  pockets 
near  the  Brillouin  zone  boundary.  However,  this  effect  may  be  destroyed  by  the 
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Fig.  4.10.  Electronic  energy  bands  for  Bi2CaSr2Cua08. 
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disorder  in  the  BiO  layers.30 

Figure  4.11  shows  the  local  densities  of  states  for  Bi,  C a,  Sr,  and  Cu,  with 
the  results  for  oxygen  shown  in  Fig.  4.12.  The  notation  for  the  oxygen  sites 
is  the  same  as  for  Bi2Sr2Cu06  in  the  previous  section.  The  Bi  d  bands  lie 
approximately  24  eV  below  the  Fermi  energy,  while  the  Bi  p  bands,  centered 
slightly  above  Ep ,  interact  strongly  to  produce  a  metallic  contribution  to  p(Ep). 
The  Ca  d  and  Sr  d  states  are  above  Ep,  and  are  quite  ionic.  The  Cu  and 
0(1)  sites  show  strong  interactions,  related  to  the  pda  bonds,  resulting  in  the 
antibonding  bands  protruding  above  Ep  in  Fig.  4.10.  The  0(2)  and  0(3) 
sites  show  primarily  p  character  just  below  the  Fermi  energy  (with  s  bands 
far  below  Ep ),  but  some  p  character  is  distributed  above  the  Fermi  energy 
because  of  mixing  with  the  s  and  p  orbitals  of  neighboring  Bi.  The  density 
of  states  at  the  Fermi  energy  is  2.4  states  eV  cell,  comparable  to  the  value 
p(Ep)  =  2.1  states  eV  cell  from  Ref.  206.  This  gives  a  p(Ejp)  larger  than  in 
Bi2Sr2Cu06-  However,  the  specific  value  of  p(Ep)  is  somewhat  sensitive  to  the 
model,203-209  and  may  be  affected  by  the  BiO  disorder. 

The  atomic  valences  An  are  presented  in  Table  4.2.  The  valences  are  quite 
similar  to  those  for  B^S^CuOg.  The  two  adjacent  Cu02  planes,  separated  by 
the  Ca  ion,  may  support  a  larger  carrier  density  in  Bi2CaSr2Cu20g  than  in 
Bi2Sr2Cu06-  Combined  with  the  larger  p(Ep ),  this  may  account  for  the  higher 
transition  temperature  of  85  K  in  this  phase. 

Photoemission  studies210-216  of  Bi2CaSr2Cu208  confirm  the  basic  occupied 
electronic  states  of  the  tight-binding  and  LDA  calculations205-209  (with  some 
discrepancies  in  detail).  The  spectral  features  near  Ep  are  again  observed  to  be 
more  tightly  bound  by  approximately  1.5  eV  with  respect  to  the  calculated 
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density  of  states.  These  copper-oxygen  hybrid  states  form  occupied  bands 
extending  from  —8  eV  to  a  clear  cutoff  at  Ep. 210-216  The  photoemission  feature 
observed  at  a  binding  energy  of  10  -  11  eV  corresponds  to  the  occupied  Bi  s 
bands  of  the  calculations,210-213  while  the  O  s  spectra  appear  near  -20  eV.210 
Some  Bi  p  character  has  been  observed  at  the  Fermi  energy.213,216 

Angle-resolved  resonant  photoemission  studies217  indicate  an  energy  band 
dispersion  of  0.2  to  O.o  eV  near  the  Fermi  energy  along  the  FX  symmetry  line 
in  Bi2CaSr2Cu208.  The  hole  carriers  are  observed  only  on  the  oxygen  sites 
because  of  strong  Coulomb  repulsion  on  the  copper  sites.217,212  Oxygen  p  states 
just  above  Ep  are  observed  with  resonant  photoemission  from  the  occupied 
oxygen  s  states.21'  These  empty  states  have  only  pt  and  py  symmetry,  so  that 
the  oxygen  holes  are  well-confined  in  the  Cu02  planes.212 

Inverse  photoemission  measurements218  of  Bi2CaSr2Cu208  reveal  unoccu¬ 
pied  states  similar  to  the  tight-binding  density  of  states,  although  the  empty  Ca 
d  and  Sr  d  bands  appear  at  higher  energies  above  Ep  than  calculated.  The  Bi 
p  states  form  a  dispersive  band  centered  =:  4  eV  above  Ep.2is 

Tl2  Ba2  C  uOe 

We  now  consider  the  electronic  structure  of  the  thallium-containing  super¬ 
conductors.  As  noted  above,  the  tight-binding  energies  e,,  ep,  and  td  for  T1 
in  Table  3.1  were  extrapolated  from  the  fitted  parameters  for  Pb  and  Bi.  Be¬ 
cause  of  the  similarity  of  the  crystal  structures  for  the  T1  and  Bi  cuprates,  we 
included  the  same  second-neighbor  Tl-Tl  interactions  that  were  required  for  the 
Bi  atoms  to  properly  model  the  electronic  structures  of  the  Bi  cuprates.  Since 
T1  and  Bi  have  similar  covalent  radii  and  the  same  type  of  valence  orbitals,  we 
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expect  their  bonding  properties  to  be  similar  in  the  tight-binding  model.  None 
of  the  T1  cuprates  appear  to  have  superlattice  modulations  as  do  the  Bi  cuprates. 
However,  there  are  well-characterized  displacements  of  the  T1  and  0  atoms  in 
the  T10  layers,  similar  to  the  BiO  disorder. 66,70-73  We  have  incorporated  these 
displacements  in  the  atomic  substructure  of  each  system.  However,  for  this  cal¬ 
culation  we  neglect  any  partial  substitutions  between  the  T1  and  Ca  sites,66  '0,72 
and  assume  full  occupation  for  these  sites. 

The  bet  crystal  structure  of  Tl2Ba2CuOg  has  the  lattice  parameters  a 
3.866  A  and  c  —  23.239  A.66  The  single  CuC>2  layer  in  the  structure  of  this 
Tc  —  83  K  superconductor  is  strictly  planar.  The  calculated  electronic  energy 
bands  along  the  symmetry  lines  of  the  Brillouin  zone172  are  shown  in  Fig.  4.13, 
with  the  zero  of  energy  shifted  to  Ep.  The  bands  are  virtually  dispersionless 
along  TZ,  perpendicular  to  the  Cu02  plane,  and  are  dominated  as  usual  by  the 
pda  antibonding  band  of  Cu(d)-0(p)  near  the  Fermi  energy.  The  pda  band 
peaks  1.6  eV  above  Ef  at  X.  The  LDA  calculations219,220  show  this  band 
peaking  slightly  higher  in  energy  (2.2  eV  in  Ref.  219  and  2.8  eV  in  Ref.  220). 
The  Tl-Tl  in-plane  interactions  cause  the  T1  p  bands  to  disperse  by  about  3  eV 
from  T  to  D  along  the  100'  direction  in  Fig.  4.13.  These  bands  do  not  dip  belowr 
the  Fermi  energy  as  do  the  Bi  p  bands  in  B^S^CuOg  and  B^CaS^CujOg,  but 
the  Tl(s)-0(p)  hybrid  states  do  cross  below  Ep  to  about  -0.5  eV  at  T.  forming 
occupied  electron  pockets.  A  smaller  dispersion  below  Ep  to  about  -0.1  eV  is 
seen  in  the  LDA  calculation  for  the  T1  s  bands.219 

The  local  densities  of  states  for  the  metal  and  the  oxygen  sites  are  shown 
in  Figs.  4.14  and  4.15.  The  notation  L  be  oxygen  sites  is  similar  to  that  for 
B^S^CuOg:  0(1)  is  in  the  Cu02  plane.  >(2)  is  in  the  BaO  region  separating 
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TABLE  4.3.  Valences  An  for  Tl-Ca-Ba-Cu-O  superconductors. 


T1 

Caa 

Ba 

Cu(l)6 

Cu(2) 

0(1) 

0(2) 

0(3) 

0(4) 

Tl2Ba2CuOg 

1.05 

1.36 

0.89 

-1.10 

-1.04 

-0.72 

Tl2CaBa2Cu208 

1.05 

1.45 

1.37 

0.86 

-1.13 

-1.03 

-0.71 

Tl2Ca2Ba2Cu30io 

1.06 

1.45 

1.38 

0.54 

1.02 

-1.27 

-1.09 

-1.08 

-0.65 

TlCa3Ba2Cu40n 

1.07 

1.45 

1.36 

0.90 

1.07 

-1.12 

-1.06 

-1.10 

-1.17 

a  Averaged  over  inner  Ca(l)  and  outer  Ca(2)  sites  for  TlCa3Ba2Cu40n. 
^Cu  site  for  Tl2Ba2CuOg  and  Tl2CaBa2Cu208. 
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the  Cu  and  T1  layers,  and  0(3)  is  in  the  T10  layer.  The  dispersive  T1  p  bands 
range  from  3  to  6  eV  above  Ep,  while  the  s  states  give  a  metallic  contribution 
to  p(Ep).  Cu  and  0(1)  show  broad  features  below  Ep  arising  from  the  in-plane 
bonding.  In  the  present  model,  the  total  density  of  states  at  the  Fermi  energy 
is  2.0  states /eV  cell. 

The  valences  for  Tl2Ba2Cu06  are  shown  in  Table  4.3.  Notice  that  the 
valences  for  the  Cu  and  0(1)  sites  are  very  close  to  the  corresponding  valences 
for  the  plane  sites  in  Bi2Sr2Cu06  and  Bi2CaSr2Cu20g  shown  in  Table  4.2.  The 
charge  associated  with  T!  is  significantly  smaller  than  that  found  on  the  Bi  sites, 
even  though  both  atoms  are  considered  to  have  a  formal  valence  of  3—  in  these 
materials.  The  Ba  atoms  appear  more  ionic  than  the  corresponding  Sr  atoms  in 
Bi2Sr2Cu06  and  Bi2C'aSr2Cu208. 

Tl2CaBa2Cu208 

The  lattice  parameters  for  the  bet  crystal  structure  of  Tl2CaBa2Cu20g  are 
a  —  3.855  A  and  c  -  29.318  A.70  The  increase  in  the  c  axis  is  consistent  with 
the  inclusion  of  a  CaCu02  layer  into  Tl2Ba2Cu08,  with  a  separation  between 
copper  layers  of  5;  3.2  A.  As  noted  above,  the  displacements  of  the  oxygens  in 
the  T10  layers  from  their  ideal  positions  are  included  in  the  tight-binding  model. 
The  energy  bands  calculated  for  Tl2CaBa2Cu208  are  shown  in  Fig.  4.16.  Two 
degenerate  Cu(d)-0(p)  states  protrude  above  the  Fermi  energy  to  a  peak  of 
1.6  eV  at  X,  corresponding  to  the  two  adjacent  Cu02  planes  in  this  material. 
The  LDA  calculations  show  similar  bands  peaking  2.2  eV221  and  3.0  eV220  above 
Ep.  The  width  of  the  Cu-0  manifold  is  ^  12  eV  at  X,  somewhat  larger  than  the 
bandwidth  of  9  eV  in  the  LDA  results.221  The  unoccupied  T1  p  bands  are  again 
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Fig.  4.16.  Electronic  energy  bands  for  Tl2CaBagCu208. 
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well-removed  from  the  Fermi  energy.  The  Tl(s)-0(p)  hybrid  bands  disperse 
below  the  Fermi  energy  to  -0.3  eV  at  the  symmetry  points  V  and  Z,  compared 
with  an  overlap  of  0.1  to  0.2  eV  in  the  LDA  calculations.220,221 

The  local  densities  of  states  for  Tl,  Ca,  Ba,  and  Cu  are  shown  in  Fig.  4.17, 
with  those  for  oxygen  presented  in  Fig.  4.18.  The  notation  for  the  oxygen 
sites  is  the  same  as  for  Tl2Ba2Cu06  in  the  previous  section.  The  Tl  d  bands 
lie  =:  12  eV  below  the  Fermi  energy,  with  Tl(s)-0(p)  bonding  states  falling 
near  —8  eV.  The  Ca  d  and  Ba  d  states  are  unoccupied  and  insulating,  with  no 
contribution  to  p(Ep).  The  Tl  p  bands  above  Ep  interact  somewhat  weakly  with 
the  neighboring  0(3)  sites.  The  Cu  and  0(1)  sites  show  strong  pdcr  interactions, 
with  the  antibonding  bands  protruding  above  Ep  at  X  in  Fig.  4.16.  The  total 
density  of  states  at  the  Fermi  energy  is  2.7  states /eV  cell.  The  increase  in  p(Ep) 
with  respect  to  that  for  Tl2Ba2Cu06  is  attributed  primarily  to  the  extra  Cu02 
plane  in  Tl2CaBa2Cu208 . 

The  valences  for  Tl2CaBa2Cu208  are  shown  in  Table  4.3,  and  are  very 
similar  to  those  for  Tl2Ba2Cu06-  The  most  significant  factor  contributing  to  the 
higher  Tr  of  112  K  in  Tl2CaBa2Cu208  may  be  th'’  increase  in  p{Ep)  associated 
with  the  double  Cu02  layers. 

Tl2Ca2Ba2Cu3Oio 

The  triple  Cu02  layers  of  Tl2Ca2 Ba2Cu3Oio  increase  the  c  axis  of  the 
bet  crystal  structure  to  35.88  A,  with  a  =  3.850  A.'2  The  increase  in  the 
c  axis  is  again  consistent  with  adding  a  C'aCu02  layer  to  Tl2CaBa2Cu208. 
The  displacements  of  the  oxygens  in  the  T10  layers  are  again  included  in  the 
atomic  substructure  for  the  tight-binding  model.  The  electronic  energy  bands 
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of  Tl2Ca2 Ba2Cu30io  are  shown  in  Fig.  4.19.  This  125  K  superconductor  has 
a  triply-degenerate  antibonding  Cu(d)-0(p)  band  that  crosses  Ep  along  the 
TX  symmetry  line.  This  band  peaks  near  1.7  eV  at  X  in  the  tight-binding 
results,  compared  with  2.5  eV220  and  2.2  eV221  in  the  LDA  calculations.  The 
Cu-0  bandwidth  at  D  (the  point  at  the  Brillouin  zone  boundary  in  the  (100] 
direction)  is  7.9  eV,  somewhat  larger  than  the  LDA  result  of  5.8  eV.221  Small 
electron  pockets  are  again  formed  by  Tl(s)-0(p)  antibonding  states  that  dip 
below  Ep  to  —0.6  eY  at  the  symmetry  points  T  and  Z.  The  LDA  calculation221 
shows  the  occupied  band  dipping  to  —0.3  eV  at  T. 

The  local  densities  of  states  for  the  metal  atoms  of  Tl2Ca2Ba2Cu30io  are 
shown  in  Fig.  4.20,  with  those  for  the  oxygen  atoms  shown  in  Fig.  4.21.  Two 
crystallographical'y  distinct  copper  sites  exist  in  this  structure:  Cu(l)  lies  in 
the  center  of  the  triple  Cu02  layers,  and  is  separated  by  Ca  ions  from  the  outer 
C’u02  planes  with  the  Cu(2)  site.  The  notation  for  the  oxygen  atoms  places 
0(1)  in  the  central  Cu02  plane,  0(2)  in  the  outer  Cu02  planes.  0(31  in  the 
BaO  layer,  and  0(4)  in  the  T10  layer.  Cu(l)  appears  more  ionic  than  Cu(2), 
with  a  narrower  local  density  of  states.  0(3)  and  0(4)  also  have  relatively 
narrow  spectra,  related  to  the  ionic  nature  of  the  bonds  formed  with  Ba  and  Tl. 
T1  p  states  form  a  broad  band  above  Ep,  with  Tl  s  states  contributing  to  p(Ep). 
The  total  density  of  states  at  Ep  for  Tl2Ca2Ba2Cu30io  is  5.3  states  eV  cell, 
considerably  larger  than  p(Ep)  for  the  other  Tl  cuprates,  but  quite  sensitive  to 
small  changes  in  Ep. 

The  atomic  valences  in  Table  4.3  show  that  Cu(2),  in  the  outer  Cu02  planes 
adjacent  to  the  BaO  layers,  has  a  higher  valence  than  that  of  Cu  in  Tl2Ba2CuOs 


and  Tl2CaBa2Cu208 .  However,  the  central  Cu(l)  site  shows  a  decreased  valence 
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Fig.  4.19.  Electronic  energy  bands  for  T1  Ca  Ba  Cu  0 
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of  only  0.54,  while  the  0(1)  site  has  substantially  more  negative  charge  than  the 
0(2)  site  in  the  outer  Cu02  planes.  The  bonds  in  the  central  Cu02  plane 
thus  seem  more  ionic  than  in  the  outer  Cu02  planes.  The  prediction  of  a 
greatly  reduced  valence  for  the  central  Cu(l)  atom  in  Tl2Ca2Ba2Cu30io  may 
be  experimentally  testable. 

Photoemission  and  inverse  photoemission  studies  of  Tl2Ca2Ba2Cu3Oio  have 
revealed  the  occupied  and  vacant  electronic  states.222  The  occupied  copper- 
oxygen  bands  within  ^  8  eV  below  Ep  are  more  tightly  bound  by  about  1.5  eV 
than  the  calculated  results.  T1  may  contribute  some  s  character  at  £f  as  found 
in  the  tight-binding  and  LDA  results.220-221  Unoccupied  bands  are  observed  at 
3.  7,  9.6.  and  14  eV.  associated  with  T1  p.  Ba  d ,  Ca  d,  and  Ba  /  empty  states.222 
These  results  are  in  reasonable  agreement  with  the  calculated  densities  of  states 
(the  Ba  /  states  were  ignored  in  the  tight-binding  model).  Oxygen  p  holes  are 
observed  near  Ep  with  resonant  inverse  photoemission.222 

TlCa3  Ba2Cu40n 

We  next  examine  TlCa3Ba2Cu40n  as  a  representative  member  of  the  T1 
cuprates  that  contain  only  a  single  TIO  layer  per  unit  cell.  This  particular  phase 
has  Tc  -  122  K.  and  has  a  simple  tetragonal  structure  with  a  =  3.85  A  and 
c  —  19.01  A.78’*9  The  atomic  positions  for  TlCa3Ba2Cu40n  were  extrapolated 
from  those  of  the  analogous  single  TIO  layer  superconductor  TlCa2Ba2Cu309 
by  inserting  an  additional  CaCu02  layer.76  Lacking  other  experimental  data, 
the  idealized  atomic  positions  for  the  TIO  layer  have  been  used  in  this  case. 
The  electronic  energy  bands  of  TlC'a3 Ba2Cu40i i  are  shown  in  Fig.  4.22  along 
the  symmetry  lines  of  the  simple  tetragonal  Brillouin  zone.1'2  The  point  X  lies 
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along  the  100  axis,  and  M  lies  at  the  Brillouin  zone  boundary  along  the  110; 
direction.  The  crystal  structure  contains  four  CuC>2  layers,  separated  by  Ca 
ions,  which  contribute  the  Cu(d)-0(p)  antibonding  bands  that  form  the  hole 
conduction  bands.  These  four  degenerate  bands  peak  1.8  eV  above  Ep  at  \1 
and  R.  The  Cu-0  bandwidth  at  \I  is  12.2  eV.  These  bands  are  essentially 
two  dimensional,  as  demonstrated  by  the  similarity  of  the  dispersion  along  the 
symmetry  lines  TM  and  ZR.  The  TIO  layer  provides  a  single  Tl(s)-0(p)  band 
that  just  crosses  below  Ep  at  T  and  Z. 

The  local  densities  of  states  for  the  metal  atoms  are  shown  in  Fig.  4.23. 
with  those  for  oxygen  shown  in  i  ig.  4.24.  Cu(l)  labels  the  copper  sites  in  the 
two  inner  Cu02  planes,  while  Cu(2)  labels  the  copper  sites  in  the  outer  planes. 
Ca(l)  lies  in  the  center  of  the  four-layer  structure,  while  Ca(2)  labels  the  two 
symmetric  sites  adjacent  to  the  outer  Cu02  planes.  The  oxygen  sites  are  labeled 
just  as  for  Tl2Ca2Ba2Cu30iu  in  the  previous  section.  The  two  inner  copper- 
oxide  planes  are  relatively  isolated  from  the  remaining  structure  by  the  two 
outer  Cu(2)-0(2)  layers,  but  the  Cu(l)  and  Cu(2)  spectra  are  similar.  The 
C'a(l)  and  C'a(2)  unoccupied  d  bands  are  almost  identical.  The  0(3)  and  0(4) 
sites  have  quite  narrow  bands,  indicating  the  ionic  nature  of  their  bonds  with 
Ba  and  Tl.  The  total  density  of  states  at  Ep  within  the  tight-binding  model  is 
7.5  states/ eV  cell,  the  largest  value  among  the  systems  examined. 

The  results  of  the  valence  calculation  for  TlCa3Ba2Cu.}Oi1  are  listed  in 
Table  4.3.  The  outer  Cu(2)  sites  have  the  largest  valence  of  all  the  Tl  cuprate 
superconductors.  The  0(1)  site  again  has  slightly  more  negative  charge  than 
the  oxygens  in  the  outer  Cu02  planes.  0(4)  is  significantly  more  ionic  than  the 
oxygens  in  j  materials  with  double  TIO  layers. 


Fig.  4.23.  Local 
atoms  in  TlCa.I 
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BaPbo.75  B10.25O3 

Let  us  now  examine  the  bismuth-oxide  superconductors.  The  material 
BaPbj-yBiyOa  has  Tc  —  13  K  for  y  =  0.25. 88  This  phase  was  considered  to 
have  an  unusually  high  Tc  for  an  oxide  prior  to  the  discovery  of  the  copper-oxide 
superconductors.  This  system  is  much  simpler  to  examine  than  the  cuprates,  and 
is  very  similar  to  the  high-temperature  superconductor  Bai_xKxBi03  (described 
in  the  next  section).  The  tight-binding  energy  parameters  used  for  Bi  and  Pb  are 
the  same  as  those  fitted  to  the  LDA  energy  bands204  of  BaBiOs  and  BaPbC>3. 
The  doping  is  treated  in  the  virtual  crystal  approximation. 

The  tight-binding  energy  bands  in  Fig.  4.25  for  BaPbo.75Bio.25O3  agree  well 
with  the  previous  LDA  calculation223  for  BaPbo.7Bio.3O3,  with  both  calculations 
using  a  model  cubic  structure  with  a0  =  4.29  A  rather  than  the  actual  tetragonal 
structure.90  The  ten-band  complex  near  the  Fermi  energy  is  composed  of 
Pb  Bi(s)-0(p)  hybrid  states,  with  the  highest  antibonding  band  protruding 
above  Ef  to  3.0  eV  at  R.  as  in  the  LDA  result.223  The  tight-binding  bandwidth 
at  R  is  16  e\\  slightly  larger  than  the  width  of  14  eV  in  the  LDA  calculation. 
The  density  of  states  at  the  Fermi  energy  is  0.6  states,  eV  cell  for  the  mode! 
structure,  contributed  by  the  Pb/Bi  s  and  O  p  states.  The  LDA  result  is 
0.3  states  eV  cell.223  The  calculated  valences  for  BaPbo.75Bio.25O3  are  listed 
in  Table  4.4;  the  effective  valence  of  the  Pb  Bi  atom  is  1.54. 

A  photoemission  study  of  BaPb0  75  Bio.2503  shows  that  the  valence  bands 
down  to  -8  eV  are  derived  from  Pb  Bi(.s)-0(p)  hybrid  states,  with  a  contribu¬ 
tion  between  -8  and  -12  eV  from  the  Pb  and  Bi  s  states.224  The  core-like  Ba 
p  states  appear  near  -14  eV,  with  no  Ba  d  or  s  states  appearing  at  Ef. 
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Fig.  4  .25.  Electronic  energy  bands  for  BaPb076Bi026O3. 
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TABLE  4.4.  Valences  An  for  BaPbo.75Bio.25O3  and  Ba0.6Ko,4Bi03. 


Ba 

K 

Bi 

Pb 

0 

Bapbo.75Bio.2503 

1.49 

2.29 

1.29 

-1.01 

Bao,6  K0.4  Bi03 

1.62 

0.62 

1.80 

-1.01 
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Bao.6Ko.4Bi03 

The  high-temperature  superconductor  Bai_xKxBi03  has  Tc  ==  34  K  at  x  — 
0.4. 93-93  The  key  to  the  higher  Tc  in  this  phase  compared  to  BaPbo.75Bio.25O3 
seems  to  be  that  the  hole  doping  is  performed  on  the  Ba  site  rather  than  on  the 
Bi  site,  so  that  the  conduction  bands  are  relatively  undisturbed,  fhe  structure 
used  for  Ba0.6K0.4BiO3  is  strictly  cubic,93  with  a0  =  4.293  A  and  an  ideal 
perovskite  atomic  substructure.  The  doping  is  again  treated  in  the  virtual  crystal 
approximation. 

The  energy  bands  in  Fig.  4.26  are  in  good  agreement  with  the  recent 
LDA  calculation223  for  Bao.sKo.sBiOs,  and  are  quite  similar  to  the  bands  for 
BaPbo.75 Bi0  25O3  in  Fig.  4.25.  The  conduction  band  is  formed  from  Bi(s)-0(p) 
antibonding  states.  The  bandwidths  at  M  and  R  are  12.9  and  17.0  eV,  compared 
with  12.7  and  15.2  eV  in  the  LDA  results.225  The  larger  bandwidth  at  R  in  the 
tight-binding  calculation  is  partially  attributed  to  the  omission  of  the  core-like 
Ba  p  states  near  —12  eV.223  The  unoccupied  Bi  p  bands  disperse  from  2  to  5  e V 
above  Ep .  while  the  empty  Ba  d  bands  form  a  flat  band  at  6.3  eV. 

The  local  densities  of  states  for  Bao.6  K0.4  BiOs  are  shown  in  Fig.  4.27.  The 
Bi  s  and  O  p  states  form  broad  spectra  from  -10  eV  to  the  Fermi  energy,  and 
contribute  the  metallic  p(Ep)  =  0.7  states  eV  cell.  The  Ba  K  d  bands  are  ionic 
and  insulating. 

The  valences  for  Bao.6Ko.4Bi03  are  listed  in  Table  4.4,  with  an  effective 
Ba  K  valence  of  1.22.  We  find  an  increase  of  0.26  holes  per  Bi02  unit  in 
Ba0  6K0.4BiO3  with  respect  to  the  Pb  Bi  oxide  region  of  BaPbo.75Bio.25O3-  The 
calculation  also  indicates  that  doping  Ba5_xKxBi03  at  .r  =  0.4  gives  a  carrier 
density  of  0.22  holes  per  Bi02  unit  with  respect  to  semiconducting  BaBiOs,  while 
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Fig.  4.27.  Local  densities  of  states  for  Ba„„K„,BiO„. 
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the  remaining  0.18  doped  holes  are  localized,  lowering  the  charge  associated  with 
the  Ba/K  site. 

Photoemission  studies226,227  of  Bao.eKo^BiOa  show  occupied  states  similar 
to  the  calculated  bands.  Oxygen  p  states  dominate  the  structure  near  the  Fermi 
energy,  with  0  3  core  levels  seen  near  —21  eV.  Bi(s)-0(p)  hybrid  states  are 
seen  down  to  —13  eV.  Inverse  photoemission  results227  place  the  unoccupied  Bi 
p  bands  near  4  eV  and  the  Ba  d  bands  near  7  eV,  in  agreement  with  the  local 
densities  of  states  for  Ba0.6K0.4BiO3  in  Fig.  4.27. 

Summary  of  Electronic  Structure 

We  find  that  a  single  tight-binding  model,  with  fully  transferable  pa¬ 
rameters,  provides  a  good  description  of  the  electronic  structure  of  the  high- 
temperature  superconductors,  in  the  sense  that  the  results  satisfactorily  agree 
with  the  LDA  calculations  and  with  certain  experimental  details  of  the  electronic 
states.  The  parameters  fitted  to  I^CuO.*,  BaPbC>3,  and  BaBi03  provide  suffi¬ 
cient  input  to  extrapolate  other  parameters  needed  for  the  model. 

All  the  copper-oxide  superconductors  have  two-dimensional  conduction 
bands  arising  from  Cu(d)-0(p)  antibonding  states.  The  metal-oxide  layers 
sandwiching  the  CuC>2  planes  (such  as  LaO,  BaO,  and  SrO)  all  display  ionic 
and  insulating  behavior,  in  the  sense  that  there  is  large  charge  transfer  and  the 
contribution  to  p(Ep)  is  small.  The  copper-free  bismuth-oxide  superconductors 
BaPbo.75Bio.25O3  and  Bao.eKo^BiOs  have  valence  bands  of  antibonding  Bi(s)- 
O (p)  states,  with  cubic  or  almost-cubic  band  symmetry.  These  materials  also 
have  ionic  BaO  in  their  electronic  structures. 
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CHAPTER  V 

ATOMIC  SUBSTITUTION  EFFECTS 

In  this  chapter  we  report  studies  of  the  effect  of  various  atomic  substitutions 
on  the  electronic  properties  of  the  high-temperature  superconductors.  For  the 
material  YBa2Cu3  07,  we  examined  the  following  atomic  substitutions:  Al,  Fe. 
Co.  NT  or  Zn  replacing  Cu;  Sr  or  La  replacing  Ba;  and  F  or  N  replacing  0.  The 
phase  Bi2CaSr2Cu208  was  investigated  for  the  substitution  of  Pb  on  the  Bi  site. 
Of  the  many  phases  of  T1  cuprates,  the  representative  material  Tl2CaBa2Cu208 
was  examined  for  replacement  of  T1  with  Hg  or  Pb.  In  each  case,  the  change  in 
the  densities  of  states,  the  modified  atomic  valences,  and  the  shift  in  the  Fermi 
energy  were  calculated. 

There  have  been  many  experimental  studies  of  changes  in  the  supercon¬ 
ducting  properties  of  YBa2Cu307  induced  by  replacement  of  the  original  atomic 
species.  The  replacement  of  Y  by  any  of  the  rare  earth  elements  Xd,  Sm,  Eu. 
Gd.  Dy,  Ho,  Er,  Tm,  Yb,  or  Lu  had  only  a  very  small  effect  on  J^26-228-233 
indicating  a  remarkable  insensitivity  of  the  relevant  electronic  and  structural 
properties  to  the  species  on  this  site.  (The  rare  earths  Ce,  Pr,  and  Tb  do 
not  form  superconducting  phases  in  f?Ba2Cu307,  perhaps  because  of  the  stable 
tetravalent  state  of  these  atoms. 26'95'234-236  )  In  particular,  the  magnetic  mo¬ 
ment  of  the  substitutional  atom  does  not  appear  to  affect  the  superconducting 
transition  temperature,232'233  indicating  that  the  Y  site  is  electronically  isolated 
from  the  superconducting  region  of  this  material. 

Replacement  of  Ba  by  Sr  was  found  to  yield  a  depression  of  Tc  that  is 
linear  in  the  Sr  concentration,  and  which  was  attributed  to  lattice  distortions 
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in  the  neighborhood  of  the  Sr  atom.237,238  Substitution  of  La  onto  the  Ba  site 
leads  to  a  more  rapid  decrease  of  Tc ,  with  a  sudden  loss  of  superconductivity 
in  YBa2-xLaxCu307  for  x  >  0.4. 33,239,240  This  effect  was  attributed  to  charge 
compensation  for  the  donated  electron  of  La  by  a  reduction  of  the  copper  valence 
or  by  an  increase  in  the  oxygen  content.240  Further  studies  of  the  substitution 
of  the  rare  earths  Nd,  Sm,  Eu,  Gd,  and  Dy  for  both  Y  and  Ba  indicate  a 
strong  suppression  of  Tr  in  /?Ba2_xf?xCu307,  with  a  loss  of  superconductivity 
for  x  >  0.5. 241 

Substitutions  for  copper  lead  to  a  particularly  dramatic  lowering  of  Tc.  The 
substitution  of  Ai  onto  the  Cu(l)  chain  site  in  single  crystal  YBa2Cu3_x Alx07 
leads  to  a  rapid  decrease  in  Tc  for  x  >  0.1,  with  complete  loss  of  superconduc¬ 
tivity  for  x  >  0.22. 242  For  the  3d  transition  metals,  substitution  in  polycrys¬ 
talline  YBa2Cu3_x.l/x07  for  x  £  0.3  yields  a  strong  depression  of  the  transition 
temperature.243-247  Note  particularly  that  nonmagnetic  Zn  suppresses  Tc  more 
than  magnetic  Fe  or  Co.  Replacement  of  Cu  by  Ag  also  leads  to  a  decreased 
transition  temperature.248,249 

The  relative  importance  of  the  Cu(l)  and  Cu(2)  .  :tes  for  high-temperature 
superconductivity  in  YBa2Cu307  has  been  examined  by  correlating  the  specific 
site  for  various  metal  substitutions  with  the  effect  on  the  superconducting  prop¬ 
erties.  Xeutron  diffraction  measurements  place  Fe  substitutions  on  the  Cu(l) 
site,250-252  although  other  experiments  indicate  that  Fe  may  also  substitute 
on  the  Cu(2)  site.253  Doping  with  Co  shows  preferential  substitution  on  the 
Cu(l)  site,251,234  along  with  an  increase  of  the  oxygen  stoichiometry  above  7.0. 
Both  Fe  and  Co  substitution  allow  superconductivity  to  persist  (at  a  reduced 
Tc)  up  to  relatively  large  dopant  levels  (x  0.4  in  YBa2Cu3_x.Ux07),  indicat- 
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ing  that  the  Cu(l)  chain  site  does  not  provide  the  primary  states  to  support 
superconductivity.250,254  However,  Ni  and  Zn  substitutions  are  found  to  occupy 
the  Cu(2)  sites  in  the  copper-oxide  planes.251,255,256  The  sharp  depression  in  Tc 
with  Zn  substitution  is  thus  directly  attributable  to  disruption  of  the  important 
Cu02  plane  region. 

Replacement  of  oxygen  with  sulfur  in  YBa2Cu306S  appears  to  enhance  the 
Meissner  effect,  while  leaving  Tc  unaffected  near  90  K.257  The  substitutions  of 
fluorine238  and  nitrogen239  into  YBa2Cu307.  presumably  onto  the  oxygen  sites, 
have  produced  conflicting  and  often  irreproducible  results.260-264 

Substitution  of  Pb  into  Bi2CaSr2Cu20a  has  been  found  to  increase  the 
onset  of  superconductivity  to  107  ft  30-60-62  This  effect  has  been  attTbuted 
to  stabilization  of  the  higher-rc  structure  Bi2Ca2Sr2Cu3Oio , 30,61,6  out  the  p- 
type  dopant  nature  of  Pb  replacing  Bi  may  also  play  a  role  in  the  enhanced 
superconductivity. 

The  T1  cuprates  display  high- temperature  superconductivity  when  Ba 
is  completely  replaced  by  Sr.26o-26‘  The  phase  with  nominal  composition 
TI0.75 Bio.2sCaSr2Cu207  displays  superconductivity  at  Tc  =  75  K.263  A  slightly 
higher  Tc  near  85  K  is  obtained  in  the  Pb-doped  phase  Tl0.5Pb0.5CaSr2Cu2O7, 
while  superconductivity  at  Tc  120  K  is  observed  for  the  material 

Tlo.5Pbo.5Ca2Sr2Cu309.267  The  phases  with  complete  substitution  of  Sr  for  Ba 
require  additional  doping  on  the  T1  site  with  Pb  or  Bi  to  stabilize  the  structure. 

Method 

I'sing  the  same  tight-binding  model  for  the  electronic  structure  described 
in  Chapter  III,  we  consider  the  effects  of  atomic  substitutions.  For  simplicity,  we 
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employ  a  periodic  substitution  scheme,  with  the  substitutional  atom  occupying 
the  same  location  within  each  unit  cell  throughout  the  material.  Systems  that 
experimentally  accommodate  full  atomic  substitutions  can  be  directly  compared 
with  these  results.  This  technique  does  not  address  partial  substitutions  or 
defect  ordering  as  a  function  of  dopant  fraction;  however,  the  changes  in  the 
electronic  properties  calculated  for  the  full  atomic  substitution  provide  an 
interpolative  description  of  the  model  system  at  all  dopant  levels. 

Starting  with  the  fitted  parameters  in  Table  2.1.  the  atomic  energies  for 
substitutional  atoms  are  extrapolated  from  the  solid  state  table  of  Ref.  165. 
maintaining  the  differences  in  atomic  energies  between  elements.  We  extrapolate 
the  energies  e3,  ep,  and  for  T1  and  Hg  from  the  fitted  Bi  and  Pb  atomic 
energies.  The  resulting  parameters  are  listed  in  Table  5.1.  We  neglect  the  small 
changes  that  occur  in  the  lattice  constants  due  to  differences  in  the  covalent  radii 
of  the  substitutional  atoms.  Although  such  changes  in  the  lattice  parameters 
are  observed,  and  may  have  a  role  in  the  superconductivity  of  these  materials, 
they  nevertheless  have  only  a  minor  effect  within  the  tight-binding  model  for 
the  electronic  structure  of  the  material. 

In  the  calculations  for  BiPbCaS^C^Og  and  TLM  CaBa2Cu2  0s  (M  =  Hg 
or  Pb),  we  include  the  same  second-neighbor  interactions  (e.g.,  Pb-Pb)  that 
were  required  to  properly  describe  the  electronic  energy  bands  of  B^S^CuOe 
and  B^CaS^C^Og  in  Chapter  IV.  Since  the  covalent  radii  of  Hg,  Tl,  Pb, 
and  Bi  are  all  nearly  equal,  we  expect  these  atoms  to  have  similar  bonding 
characteristics  within  the  same  crystal  structure. 

The  notation  in  labeling  the  substitutional  atoms  parallels  the  form  for 
the  undoped  materials  as  given  in  Chapter  IV.  For  example,  A 1  ( 1 )  refers  to 
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TABLE  5.1.  Tight-binding  parameters  for  atomic  substitution  calculations. 


(eV) 

eP  (eV) 

^  (eV) 

rd  (A) 

N 

-23.0 

-11.5 

0 

-29.0 

-14.0 

F 

-36.0 

-17.0 

A1 

-10.1 

-4.9 

Fe 

-11.6 

-10.4 

0.95 

Co 

-11.7 

-11.6 

0.95 

Xi 

-11.9 

-12.8 

0.95 

Cu 

-12.0 

-14.0 

0.95 

Zn 

-13.5 

-15.2 

0.95 

Sr 

-5.0 

-6.8 

1.6 

Ba 

-4.5 

-6.6 

1.6 

La 

-4.9 

-6.6 

1.6 

Hg 

-11.6 

-7.2 

-17.0 

1.0 

T1 

-14.8 

-8.3 

-23.0 

1.0 

Pb 

-18.0 

-9.4 

-29.0 

1.0 

Bi 

-21.2 

-10.5 

-35.0 

1.0 
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an  aluminum  atom  that  replaces  the  copper  atom  on  the  chain  site  Cu(l)  in 
YBa2Cu307,  while  Zn(2)  refers  to  a  zinc  atom  that  replaces  a  plane  site  Cu(2). 

YBa2Cu24/07,  M  =  Al,  Fe,  Co,  Ni,  or  Zn 

The  substitution  of  Al  onto  the  Cu(l)  site  has  a  dramatic  effect  on  the 
band  structure  of  YBa2Cu2A107,  shown  in  Fig.  5.1,  as  compared  with  that  of 
\Ba2Cu307  in  Fig.  4.4.  In  particular,  the  conduction  state  associated  with 
the  Cu(d)-0(p)  antibonding  band  of  the  one-dimensional  chain  is  completely 
removed  from  the  structure,  indicating  a  complete  loss  of  metallic  conductivity 
for  this  region.  However,  the  two-dimensional  nature  of  the  electronic  energy 
bands  is  maintained,  as  demonstrated  in  Fig.  5.1  by  the  similarity  of  the  two 
wave-vector  paths  TXMY  and  ZX'M'Y'.  which  trace  the  edge  of  the  Brillouin 
zone  in  the  F  and  Z  planes.  The  total  density  of  states  of  the  Al-doped  material  is 
shown  in  Fig.  5.2.  Note  that  the  Al  p  bands  lie  about  6  eV  above  Ep ,  leading  to 
weak  interactions  with  the  neighboring  chain  oxygens.  The  calculated  valences 
in  Table  5.2  show  only  minor  changes  in  the  Cu02  plane  regions,  but  A 1  ( 1 ) 
appears  much  more  ionic  than  Cu(l)  in  YBa2Cu3C>7.  Al(l)  substitution  thus 
leads  to  a  localized  destruction  of  conduction  states  within  the  one-dimensional 
chains,  and  substantially  increases  the  ionic  character  of  the  chain  region.  The 
shift  in  the  Fermi  energy  due  to  the  A 1  ( 1 )  substitution,  along  with  the  total 
density  of  states  at  Ep,  is  given  in  Table  5.3.  The  Fermi  energy  shifts  upwards, 
and  the  removal  of  the  Cu(l)-0(1)  chain  band  results  in  a  decrease  in  p(Ep). 

For  the  transition-metal  substitutions,  we  consider  both  Cu(l)  and  Cu(2) 
sites  within  YBa2Cu307,  and  perform  the  calculation  independently  for  each 
substitution  site.  The  total  densities  of  states  for  the  replacement  of  Cu  with 
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Fig.  5.2.  Total  density  of  stales  for  YBa2Cu2A107. 


TABLE  5.2.  Valences  An  for  YBajC^A/Oj. 


Cu(l) 

Cu(2) 

M 

0(1) 

0(2)a 

0(3)“ 

0(4) 

Al(l)" 

0.94 

2.46 

-1.59 

-1.10 

-1.11 

-1.50 

Fe(l)‘ 

0.38 

3.24 

-1.13 

-1.30 

-1.31 

-1.19 

Fe(2) 

0.67 

0.36 

2.85 

-1.35 

-1.15 

-1.16 

-1.30 

Co(l)b 

0.67 

2.03 

-1.10 

-1.17 

-1.17 

-1.14 

Co(2) 

1.03 

0.76 

1.44 

-1.18 

-1.04 

-1.06 

-1.25 

Ni(l) 

0.99 

1.12 

-1.12 

-1.07 

-1.07 

-1.18 

Xi(2)k 

1.31 

0.96 

0.82 

-1.12 

-1.04 

-1.05 

-1.23 

Cu 

1.34 

0.98 

-1.12 

-1.09 

-1.11 

-1.21 

Zn(  1) 

1.02 

1.41 

-1.33 

-1.07 

-1.07 

-1.25 

Zn(2)f> 

1.37 

0.99 

1.20 

-1.11 

-1.15 

-1.17 

-1.23 

'‘Averaged  over  sites  in  unit  cell. 

6  Experimentally-determined  substitution  site. 
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TABLE  5.3.  Shift  in  Ep  and  density  of  states  for  YBa2Cu2d/07. 


M 

A  Ef 

(eV) 

p(Ef) 

(states/e V  cell) 

Al(l)a 

0.15 

2.4 

Fe(l)a 

0.79 

6.3 

Fe(2) 

0.78 

9.3 

Co(l)“ 

0.22 

7.5 

Co(2) 

0.16 

6.3 

Xi(l) 

-0.12 

4.2 

Ni(2)a 

-0.08 

3.5 

Cu 

3.2 

Zn(l ) 

0.16 

3.0 

Zn(2)a 

0.08 

3.3 

“Experimentally-determined  substitution  site. 
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Fe,  Co,  Ni,  or  Zn  are  shown  in  Figs.  5. 3-5. 6.  The  results  for  Fe  substitution  in 
Fig.  5.3  show  a  large  contribution  from  the  Fe  d  bands  at  Ep,  while  the  results  for 
Co  in  Fig.  5.4  have  Co  d  character  at  and  below  the  Fermi  energy.  The  Ni  d  and 
Zn  d  bands  overlap  the  Cu(d)-0(p)  bands  below  Ef .  In  general,  substitutions 
on  the  Cu(l)  site  give  sharper  peaks  in  the  density  of  states,  arising  from  the  d 
bands  of  the  substitutional  atom.  The  d  peaks  for  Cu(l)  substitutions  also  lie 
at  slightly  higher  energies  than  those  for  substitutions  on  the  Cu(2)  site. 

The  valences  for  VBa2Cu2-l/07  are  presented  in  Table  5.2,  with  the  valences 
for  M  =  Cu  from  Table  4.1  shown  for  comparison.  The  experimentally- 
determined  substitution  sites230 _236  are  also  indicated  in  Table  5.2.  There  is 
a  strong  dependence  on  the  substitution  site  for  both  the  metal  and  oxygen 
valences.  Substitution  on  the  Cu(l)  site  always  results  in  a  larger  valence  of  the 
doped  metal  than  for  substitution  on  the  Cu(2)  site.  The  large  dopant-metal 
valences  for  Fe  and  Co  substitution  may  be  due  to  the  neglect  of  correlation 
effects  within  the  d  orbitals,  which  might  alter  the  occupation  of  these  bands. 
Both  Fe(l,2)  and  Co(1.2)  substitution  give  a  substantial  decrease  in  the  valence 
of  the  copper  atom  on  site  Cu(2),  but  Ni  or  Zn  substitution  results  in  little 
change  in  the  valence  of  the  copper  atom  in  the  CuC>2  plane. 

The  density  of  states  results  for  Fe  and  Co  substitution  in  Table  5.3  show 
large  changes  in  p(Ef )  for  both  Cu(l)  and  Cu(2)  dopant  sites.  The  Fermi  energy 
is  pinned  near  the  peak  of  the  Fe  or  Co  d  state  energy,  so  that  p(Ep)  is  insensitive 
to  small  variations  of  the  parameters  e^(Fe)  or  ej(Co)  in  Table  5.1,  but  for  the 
same  reason  AEp  is  somewhat  sensitive  to  these  values.  The  Fermi  energy 
shifts  upwards  for  both  Fe  and  Co  substitution,  with  only  a  small  dependence 
on  the  particular  dopant  site.  These  substitutions  do  disrupt  the  character  of  the 
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Fig.  5.6.  Total  densities  of  states  for  YBa2Cu2Zn07. 
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Cu(d)-0(p)  antibonding  band  at  Ef,  even  when  located  on  the  Cu(l)  site,  which 
may  explain  the  decrease  in  Tc  seen  for  Fe  and  Co  substitution. 243,244,250-254 
Replacement  of  Cu  by  Ni  shows  a  smaller  change  in  p{Ep),  as  well  as  a  slight 
lowering  of  the  Fermi  energy.  This  may  explain  the  moderate  decrease  in  Tc 
observed  for  this  case,243-245  even  though  Ni  dopes  the  Cu(2)  site.251,256 

The  shift  in  Ep  in  Table  5.3  is  relatively  site  independent  for  all  the 
transition-metal  substitutions,  indicating  the  dopant  origin  of  the  shift,  while 
the  density  of  states  at  Ep  is  strongly  site  dependent,  except  for  the  case  of  Zn. 
Substitution  of  Zn  for  Cu  gives  only  small  changes  in  the  density  of  states  at  Ep , 
along  with  a  slight  increase  in  the  Fermi  energy,  in  contrast  to  the  large  decrease 
in  Tc  seen  experimentally.243,230  The  lack  of  resonance  of  the  Zn  d  states  with 
the  0  p  states  leads  to  a  smaller  fraction  of  Zn  d  character  at  the  Fermi  energy 
for  this  case.  Nevertheless,  the  present  results  indicate  that  the  direct  electronic 
effect  of  Zn  should  be  less  than  that  of.  e.g.,  Fe.  The  Zn  dopants  on  the  Cu(2) 
site  will  interfere  with  the  C'u02  conduction  bands,  but  additional  calculations 
(perhaps  including  correlation  effects)  may  be  needed  to  properly  describe  this 
case. 

YBa.\/Cu307,  M  =  Sr  or  La 

Substitutions  on  the  Ba  site  were  examined  in  YBa.UCusO?  for  M  —  Sr 
and  M  =  La.  For  the  case  of  Sr,  little  change  is  seen  in  direct  electronic  effects, 
as  shown  by  the  valences  in  Table  5.4.  The  shift  in  the  Fermi  energy  is  only 
-0.1  eV.  as  expected  for  the  substitution  of  divalent  Sr  for  divalent  Ba.  The 
total  density  of  states  at  the  Fermi  energy  is  3.2  states/eV  cell,  identical  within 
the  model  to  the  result  for  YBa2Cu307.  This  is  consistent  with  the  experimental 
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TABLE  5.4.  Valences  An  for  YBa-l/Cu307. 


Y 

Ba 

M 

Cu(l) 

Cu(2)° 

0(1) 

0(2)“ 

0(3)“ 

0(4)' 

M  = 

Sr 

1.77 

1.44 

1.41 

1.34 

0.98 

-1.12 

-1.09 

-1.11 

-1.20 

M  = 

Ba 

1.77 

1.44 

1.34 

0.98 

-1.12 

-1.09 

-1.11 

-1.21 

M  = 

La 

1.76 

1.44 

2.43 

1.30 

0.71 

-1.12 

-1.19 

-1.20 

-1.23 

“  Averaged  over  sites  in  unit  cell  for  M  =  Sr  and  La. 
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results  and  would  support  a  purely  structural  effect  for  the  observed  decrease  in 
the  superconducting  transition  temperature  of  YBa2_xSrxCu307. 237,238 

Substitution  of  La  for  Ba  does  change  the  atomic  valences,  as  shown  in 
Table  5.4,  resulting  in  0.46  fewer  holes  within  each  Cu02  plane  compared  with 
the  undoped  material  YBa2Cu307-  The  additional  electron  of  La  raises  the 
Fermi  energy  by  A Ef  =  0.31  eV.  There  is  also  a  slight  decrease  in  p(£V)  to 
2.8  states/eV  cell.  We  propose  that  the  dominant  effect  of  La  substitution  on 
the  superconducting  properties  of  YBa2_xLaxC'u307  is  the  destruction  of  the 
hole  charge  carriers  within  the  Cu02  planes. 

YBa2Cu306^,  L  =  X  or  F 

Substitution  of  X  or  F  for  O  is  complicated  by  the  four  crvstallographicallv 
distinct  oxygen  sites,  along  with  the  many  possible  combinations  for  fractional 
or  multiple-site  substitutions.  These  atoms  might  also  be  accommodated  in  the 
structure  at  the  normally- vacant  sites  that  separate  the  CuO  chains.264  However, 
we  find  that  sites  0(2)  and  0(3)  in  the  Cu02  planes  behave  quite  similarly  in 
their  electronic  structure,  as  expected  from  the  minor  orthorhombic  distortion 
of  the  structure.17-21  We  consider  only  single-site  substitution  of  nitrogen  or 
fluorine  for  oxygen,  and  neglect  the  possibility  of  interstitial  sites  or  fractional 
occupations  for  the  substitutional  atom. 

Table  5.5  gives  the  valences  for  the  oxygen  replacements,  with  the  valences 
for  YBa2Cu307  listed  for  comparison.  The  largest  valence  effect  occurs  on  the 
copper  site  nearest  the  substitutional  atom.  The  substitutions  X(l)  and  F(l) 
both  lead  to  a  significant  decrease  in  the  valence  of  copper  on  the  Cu(l)  site, 
while  F(2)  or  F(4)  substitutions  cause  a  large  decrease  in  the  Cu(2)  valence. 


TABLE  5.5.  Valences  An  for  YBa2Cu306L. 


L 

Cu(l) 

Cu(2)a 

0(1) 

0(2)a 

0(3)° 

0(4)“ 

N(l) 

1.18 

0.97 

-0.89 

-1.06 

-1.08 

-1.25 

N(2) 

1.24 

1.06 

-1.13 

-1.05 

-1.05 

-1.25 

\(4) 

1.30 

0.98 

-1.14 

-1.08 

-1.09 

-1.14 

0 

1.34 

0.98 

-1.12 

-1.09 

-1.11 

-1.21 

F(l) 

0.81 

1.00 

-0.67 

-1.08 

-1.07 

-1.28 

F(2) 

1.36 

0.73 

-1.11 

-0.86 

-1.16 

-1.23 

F(4) 

1.26 

0.78 

-1.15 

-1.16 

-1.18 

-0.87 

Averaged  over  sites  in  unit  cell. 
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TABLE  5.6.  Shift  in  Ep  and  density  of  states  for  YBa2Cu30eL. 


L 

(eV) 

p(ef) 

(states/eV  cell) 

N(l) 

0.02 

5.0 

N(2) 

-0.06 

3.7 

N(4) 

-0.08 

4.8 

0 

3.2 

F(l) 

0.11 

3.5 

F(2) 

0.14 

3.0 

F(4) 

0.25 

2.8 
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The  results  for  N(2)  substitution  do  show  a  slight  increase  in  the  number  of 
holes  in  the  Cu02  planes.  However,  we  do  not  expect  these  small  shifts  to  give 
a  significant  enhancement  of  the  superconducting  properties  for  YBa2Cu30eN. 

Table  5.6  gives  the  shift  in  the  Fermi  energy  and  p{Ep)  for  the  replacement 
of  oxygen.  We  note  that  X  acts  primarily  to  lower  the  Fermi  energy,  while  F 
raises  it;  i.e.,  X  acts  as  an  acceptor  and  F  as  a  donor  when  replacing  oxygen. 
Fluorine  is  found  to  have  a  larger  effect  on  Ep  than  nitrogen,  but  F  substitution 
has  little  effect  on  p(Ep).  The  X  substitutions  all  show  increases  in  p(Ep)  within 
the  present  model. 

BiPbCaSr2Cu208 

We  now  consider  the  effect  of  periodic  substitution  of  Pb  for  one  Bi  atom 
within  the  crystal  structure  of  Bi2CaSr2Cu20g.  We  again  neglect  the  incom¬ 
mensurate  structural  modulation48-30  and  the  disorder  in  the  BiO  layers48-32-34 
associated  with  this  phase.  Although  Pb  substitution  may  also  stabilize  a  dif¬ 
ferent  crystal  structure,30,6 1,62  we  confine  this  calculation  to  the  same  model 
structure  that  was  used  for  Bi2CaSr2Cu20s  in  Chapter  IV. 

The  electronic  energy  bands  for  BiPbCaSr2Cu208  are  shown  in  Fig.  5.7. 
The  Fermi  energy  shifts  by  —0.19  eV  compared  with  that  of  Bi2CaSr2Cu20s, 
consistent  with  the  removal  of  one  electron  per  formula  unit  with  Pb  substitu¬ 
tion.  The  Pb  p  bands  lie  well  above  Ep,  while  the  remaining  Bi  p  band  crosses 
below  Ep  to  only  -0.3  eV  at  D,  so  that  part  of  the  electron  reservoir  is  forced 
back  into  the  remaining  structure.  However,  the  calculated  valences  in  Table  5.7 
show  that  the  hole  donated  when  Pb  substitutes  for  Bi  more  than  compensates 
for  this  addition  of  electrons,  resulting  in  a  net  increase  of  0.29  holes  per  Cu02 
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Fig.  5.7.  Electronic  energy  bands  for  BiPbCaSrgCu208. 
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TABLE  5.7.  Valences  for  BLVCaS^Ci^Og. 


Bi 

Pb 

Ca 

Sr" 

Cua 

0(1)- 

0(2)a 

0(3)a 

M  = 

=  Bi 

1.80 

1.51 

1.09 

0.84 

-1.09 

-1.14 

-1.18 

M  = 

=  Pb 

1.79 

1.23 

1.51 

1.09 

1.01 

-1.03 

-1.14 

-1.17 

"Averaged  over  sites  in  unit  cell  for  M  =  Pb. 
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Fig.  5.8.  Local  densities  of  states  for  the  meta! 

atoms  in  BiPbCaSroCuo0a. 
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Fig.  5.9.  Local  densities  of  states  for  the  oxygei 
atoms  in  BiPbCaSr.Cu.0_. 
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Fig.  5.10.  Total  densities  of  states  for  BiAf^aSr3Cu208. 
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layer  in  each  unit  cell. 

The  local  densities  of  states  for  BiPbCaSr2Cu20g  are  shown  in  Figs.  5.8  and 
5.9.  The  results  are  very  similar  to  those  for  Bi2CaSr2Cu20g  in  Figs.  4.11  and 
4.12,  with  an  overall  shift  in  the  Fermi  energy  of  i\Ep  =  —0.19  eV  and  an  upward 
shift  of  the  Pb  local  density  of  states  with  respect  to  Bi.  The  total  densities  of 
states  for  Bi2CaSr2Cu20g  and  BiPbCaSr2Cu20g  are  shown  in  Fig.  5.10.  Within 
the  present  model,  the  density  of  states  at  the  Fermi  energy  increases  with  Pb 
doping,  from  2.4  states /eV  cell  to  4.1  states/eV  cell.  We  propose  that  the 
primary  effect  of  Pb  substitution  is  to  increase  the  carrier  density  within  the 
superconducting  copper-oxide  regions. 

TU/CaBa2Cu208,  M  =  Hg  or  Pb 

Let  us  now  consider  the  modification  of  the  electronic  structure  when 
T1  is  partially  replaced  by  Hg  or  Pb,  focusing  on  the  representative  material 
Tl2CaBa2Cu20g.  The  electronic  energy  bands  for  TlHgCaBa2Cu20g  are  shown 
in  Fig.  5.11.  Only  the  hole  conduction  bands  cross  Ep  in  this  material,  with 
the  T1  s  and  Hg  s  states  remaining  unoccupied.  The  local  density  of  states  for 
Hg  is  shown  in  Fig.  5.12,  with  the  d  states  at  about  —6  eV.  The  local  densities 
of  states  for  the  other  atoms  in  TlHgCaBa2Cu20g  are  very  similar  to  those  in 
Figs.  4.17  and  4.18,  but  with  a  shift  in  the  Fermi  energy  of  A Ep  =  -0.23  eV. 
The  valences  given  in  Table  5.8  indicate  that  Hg  is  much  less  ionic  than  T1  in 
this  material.  The  added  hole  of  Hg  leads  to  a  lowering  of  Ep ,  and  in  the  present 
model  there  is  an  increase  in  p{Ep )  from  2.7  states/eV  cell  to  4.2  states  eV  cell. 
As  can  be  seen  in  Table  5.8,  there  is  a  positive  shift  in  the  valence  of  both  Cu 
and  0(1)  in  the  charge-carrying  Cu02  planes,  indicating  that  hole  carriers  are 
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Fig.  5.11.  Electronic 
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TABLE  5.8.  Valences  An  for  Tl„l/CaBa2Cu2  08. 


T1 

M 

Ca 

Baa 

Cua 

0(1)“ 

0(2)a 

0(3)a 

M  = 

Hg 

1.02 

0.50 

1.45 

1.38 

1.14 

-1.03 

-1.10 

-0.83 

M  = 

T1 

1.05 

1.45 

1.37 

0.86 

-1.13 

-1.03 

-0.71 

M  = 

Pb 

1.06 

1.46 

1.45 

1.37 

0.66 

-1.21 

-0.98 

-0.62 

a  Averaged  over  sites  in  unit  cell  for  M  =  Hg  and  Pb. 
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doped  into  the  superconducting  region  by  Hg  substitution. 

A  similar  calculation  for  Pb  replacing  T1  finds  that  Pb  raises  the  Fermi 
energy  by  A Ep  =  0.19  eV.  The  valences  of  Table  5.8  for  TlPbCaBa2Cu20a 
demonstrate  that  doping  with  Pb  significantly  decreases  the  valences  of  the  Cu 
and  0(1)  sites  in  the  Cu02  planes.  The  local  density  of  states  for  Pb  is  shown 
in  Fig.  5.12.  The  Pb  d  states  lie  about  18  eV  below  Ep ,  while  the  dispersive  p 
states  are  unoccupied.  The  total  density  of  states  at  the  Fermi  energy  is  lower 
for  this  case:  p(Ep )  =  1.9  states /eV  cell. 

Summary  of  Substitution  Effects 

In  summary,  calculations  for  the  replacement  of  Cu,  Ba,  and  0  with  other 
elements  have  been  performed  for  the  system  YBa2Cu307.  The  results  are  in 
reasonable  agreement  with  expected  chemical  trends.  The  suppression  of  Tc 
by  A1  substitution  on  the  Cu(l)  site  is  associated  with  the  destruction  of  the 
conduction  band  within  the  CuO  chain  and  a  more  ionic  character  for  A1  and 
0  within  the  chain.  The  substitutions  of  the  transition  metals  Fe,  Co,  Ni,  and 
Zn  for  Cu  also  give  results  that  are  compatible  with  experiment,  except  in  the 
case  of  Zn.  However,  the  experimentally-observed  changes  in  superconducting 
properties  may  result  from  structural  changes  that  are  not  considered  here. 
Substitutions  on  the  Cu(2)  site  disrupt  the  CuC>2  plane  conduction  bands 
more  than  do  substitutions  on  the  Cu(l)  site.  Results  for  the  replacement 
of  Ba  by  Sr  give  support  for  a  purely  structural  suppression  of  Tc ,  while  we 
suggest  that  La  substitution  leads  to  destruction  of  the  hole  charge  carriers 
by  the  donor  electrons,  resulting  in  suppression  of  superconductivity.  The 
results  for  F  and  X  substitutions  on  the  oxygen  sites  show  that  the  electronic 
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structure  of  YBa2Cu3  07  may  be  affected  by  such  replacements,  but  no  obvious 
enhancements  of  the  superconducting  properties  are  found. 

Substitution  of  Pb  for  Bi  in  BiA/CaSr2Cu2C>8  gives  an  increase  in  the 
hole  carriers  within  the  Cu02  planes,  assuming  an  isostructural  form  for 
BiPbCaSr2Cu208  and  Bi2CaSr2Cu20s.  The  most  obvious  effect  of  doping  with 
Hg  or  Pb  in  Tl.\/CaBa2Cu208  is  the  shift  in  the  Fermi  energy,  with  Hg  creating 
holes  and  Pb  destroying  them. 
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CHAPTER  VI 

VACANCY  EFFECTS 

Using  the  tight-binding  model  of  Chapter  III,  we  now  consider  the  electronic 
effects  of  isolated  oxygen  vacancies  in  Lai .ssSro.isCuC^-y  and  YBa2Cu307_y, 
and  of  lanthanum  vacancies  in  La2_xCu04.  We  neglect  local  strains  and  the 
small  changes  in  the  lattice  constants  that  occur  as  vacancies  are  introduced. 
For  simplicity,  we  use  the  bet  structure9  for  La!  85Sr0.i5CuO4,  neglecting  the 
small  orthorhombic  distortion.  For  YBa2Cu307,  the  observed  orthorhombic 
structure18  is  used.  The  change  in  the  density  of  states  is  computed  with  a 
Green’s  function  method. 

As  described  in  Chapter  II,  oxygen  vacancies  have  a  significant  effect 
on  the  electronic  properties  of  La2_xSrxCu04_y  and  YBa2Cu307_y.  In  the 
simplest  picture  of  a  rigid-band  model,  the  removal  of  oxygen  is  expected  to 
add  electrons,  giving  the  opposite  effect  of  doping  with  holes.  Indeed,  both 
La2 -xSrxCu04_y  and  YBa2Cu307_y  exhibit  a  depression  in  Tc  with  increasing 
oxygen  vacancies.3, 6,7,26,27  However,  the  lead  cuprates  Pb2-4i_xCaxSr2Cu308-ty 
(A  =  Y  or  a  rare  earth  element)  show  little  change  in  electronic  properties  as 
oxygen  is  added  to  or  removed  from  the  Cu  layer  intercalated  between  the 
PbO  layers,83  while  the  electron-doped  materials  such  as  Nd2_xCexCu04_y 
have  improved  superconducting  properties  with  the  introduction  of  oxygen 

*  ft  ft 

vacancies. 

The  magnetic  properties  of  these  materials  are  also  strongly  affected  by  oxy¬ 
gen  vacancies.142,150  Electronic  structure  calculations  within  the  local-density 
approximation  for  La2Cu04_y  have  been  interpreted  as  indicating  that  oxy- 
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gen  vacancies  may  lead  to  antiferromagnetism  through  narrowing  of  the  valence 
bands  in  the  neighborhood  of  Ep,26S  although  the  LDA  may  be  inadequate  for 
a  proper  treatment  of  magnetic  effects.269  LDA  calculations  for  Yl^CusOe  do 
not  properly  describe  the  antiferromagnetism  in  this  material.270,271 


Green’s  Function  Technique 

The  unperturbed  Green’s  function  from  (3.7)  is 

v(k.n)  v^(k, n) 


Go  (£•)  =  y>f 


k  ,n 


E  -  E(k,n)  -  iS 


(6.1) 


In  the  tight-binding  representation,  the  Green’s  function  is  a  .V  x  _V  matrix, 
Go(E),  over  the  .Y  basis  states  of  the  system.  For  a  perturbing  potential  V,  the 
total  density  of  states  p(E)  is  changed  by272,273 


-  V-W • 


(6.2) 


where  q(£)  is  the  phase  shift 


'?(£) 


Im  Log  det  Y  -  VGo(-E')!  • 


(6.3) 


Here  I  is  the  .V  x  .V  identity  matrix.  Now  suppose  the  perturbing  potential  V 
acts  only  on  a  particular  n  x  n  subspace  of  states  in  the  system,  e.g.,  the  n  states 
associated  with  an  isolated  atomic  vacancy.  Using  the  subspace  Hamiltonian 
technique,2'4  we  partition  the  Green’s  function  such  that  the  perturbation  V 
has  nonzero  elements  only  in  the  n  x  n  subspace.  The  determinant  of  (6.3)  is 
then 


det  T  -  VG0(£); 
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and  thus 


det  T  -  VG0(£)]  =  det  [1  -  l'G,s0ub(£)j  . 


The  matrices  1,  V',  and  G'gub(£')  are  thus  reduced  in  size  tonxn  for  the  n  states 
in  the  perturbation  subspace,  substantially  reducing  the  numerical  effort  for  the 
problem. 

For  the  particular  case  of  an  atomic  vacancy  calculation,  the  perturbing 
potential  V  may  be  approximated2 13  by  letting  the  diagonal  elements  of  V 
approach  —  oc,  so  that  for  a  full  atomic  vacancy  the  perturbed  states  are 
completely  unoccupied.  The  elements  of  the  identity  matrix  may  then  be 
neglected,  so  that  (6.5)  reduces  to 

det  I  -  VG0(£)i  =s  det  -FG's0ub(£)] 

(6.6) 

*  det  (-V)  det  Gs0ub(£)  . 

Within  this  approximation  the  perturbation  V  is  independent  of  the  energy,  so 
substituting  (6.6)  into  (6.2)  and  (6.3)  yields  the  change  in  the  total  density  of 


states 


2  d 


A p(E)  =  -  -  —  Im  Log  det  GJub(£) 


2  d  ,  Im  det  Gs0  £  \ 

= - —  arctan  - — -  . 

7 T  dE  \  Re  det  G'oub(£)  ’ 

The  second  form  in  (6.7)  is  used  for  the  calculations.  The  principal  value  of 


the  arctangent  is  used,  but  with  shifts  of  ±27r  added  to  preserve  continuity  in 
the  physical  result.  The  energy  derivative  is  calculated  numerically.  For  the  full 
vacancy  result,  a  simple  check  of  the  calculation  is  that  the  integral  of  A p(E) 
gives  the  number  of  states  (for  both  spins)  associated  with  the  vacancy  site: 


—  -X) 

J  ME) 


d  E  =  2  n  . 
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The  modified  total  density  of  states  for  y  atomic  vacancies  is  given  by 

p{E)  =  p0(E)  +  yAp(E)  ,  (6.9) 

with  the  unperturbed  total  density  of  states  po(E)  given  by  (3.6).  The  Fermi 
energy  is  calculated  by  integrating  p(E)  up  to  the  number  of  valence  electrons 
for  a  given  concentration  of  vacancies.  This  calculation  neglects  the  interaction 
between  vacancies,  and  is  strictly  valid  only  for  a  small  fraction  of  vacancies. 
However,  for  y  =  1  we  have  performed  an  independent  calculation  with  one 
oxygen  vacancy  per  formula  unit  for  Lai.gsSro.isCuOs  and  for  YBa2Cu306.  and 
the  change  in  the  density  of  states  is  approximately  the  same.276 

Oxygen  Vacancies  in  La!  85Sr0.i5CuO4_y 

We  consider  both  oxygen  sites  in  Laj.gsSro.isCuO*,  with  0(1)  in  the  C’u02 
plane  and  0(2)  in  the  LaO  layer.  The  perturbation  affects  only  the  4x4  subspace 
associated  with  the  s,  pr,  py,  and  p.  states  of  the  oxygen  atom.  Figure  6.1 
shows  the  change  in  the  density  of  states  Ap(E)  for  a  single  oxygen  vacancy  in 
Lat.soSro.isCuO^y .  Notice  that  Ap  has  more  structure  near  the  Fermi  energy 
for  the  0(1)  site  than  for  the  0(2)  site,  because  of  the  strong  pdcr  interactions 
of  0(1)  p  orbitals  with  copper  d  orbitals.  For  both  sites,  A p  is  negative  at  the 
unperturbed  Fermi  energy.  The  larger  magnitude  at  E°F  for  the  0(1)  site  is 
consistent  with  the  observation  of  substantial  in-plane  0  pz  and  py  character 
at  the  Fermi  energy.168  Figure  6.2  shows  the  resulting  total  density  of  states  for 
y  oxygen  vacancies  per  formula  unit  on  the  0(1)  site.  A  substantial  decrease 
is  seen  in  the  total  density  of  states  within  2  eV  below  Ep  as  y  increases.  As 
shown  in  Table  6.1,  the  density  of  states  at  Ep  is  smaller  than  the  y  =  0.0  value 


Energy  (eV) 


Fig.  6.1.  Change  in  density  of  states  for  an  isolated 
oxygen  vacancy  in  La}  BSSrQ  19Cu04_y.  Oxygen  site  0(1) 
is  in  the  Cu02  plane,  and  0(2)  is  in  the  La/Sr  layer. 
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TABLE  6.1.  Shift  in  Ep  and  density  of  states  for  y  oxygen  vacancies 
per  formula  unit  on  site  0(1)  in  Lai.ssSro.isCuC^-y . 


y 

A  Ef 

(eV) 

p(Ef) 

(states/eV  cell) 

0.0 

0.00 

1.9 

0.1 

0.07 

1.3 

0.2 

0.20 

0.7 

0.3 

0.44 

0.6 

0.4 

0.65 

1.0 

0.5 

0.76 

1.5 

0.6 

0.83 

2.0 

0.7 

0.88 

2.3 

0.8 

0.93 

2.3 

0.9 

0.97 

2.3 

1.0 


1.01 


2.3 
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TABLE  6.2.  Shift  in  Ep  and  density  of  states  for  y  oxygen  vacancies 
per  formula  unit  on  site  0(2)  in  Lai.85Sro.isCuC>4_y . 


y 

A  EF 

(eV) 

p(ef) 

(states/eV  cell) 

0.0 

0.00 

1.9 

0.1 

0.12 

1.3 

0.2 

0.31 

0.7 

0.3 

0.60 

0.6 

0.4 

0.87 

0.9 

0.5 

1.03 

1.4 

0.6 

1.14 

1.8 

0.7 

1.26 

1.5 

0.8 

1.41 

0.9 

0.9 

1.66 

0.6 

1.0 


1.94 


0.8 
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for  y  <  0.5,  and  Ep  monotonically  increases  with  y.  The  0(1)  vacancies  thus 
behave  as  donors,  filling  the  hole  states  near  the  Fermi  energy. 

For  vacancies  on  the  0(2)  site,  Table  6.2  indicates  that  p(Ep )  is  smaller 
than  the  y  =  0.0  value  for  all  y  >  0.  The  Fermi  energy  increases  somewhat 
faster  for  0(2)  vacancies  than  for  0(1)  vacancies,  since  the  more  ionic  nature  of 
the  Cu-0(2)  and  the  La-0(2)  bonds  leaves  most  of  the  0(2)  character  below 
Ep.  These  0(2)  vacancies  are  again  electron  donors. 

Lanthanum  Vacancies  in  La2_xCu04 

The  perturbation  subspace  for  a  La  vacancy  contains  the  six  states  (one 
s  an  I  five  d)  associated  with  lanthanum.  The  change  in  the  density  of  states 
of  La2_xCu04  for  a  single  La  vacancy  is  shown  in  Fig.  6.3.  Little  change  is 
seen  at  Ep,  because  of  the  ionic  and  insulating  nature  of  the  La O  layer.  The 
large  contribution  near  5  eV  arises  from  the  unoccupied  La  d  bands.  The 
total  density  of  states  for  ,r  lanthanum  vacancies  per  formula  unit  is  shown  in 
Fig.  6.4.  The  density  of  states  near  Ep  is  relatively  unchanged,  with  only  a 
rigid  shift  of  the  bands  observed  with  increasing  x.  As  shown  in  Table  6.3,  Ep 
is  lowered  as  La  vacancies  are  introduced.  Lanthanum  vacancies  thus  donate 
holes  to  the  material,  just  as  divalent  Sr  does  when  substituted  for  trivalent  La. 
This  may  explain  the  enhancement  of  superconductivity11,12  observed  in  some 
samples  of  La2Cu04,  in  that  deficiencies  of  La  within  the  undoped  material  can 
provide  the  required  hole  carriers.  In  our  calculations,  there  is  an  increase  in 
p(Ep)  for  x  <  0.3. 

A  recent  set  of  experiments  has  shown  that  excess  oxygen  can  be 
incorporated  into  the  structure  of  La2Cu04+y  while  maintaining  full  La 
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Fig.  6.3.  Change  in  density  of  states  for  a  lanthanum 
vacancy  in  La  CuO  . 
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Fig.  6.4.  Density  of  states  for  x  lanthanum  vacancies 
in  La„  CuO„. 

2-i  4 
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TABLE  6.3.  Shift  in  Ep  and  density  of  states  for  x  lanthanum 
vacancies  per  formula  unit  in  La2_xCu04. 


X 

A  EF 

(eV) 

p(ef) 

(states/eV  cell) 

0.0 

0.00 

1.4 

0.1 

-0.11 

1.9 

0.2 

-0.20 

1.9 

0.3 

-0.29 

1.7 

0.4 

-0.42 

1.2 

0.5 

-0.56 

1.2 

0.6 

-0.69 

1.5 

0.7 

-u.79 

2.0 

0.8 

-0.86 

2.4 

0.9 

-0.92 

2.9 

1.0 

-0.98 

3.2 
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stoichiometry.277'279  The  excess  oxygen  ( y  ^  0.1)  appears  to  cluster  in  the 
superoxide  ion  0^  within  the  material.278  This  oxygen-rich  material  forms  a  sec¬ 
ond  orthorhombic  phase  below  320  K  (with  a  slightly  larger  c  lattice  constant), 
separating  from  the  non-superconducting  phase  of  stoichiometric  La2Cu04.279 
The  oxygen-rich  phase  exhibits  bulk  superconductivity  at  Tc  —  34  K  in  30% 
of  the  sample  volume  for  y  =  0.13. 278  These  experiments  indicate  that  doping 
with  oxygen  does  in  fact  lead  to  bulk  superconductivity  in  La2Cu04.  However, 
to  test  the  results  of  our  calculations,  further  experiments  that  in+mduce  La 
vacancies  while  fixing  the  oxygen  stoichiometry  are  needed. 

Oxygen  Vacancies  in  YBa2Cu30r_y 

We  now  consider  the  electronic  effects  of  oxygen  vacancies  in  YBa2Cu307. 
We  examine  only  the  oxygen  chain  site  0(1),  which  is  experimentally  deter¬ 
mined  to  be  the  primary  site  for  the  introduction  of  oxygen  vacancies.24  The 
perturbation  subspace  is  again  4x4  for  oxygen  vacancies. 

The  change  in  the  density  of  states  of  YBa2Cu3  07_y  for  an  isolated 
oxygen  vacancy  introduced  on  the  0(1)  site  is  shown  in  Fig.  6.5.  The  largest 
contributions  to  A p{E)  are  well-removed  from  the  Fermi  energy,  but  a  negative 
contribution  is  found  at  E°F.  Figure  6.6  shows  the  resulting  densities  of  states 
for  y  oxygen  vacancies  per  formula  unit.  We  find  peaks  in  the  density  of  states 
that  correspond  to  valence  band  photoemission  results  for  y  <  0.3,  although 
the  experimental  binding  energies  are  about  0.5  eV  higher. 167,180,189,191  As 
the  number  of  oxygen  vacancies  y  increases,  p(Ep)  decreases  monotonically, 
consistent  with  a  calculation  within  the  LDA  for  YBa2Cu306.28°  Table  6.4 
shows  the  calculated  shifts  in  Ep  and  p(Ef)  for  0.0  <  y  <  1.0.  The  removal  of 


126 


Energy  (eV) 


Fig.  6.5.  Change  in  density  of  states  for  a  single 
oxygen  vacancy  on  the  0(1)  chain  site  in  YBa8Cu307  y. 
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TABLE  6.4.  Shift  in  EF  and  density  of  states  for  y  oxygen  vacancies 
per  formula  unit  on  site  0(1)  in  YBa2Cu307_y. 


y 

A  Ep 

(eV) 

p(Ef) 

(states/ eV  cell) 

0.0 

0.00 

3.2 

0.1 

0.01 

3.1 

0.2 

0.03 

3.0 

0.3 

0.05 

2.8 

0.4 

0.07 

2.7 

0.5 

0.10 

2.4 

0.6 

0.13 

2.2 

0.7 

0.16 

2.0 

0.8 

0.21 

1.8 

0.9 

0.26 

1.6 

1.0 

0.31 

1.5 
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oxygen  monotonically  raises  the  Fermi  level,  so  the  present  calculation  confirms 
the  expectation  that  oxygen  vacancies  act  as  electron  donors  in  this  system. 
The  decrease  in  the  density  of  states  at  Ep  is  consistent  with  shifts  of  the 
valence  bands  seen  experimentally  in  photoemission  studies  of  oxygen-deficient 
YBa2Cu307_y.194 

Within  this  model,  we  do  not  find  a  plateau  in  p(Ep )  with  respect  to  y  that 
might  correspond  to  the  change  in  Tc  with  oxygen  vacancies.26'27  However,  we 
note  that  we  have  neglected  ordering  of  the  oxygen  vacancies,  which  may  play 
a  role  in  this  efFect. 28,29  Also,  recent  experiments  show  that  the  electrons  added 
with  oxygen  vacancies  first  localize  on  the  CuO  chains,  with  additional  vacancies 
causing  charge  transfer  from  the  chains  to  the  Cu02  planes.32,281  The  details 
of  this  charge  transfer  may  be  responsible  for  the  plateau  in  Tc  with  oxygen 
vacancy  concentration.282  Oxygen-deficient  YBa2Cu307_y  has  been  shown  to 
be  single  phase  for  0.0  £  y  <  1.0,  so  that  a  single  Tc  is  associated  with  each 
oxygen  stoichiometry.283 

Summary  of  Vacancy  Effects 

Using  the  tight-binding  model  of  Chapter  III,  we  have  calculated  the 
electronic  effects  of  oxygen  vacancies  in  La185Sr0.i5CuO4_y  and  YBa2Cu307_y, 
and  of  lanthanum  vacancies  in  La2_xCu04.  We  find  that  oxygen  vacancies  act  as 
electron  donors  in  both  systems,  raising  the  Fermi  energy  and  decreasing  p(Ep ) 
for  y  £  0.5.  The  concentration  of  holes  is  thus  decreased,  consistent  with  the 
observations  that  oxygen  vacancies  degrade  the  superconductivity  and  metallic 
conductivity  in  these  materials. 

Lanthanum  vacancies  in  La2_xCu04  donate  holes,  lowering  the  Fermi 
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energy  and  increasing  p(Ep)  for  x  <  0.3.  Small  concentrations  of  La  vacancies 
may  thus  have  the  same  effect  electronically  as  divalent  metal  dopants  such 
as  Sr,  leading  to  enhancement  of  superconductivity  in  nominally  stoichiometric 
La2Cu04- 
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CHAPTER  VII 

EXCITONIC  MECHANISM  FOR  SUPERCONDUCTIVITY 

Although  many  mechanisms  have  been  proposed  to  explain  high-Tc  super¬ 
conductivity,  no  theory  has  yet  successfully  described  the  properties  of  these 
materials.284  In  this  chapter  we  examine  a  specific  version  of  the  generic  ex- 
citonic  mechanism  for  superconductivity  originally  proposed  by  Little,285  and 
elaborated  bv  Ginzburg286  and  by  Allender,  Bray,  and  Bardeen.287  Our  mecha¬ 
nism  requires  a  specific  structural  and  electronic  configuration,  namely  metallic 
planes  (containing  the  charge  carriers)  adjacent  to  insulating  metal-oxide  layers 
such  as  LaO,  BaO,  or  NdO.  The  insulating  layers  must  support  virtual  excita¬ 
tions  from  oxygen  p  states,  below  the  Fermi  energy,  to  metal  d  states,  above 
Ep.  The  effective  pairing  potential  of  the  charge  carriers  is  mediated  by  the 
interaction  of  the  carriers  with  these  excitations. 

One  advantage  of  the  present  mechanism  is  its  applicability  to  both  the 
copper-oxide  and  the  bismuth-oxide  superconductors,  since  no  specific  mag¬ 
netic  properties  are  required.  For  example,  the  systems  Bai_xKxBi03  and 
BaPbi_yBiy03  are  diamagnetic,160  while  both  the  hole-doped  and  the  electron- 
doped  cuprates  generally  display  antiferromagnetic  order  in  their  nonsuper¬ 
conducting  stoichiometries. 85,142-145  150, 157-159  It  is  thus  unlikely  that  a  mag¬ 
netic  mechanism  for  high-temperature  superconductivity  could  apply  to  both 
the  copper-oxide  and  the  bismuth-oxide  materials.  Also,  as  found  below,  the 
exciton-mediated  pairing  potential  is  independent  of  the  sign  of  the  charge  car¬ 
riers,  so  that  the  present  mechanism  can  describe  both  the  n-type  and  p-type 
cuprates.  Other  excitonic  mechanisms  (with  substantial  differences  from  each 
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other  and  from  the  present  mechanism)  have  been  recently  proposed  to  describe 
the  high-temperature  superconductors.288-298 

Electronic  and  Structural  Requirements 

The  present  excitonic  mechanism  for  high-temperature  superconductivity 
was  suggested  by  our  electronic  structure  calculations  for  85Sro.i5Cu04  and 
for  YBa2Cu307,  presented  in  Chapter  IV,  and  strengthened  by  subsequent 
calculations  for  the  other  high- 1),  materials.  In  every  material  examined,  there 
are  metallic  layers,  such  as  the  CuC>2  planes  in  all  the  cuprates  and  the  BiC>2 
layers  in  Ba0.6 Ko.4 BiC>3  and  BaPbo.75Bio.25O3.  Surrounding  these  metallic 
regions  are  metal-oxide  layers  such  as  LaO,  BaO,  SrO,  PrO,  NdO,  or  SmO. 
In  Fig.  7.1  we  show  the  local  densities  of  states  for  these  metal-oxide  layers 
found  in  the  results  of  Chapter  IV  for  the  hole-doped  superconductors.  The 
dashed  curve  peaking  below  Ep  shows  the  local  density  of  states  for  the  oxygen 
atom  in  the  metal-oxide  layer  e.g.,  0(4)  in  YBa2Cu3  07j.  The  region  within 
^  8  eV  below  Ep  is  primarily  composed  of  0  p  states.  The  solid  curve  peaking 
above  Ep  gives  p(E)  for  the  metal  atom  in  this  layer,  with  the  largest  peak 
arising  from  the  unoccupied  d  states.  In  every  case,  these  layers  are  insulating, 
in  the  sense  that  p(Ep )  is  very  small.  Similar  local  densities  of  states  have  been 
calculated  for  the  NdO  layers  in  the  electron-doped  material  Nd2-xCexCu04.297 
The  atomic  valences  presented  in  Chapter  IV  show  that  the  metal-oxide  layer 
is  also  quite  ionic  in  each  material,  with  significant  charge  transfer  between  the 
metal  and  oxygen  atom. 

These  ionic  and  insulating  layers  thus  can  support  electronic  excitations 
from  the  occupied  oxygen  p  states  to  the  vacant  metal  d  states,  with  character- 
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Fig.  7.1.  Local  densities  of  states  for  oxygen  (dashed  curves) 
and  metal  atoms  (solid  curves)  in  the  insulating  layers. 
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istic  excitational  energies  of  5  —  10  eV.  The  possibility  of  these  band  excitations 
forming  bound  electron-hole  pairs  (excitons)  is  described  in  the  next  section.  As 
in  Refs.  285-287,  the  excitations  are  spatially  separated  from  the  charge  carri¬ 
ers.  However,  the  insulating  layers  must  lie  in  close  proximity  to  the  metallic 
planes  (as  in  the  sandwich  structures  of  the  high-temperature  superconductors) 
to  produce  a  pairing  potential  sufficient  for  superconductivity. 


Formation  of  the  Exciton 


One  approach  is  to  treat  an  exciton  within  an  insulating  metal-oxide  layer  as 
a  Frenkel-like,  tightly-bound  charge  transfer  exciton,  using  the  Green’s  function 
method  that  was  applied  in  the  atomic  vacancy  calculations  in  Chapter  VI.  We 
model  the  exciton  as  a  state  that  is  split  off  from  the  one-electron  band  structure 
by  a  perturbation  Veh  given  by  the  electron-hole  Coulomb  interaction 


f'e  h  =  - 


ed0 


(7.1) 


where  do  is  the  metal-oxygen  bond  length  and  e  is  a  dielectric  constant  that 
characterizes  the  region  between  the  metal  and  oxygen  atoms.  Since  the  metal 
d  states  are  initially  unoccupied,  there  is  no  additional  electron-electron  repulsive 
term  U  as  might  be  associated  with  occupied  metal  states. 

Beginning  with  the  unperturbed  Green’s  function  Gq(E)  in  (3.7),  we 
consider  the  perturbation  to  act  only  on  the  subspace  of  states  associated  with 
the  metal  site,  so  that  the  perturbed  Green’s  function  G(E)  is  obtained  from 
the  one-electron  Dyson’s  equation 


C, 1(E)  =  Glub(E)  -  (E)VG(E)  , 


(7.2) 
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where  all  the  matrices  are  5x5  (representing  the  five  d  states  on  the  metal 
site),  and  V  is  a  diagonal  matrix  with  elements  Veh  applied  to  the  d  states  of 
the  metal  site.  We  invert  (7.2)  to  the  form 

G{E)  =  [1  -  G'pub(  E)V]~1Gs0'lh(E)  (7.3) 


for  the  numerical  calculations.  The  density  of  states  for  the  exciton  is  then  found 
from  (3.6) 


p{E)  =  -  -Trim  G(E)  .  (7.4) 

7 r 

where  the  trace  includes  only  the  d  states  on  the  metal  site. 

The  Frenkel  exciton  states  are  calculated  for  the  representative  materi¬ 
als  Lai.85Sr0.15CuO4  and  YBa2Cu3C)7.  From  the  crystal  structures  of  these 
materials,918  the  metal-oxygen  bond  lengths  are  2.73  A  for  La-0(2)  and  2.76  A 
for  Ba-0(4).  A  simple  estimate  for  the  static  dielectric  constant  e  is  obtained 
from  the  relation298 


a;2 

1  9  ’ 


(7.5) 


where  jjp  is  the  bulk  plasma  frequency  and  foMg  is  the  energy  gap  of  the  insulating 
region.  From  the  local  densities  of  states  in  Fig.  7.1,  the  energy  gap  in  the  metal- 
oxide  layer  is  =  8.3  eV  for  Lai.ssSro.isCuC^  and  7.8  eV  for  YBa2Cu3C>7.  Using 
an  estimate  hujp  ^  10  eV',  the  electron-hole  interaction  in  (7.1)  is  Veh  =  —2.1  eV 
for  the  LaO  exciton  and  Veh  =  -2.0  eV  for  the  BaO  exciton.  The  resulting 
exciton  states  are  shown  as  the  solid  curves  in  Fig.  7.2,  with  the  local  densities 
of  states  for  the  oxygen  atoms  again  shown  by  the  dashed  curve  for  reference. 
Although  the  Green’s  function  Gq(E)  in  (3.7)  is  broadened  by  a  Lorentzian 
lineshape  with  6  =  0.2  eV,  the  exciton  states  in  Fig.  7.2  have  a  larger  width  of 
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Fig.  7.2.  Densities  of  states  for  a  tightly— bound 
Frenkel  exciton  in  Lat  85Sr0 18CuO+  and  YBa2Cu30?. 
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order  1  eV,  related  to  the  bandwidth  of  the  unoccupied  d  orbitals  that  form  the 
final  state  of  the  charge  transfer  exciton. 

These  exciton  states  should  be  observable  in  electron-energy-loss  spec¬ 
troscopy  (EELS)  or  other  spectroscopies  that  create  holes  and  excited  electrons. 
Recent  EELS  measurements  on  YBa2Cu3C>7_y  show  a  rather  sharp  peak  ap¬ 
proximately  4  eV  above  the  Fermi  energy,  which  has  been  attributed  to  a  BaO 
exciton.299  The  agreement  with  the  exciton  state  in  Fig.  7.2  is  somewhat  for¬ 
tuitous  considering  the  approximations  made  for  the  binding  interaction  \eh- 
but  this  experimental  observation  appears  to  confirm  that  real  excitons  can  be 
formed  in  the  insulating  layers  of  the  high-temperature  superconductors. 

Two-Dimensional  Exciton  States 

A  more  sophisticated  treatment  of  the  exciton  states  is  needed  to  properly 
evaluate  the  matrix  elements  for  the  interaction  between  the  excitons  and  the 
charge  carriers.  Rather  than  limiting  the  exciton  to  the  tightly-bound  Frenkel 
states  of  the  previous  section,  we  now  allow  the  exciton  to  form  on  arbitrary 
metal  and  oxygen  sites  within  the  insulating  layers,  and  then  proceed  to  calculate 
the  dominant  contributions  to  the  effective  pairing  interaction  of  the  charge 
carriers. 

The  crystal  environment  of  the  insulating  metal-oxide  layers  is  essentially 
planar,  with  each  metal  or  oxygen  atom  four-fold  coordinated  with  its  nearest 
in-plane  neighbors.  Neglecting  the  displacements  along  the  c  axis  of  the  metal 
and  oxygen  atoms  in  one  layer,  we  model  each  layer  as  a  two-dimensional 
system.  This  model  is  valid  for  the  cubic  bismuth-oxide  systems  and  for  the 
hole-doped  cuprates,  which  have  small  c-axis  displacements  in  the  insulating 
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layers  e.g.,  0.57  A  for  La-0(2)  in  Lai ssSro.isCuC^  (Ref.  9)],  but  is  more 
approximate  for  the  electron-doped  systems  such  as  Nd2-xCexCu04,  which  have 
a  three-dimensional  bonding  coordination  for  the  metal  and  oxygen  atoms  in  the 
insulating  regions.86 

The  eigenstates  <f>(r,o)  for  the  two-dimensional  Coulomb  potential  may  be 
found  from  the  Schrodinger  equation 


1  did 


2 meff  I  r  dr  \  dr 


i  d2 


<T  -  -T7TT - =  E*(r,o)  ,  (7.6) 


where  meff  is  the  reduced  effective  mass  of  the  electron-hole  pair,  is  the 
dielectric  constant  characterizing  the  insulating  region,  and  (r,<p)  are  the  polar 
coordinates.  Separating  variables  by 


<E>(r,o)  =  R(r) 


gives  the  radial  equation 
'  _  1  d_ 

dr2  r  dr 


2 n?effe2  ^  2 rn,ff£ 
h2  62 r  h2 


R(r)  =  0 


where  m  =  0.  1,  2,  ....  This  equation  has  a  straightforward  solution  for  the 
bound  states  with  E  <  0. 

We  first  find  the  energy  eigenvalues  E  for  (7.8).  Defining  the  length  scale 
a  =  h2 t2/‘imeffe2  and  substituting  u(r)  =  yfr  f?(r)  gives  the  modified  radial 


equation 


(^2  -  |) 


1  2meffE 


Changing  variables  with 


u(r)  =  0  . 


(7.10) 


and  defining 
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A  =  - 


h 2  \1/2  1 


2m.fr E  /  2 a 


gives  the  differential  equation 


'  d2  {m2-\)  A  1]  ,  , 

— —  —  - - — —  +  —  —  -  ulp)  =  U  . 

dp 2  p2  2p  4j 

The  asymptotic  solution  to  (7.12)  is 

lim  i t(p)  x  e~p'2 

p— 30 

for  the  bound  states.  We  thus  seek  a  general  solution  of  the  form 

u(p)  = 

where  f(p )  satisfies 

'  d2  d  (m+  |)(m  -  |)  A  ]  n 

[ip  ~  TP - 7 - Tp\m  =  °- 

Assuming  a  polynomial  form 


(7.11) 


(7.12) 


(7.13) 


(7.14) 


(7.15) 


f(P)  = 

fe= o 


and  substituting  into  (7.15)  gives 


■  d2 

d 

/  A  -  1 

pdp 2 

-4-  (2m  -  1  -  p)—  + 

(  2  HJ 

Solving  for  the  recurrence  relation  of  the  coefficients  gives 


(7.16) 


(7.17) 


1 


(k  -  l)(Jfc  -  2m  -  1) 


(7.18) 
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For  finite  solutions  the  series  in  (7.16)  must  terminate  at  some  Armax,  so  that 
(7.18)  requires 

A  =  2 (fcmax  +  m)  4-  1  .  (7.19) 

Since  A;max,  m  =  0,  1,  2,  . . . ,  then  A  may  only  have  odd  integral  values. 
Renaming  A  as  the  principal  quantum  number  n  gives  the  energy  eigenvalues 
from  (7.11) 

En  =  -4  (  — ,ef  V  )  ~4  ,  n  =  1.  3,  5 .  (7.20) 


v  2 h2e\  )  n2 

The  two-dimensional  Coulomb  potential  thus  gives  binding  energies  four  times 
larger  than  the  three-dimensional  case. 

To  solve  for  the  eigenstates,  we  note  that  (7.17)  is  in  the  form  of  the 
differential  equation  for  the  associated  Laguerre  polynomials 


d 2  d 

zdT2  +  (p  +  1  ~:)T~  +  {q~p) 


LPS~~)  = 


dz 


dP 

dzP 


L?q{ -')  =  0  , 


(7.21) 


-  dq  ,  n  --X 

e~  -  (~?e  ') 

dz ?  v  ' 


Comparing  with  (7.17)  gives  p  =  2m  and  q  =  m  -t-  (n  —  l)/2.  The  properly 
normalized  radial  wave  function  satifying  (7.8)  is  thus  found  to  be 


Rn,m{r)  =  2 


(Sfi-m)' 


n—  1 


m 


ll3 


)•] 


{no)  6XP(  2 nJ  L^n-D/2  (na) 


n  =  1,  3,  5,  . . . 
m  =  0.  1,  2 . (n  -  l)/2  . 


(7.22) 


For  these  two-dimensional  wave  functions,  the  principal  peak  of  the  radial 


probability  density  occurs  at  a  distance  rn  —  an2  =  \(h2t2 /meffe2)nz,  which  is 
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four  times  smaller  than  the  equivalent  “Bohr  radius”  for  the  three-dimensional 
Coulomb  potential. 

To  complete  this  treatment  of  the  excitons,  we  examine  the  parameters 
for  the  model.  The  electron  in  the  exciton  is  transferred  into  an  unoccupied 
metal  d  state,  and  is  regarded  as  localized  because  of  the  large  effective  mass 
for  the  relatively  flat  d  bands.  The  hole  in  the  exciton  is  found  on  some  linear 
combination  of  oxygen  sites  neighboring  the  metal  atom.  The  oxygen  p  bands 
are  more  dispersive  than  the  metal  d  bands,  so  that  the  reduced  effective  mass 
me ff  is  essentially  the  mass  of  the  hole,  which  we  take  to  be  simply  me.  (Notice 
that  only  the  ratio  e 2/'mett  occurs  in  rn  above.)  The  exciton  wave  functions 
found  above  thus  describe  the  hole,  with  the  origin  located  on  the  metal  site 
containing  the  electron  in  a  localized  d  orbital.  The  probability  distribution 
gives  an  envelope  function  for  the  hole,  weighting  the  possible  occupation  of  the 
neighboring  oxygen  sites. 

Figure  7.3  shows  the  radial  probability  distribution,  rR\  0,  for  a  range  of 
dielectric  constants  ?2,  using  the  crystal  structure  of  Lai.gsSro.isCuCU.  Even  for 
a  dielectric  screening  as  large  as  £2  =  20,  the  peak  of  the  probability  distribution 
occurs  inside  the  first  unit  cell,  indicated  by  the  La-0(2)  bond  length  (2.73  A)  in 
Fig.  7.3.  The  s-like  exciton  ground  state  thus  gives  equal  occupation  weight  to 
the  four  nearest-neighbor  oxygen  sites  in  the  LaO  layer,  with  negligible  weight  on 
more  distant  neighbors.  Since  the  first  excited  state  has  n  =  3,  the  distribution 
for  all  but  the  ground  state  will  be  significantly  spread  to  more  distant  neighbors. 
As  described  in  the  next  section,  the  exciton  transition  density  is  negligible  for 
all  but  nearest-neighbor  atoms,  so  we  conclude  that  only  the  ground  state  of  the 
exciton  will  contribute  to  the  effective  pairing  interaction  between  the  carriers 
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Fig.  7.3.  Radial  probability  distribution  for 
two-dimensional  excitons  in  La}  MSr0 15Cu04. 
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in  the  adjacent  metallic  plane. 


Interaction  Hamiltonian  and  Matrix  Elements 


To  calculate  the  effective  interaction  between  the  charge  carriers  in  the 
metallic  plane,  we  begin  with  the  second-quantized  Hamiltonian300 

tfmet  — ins  =  \  j  d£d*  {^U'ins  (**  )' *  )^ins(®'  )0met  ( ®)  .  (7.23) 


Here  V  (x,  x1)  is  limited  to  interactions  which  occur  between  an  electron  at  x  in 
the  metallic  plane  and  one  in  the  insulating  layer  at  x1 .  Specifically,  the  Coulomb 
repulsion  between  charge  carriers  within  the  metallic  planes  has  been  excluded, 
as  well  as  interactions  between  charges  in  the  insulating  region.  The  field 
operators  in  (7.23)  are  given  in  the  metallic  region  by  the  standard  expressions 


*W*)  = 

ki  ,n 


(7.24) 


where  at  and  dr  are,  respectively,  the  creation  and  annihilation  operators 

ki ,n  *1  <n 

for  an  electron  with  momentum  k\  in  the  nth  band.  The  basis  states  can  be 
written  as  Bloch  sums 


(7.25) 


The  basis  states  are  normalized  over  the  N  unit  cells: 


/ i!  =  i  • 


(7.26) 
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In  the  following,  we  approximate  the  0n(x)  in  (7.25)  by  the  atomic  orbitals 
associated  with  a  metallic  plane. 

The  operators  in  (7.23)  that  act  on  the  insulating  region  are  written  as 


k  L  m*  771 

$  ( X ')  =  VfVv?*  (/)  ct  +  Vv^*(x')ct 

■nsv  '  ^Xk+q,m’^  >  k+q,m‘  ’  k,r 


(7.27) 


The  operators  and  c  act  only  in  the  insulating  layer,  creating  or  destroying 
electrons  in  the  basis  states 


*Lf, m.(?)  = 

R 


)  -  %/v 


(7.28) 


where  am-(x')  represents  one  of  the  five  d  atomic  orbitals  associated  with  the 
metal  atom  in  the  insulating  layer,  and  Jrn(if')  represents  one  of  the  linear 
combinations  of  p  orbitals  on  the  neighboring  oxygen  atoms.  The  momentum 
q  introduced  explicitly  in  (7.27)  and  (7.28)  will  be  found  below  to  be  the 
momentum  of  the  exciton. 

The  operators  in  the  metallic  plane  anticommute  with  those  acting  in  the 
insulating  layer.  We  can  thus  examine  the  product  of  field  operators  in  (7.23) 
for  each  region.  We  limit  the  calculation  to  the  material  Laj.gsSro.isCuO^,.  As 
mentioned  in  Chapter  IV,  EELS  data  indicate  that  the  conduction  bands  have 
primarily  oxygen  px  and  py  character  within  the  CuC>2  planes.168  We  consider 
only  the  cr-bonded  0  p  band  as  a  carrier  state,  so  that  4>n{£)  =  <p{x)  is  simply 
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an  oxygen  p  orbital.  The  operator  product  in  the  metallic  region  is  then 
t'met(^)V’met('?)  =  EE  <P  {x  ~  R\  )<t>(x  —  R\ ) 


kltk[  flt.R; 


(7.29) 


x  e 


fc;  *1 


Neglecting  the  overlap  of  atomic  orbitals  on  different  sites  reduces  (7.29)  to 


(•?)<►’««(*)  =  y  EE  IE 


i(ki-k[)'Ri  at  - 

k',  *l 


(7.30) 


ki.k[  /?! 

The  product  of  operators  in  the  insulating  region  is 

)</’ins(*  )  =  E  E  E 


k,k'  m'm'  R,R' 


i[k-R-(k'+q)-R'}  -t 

k 1 4-  <f,  m 1  ^  < m 


-  J^x1  -  R)am'(x'  -  R'y'.Ck-<n-R'-k'  R\  dt 


(7.31) 


k',m  k-q,m'  [ 


Here  we  have  dropped  the  terms  with  initial  and  final  states  on  the  same  atom 
(i.e.,  a*a  and  3" 3),  since  these  terms  do  not  lead  to  the  creation  or  destruction 
of  a  metal-oxide  exciton.  The  integrated  transition  density  for  the  exciton  is 


/ 


a;,(/)im(x') 


(7.32) 


The  oxygen  p  and  lanthanum  d  states  are  treated  as  hydrogenic  wave  functions 
with  the  Thomas-Fermi  scaling301  Zeffect,ve  =  Z1^3.  Within  this  approximation, 
all  but  nearest  neighbors  will  give  a  negligible  result  for  (7.32).  The  exciton 
states  are  therefore  limited  to  one  unit  cell,  requiring  R'  =  R  in  (7.31).  As 
found  in  the  previous  section,  the  envelope  function  for  the  exciton  ground  state 
gives  equal  weight  to  the  nearest  neighbors.  Since  the  LaO  layers  are  essentially 
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planar,  we  expect  the  out-of-plane  0  pz  orbitals  to  have  small  transition  densities 
with  the  La  d  orbitals,  so  they  are  excluded  from  the  combinations  making  up 
3m(xl).  There  are  thus  40  types  of  excitons,  arising  from  the  five  lanthanum  d 
orbitals  and  the  eight  linear  combinations  of  pz  and  py  orbitals  on  each  of  four 
adjacent  oxygen  sites.  The  sums  over  m  and  m'  in  (7.31)  can  then  be  relabeled 
to  sum  over  the  types  of  excitons  (a  =  1,  2, . . . ,  40),  giving 


!-'iTns(f'H'ins(-?') 


1 

y 


el(k-k‘-q)-R 


X  |q;(x  -  R)3&{x  -  R)  ct  +-acf  a 
-  j;(?  -  R)aa(x  -  R)  • 


(7.33) 


To  this  point,  q  is  arbitrary,  so  we  may  set  k'  =  k  and  sum  over  k  and  q.  We 
define  the  exciton  operators 
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1 
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y  4 . 

L _ -  K-rqy& 


'fc,a 


4 


C  -  Cr  - 
k,n  k-q,a. 


(7.34) 


where  it  creates  an  exciton  of  momentum  q  and  type  a  in  the  insulating  region, 

9, a 

and  6_9-a  destroys  an  exciton  of  momentum  —q  and  type  a.  We  then  write  (7.33) 
as 

=  IE  ax?  -  am?  - 1) 'b 

R 

-  j;(£,-«)a4(f‘-fl)6_,-a|  . 


(7.35) 
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The  interaction  Hamiltonian  in  (7.23)  may  now  be  written  as 


1 


*.«- 1-  =5  E  E  E  •' 


<?-a  R^R 


x  ^  j  dxdx  \<j)(x-  R\ ) |2  V{x,x')  at, 
x  <>  a;(*»  -  m(£'  -  R)  b\A  +  £(*»  -  R)a :.(*■  -  A) 


(7.36) 


To  proceed  further,  we  limit  the  form  of  the  interaction  to  V(x,  £')  =  V(x  —  x '). 
We  then  use  the  change  of  variables 

x0  =  x  -  Rx 
4  =  x  -  Rx 
l  =  R  —  R\  , 

so  that  xq  —  Xq  =  £  —  x'  ■  Substituting  into  (7.36),  we  find 


(7.37) 


=!  EEE  «-*r 


fc».*i  «’a  J?i.i 


^  J  dx0dx0  |<p(x0)|2  V(£0  -*o)  (7.38) 


x  <(4  -  O0a(*O  -  0  4, a  +  #(*0  -  Oa«(®0  “  0  b- 


9>j  * 


The  sum  over  Rx  gives  NSj- where  K  is  any  reciprocal  lattice  vector. 
If  we  neglect  Umklapp  processes,  then  K  =  0  and  k[  =  k i  —  q.  The  final 
Hamiltonian  is  then 


H, 


met 


-ins  -  2 


fc.  <?.* 


«.(?!  4..  +  i-f- 


(7.39) 


The  carrier-exciton  interaction  matrix  element  is  given  by 

Mi{q)  =  e~tq~r  J  dxdx1  \<p{x)\2V{x  -  x,)nl(x'  -  f).T(x'  -  l)  .  (7.40) 
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In  (7.40)  the  sum  over  excitons  centered  at  different  lattice  points  l  is  limited  to 
third-neighbor  unit  cells  with  respect  to  the  electron  at  x  in  the  metallic  plane, 
assuming  that  interactions  with  more  distant  excitons  will  be  fully  screened. 
The  interaction  V’(x  —  x1)  is  taken  to  be  the  Coulomb  interaction 


V(x-f)  = 


(i\X  —  x'\  ’ 


(7.41) 


where  ex  is  the  effective  dielectric  constant  that  characterizes  the  region  between 
the  metallic  and  insulating  layers.  A  bulk  static  dielectric  constant  of  ==  20  has 
been  reported302  for  insulating  La2Cu04,  which  may  be  used  as  a  crude  estimate 
for  €\ . 

The  interaction  Hamiltonian  (7.39)  is  now  in  precisely  the  form  given  by 
the  left  side  of  the  Feynman  diagram  in  Fig.  7.4.  The  exchange  of  the  exciton 
with  four-momentum  q  =  (q,  g0)  between  the  two  carriers  leads  to  the  exciton- 
mediated  effective  interaction303 


^ex(?»9o)  —  2  N  ' 

a 


'-!/a(<f)|2 

90  -  (^?>)2 


(7.42) 


For  energy  transfer  qo  <C  and  for  characteristic  excitational  energies 

hvtjq, a  ~  ^Oi  (7.42)  reduces  to 


KM 


„  r.  -V»(9ll2 
fio-o 


(7.43) 


This  is  exactly  the  form  originally  given  by  Little  for  a  generic  exciton  coupling 
mechanism,280  but  here  evaluated  specifically  for  the  oxygen  p  to  metal  d 
excitons  of  the  high-temperatui’e  superconductors.  A  significant  addition  in 
this  development  is  that  the  transition  densities  in  (7.32)  can  be  large  because 
of  the  nearest-neighbor  orbital  overlap,  so  that  a  substantial  LeX(<f)  can  be  found 
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even  for  large  excitational  energies  h<M0.  We  also  see  that  the  exciton-mediated 
pairing  potential  is  independent  of  the  sign  of  the  charge  carriers  because  of 
the  factor  -l/a(<?)|2  in  (7.43),  so  the  present  mechanism  can  be  applied  to  both 
n-type  and  p-type  carriers. 

The  six-dimensional  integral  of  (7.40)  was  evaluated  using  Gaussian  quadra¬ 
tures,  with  afactor  1  —exp  (  —  Ajx  —  5?  |2)  applied  to  (7.41)  to  avoid  the  divergence 
at  x  —  x1 .  A  value  A  =  10  A-2  was  used,  creating  a  correlation  screening  radius 
of  0.3  A  about  each  charge.  The  results  were  unchanged  for  larger  values  of  A 
(corresponding  to  shorter  screening  lengths).  The  results  for  V'e!C(<f)  are  given  in 
Table  7.1  for  q  parallel  to  the  Cu02  planes,  using  the  parameters  ei  =  20  and 
fru>o  =  7  eV.  The  largest  contributions  occur  for  q  =  0,  as  expected  from  the 
Coulomb-like  interaction  between  the  carrier  and  the  exciton. 

For  the  crystal  structure  of  La1.85Sro.i5Cu04,  there  are  two  identical  LaO 
layers  sandwiching  the  Cu02  plane,  so  that  \l&(q)  in  (7.40)  must  also  include 
a  factor  2 cos  (qzdz)  for  the  coherent  sum  over  these  two  layers,  where  dz  is 
the  distance  in  the  r  direction  between  the  metallic  and  insulating  layers  (e.g., 
1.9  A  between  the  Cu02  and  LaO  layers).  The  resulting  Vex(<f)  is  then  simply 
4[cos {qzdz)\2  times  the  values  in  Table  7.1.  We  note  that  in  the  materials  with 
multiple  Cu02  layers  per  unit  cell  (such  as  Tl2Ca2Ba2Cu30io),  there  is  only  a 
single  insulating  metal-oxide  layer  adjacent  to  the  outer  Cu02  plane,  so  that  no 
summation  occurs  for  excitons  in  different  insulating  layers. 

The  V'eX(g)  found  within  this  model  may  be  used  to  further  develop  the 
details  of  this  mechanism  for  high-temperature  superconductivity.  For  example, 
a  crude  estimate  for  the  transition  temperature  can  be  made  using  the  expression 
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TABLE  7.1.  Exciton-mediated  effective  interaction  Vex(q) 
for  La1.83Sro.i5Cu04,  in  eV. 

qx  =  0  7r/6a  W4a  7r/3a  7r/2a  2tt/3  a  3t/4  a  5ir/6a  ir/a 

qy  =  0  -1.00  -0.84  -0.67  -0.49  -0.19  -0.05  -0.02  -0.01  0.00 

=  7T  '6a  -0.88  -0.73  -0.58  -0.41  -0.15  -0.03  -0.01  0.00  0.00 

qy  =  7T  '4a  -0.74  -0.61  -0.47  -0.33  -0.11  -0.02  -0.01  0.00  0.00 

qy  =  rr  3a  -0.58  -0.47  -0.36  -0.24  -0.07  -0.01  0.00  0.00  0.00 

qy  =  ir>2a  -0.27  -0.20  -0.14  -0.08  -0.01  0.00  0.00  0.00  0.00 

qy  =  2tt  3a  -0.08  -0.05  -0.03  -0.01  -0.01  -0.02  -0.02  -0.01  0.00 

qy  =  3rr/4a  -0.04  -0.02  -0.01  -0.01  -0.03  -0.04  -0.04  -0.02  0.00 

qy  =  ore /6a  -0.02  -0.01  -0.01  -0.01  -0.05  -0.06  -0.05  -0.03  0.00 


of  Emery  and  Reiter304 


re  -  Ef  e-*'  x"  ,  (7.44) 

where  Aex  =  .V(0)V'ex  .  This  expression  is  appropriate  for  nonretarded  electronic 
interactions.304  For  .Y(0)  ~  0.5  states  ,/eV  cell,  an  average  Vex(q)  of  roughly 
—  1  eV  (including  the  q,  dependence  times  the  results  listed  in  Table  7.1), 
and  an  estimate  of  Ef  =  0.15  eV  for  Lai  gsSro.isCuC^,304  we  obtain  Tc  -~ 
200  K.  This  is  certainly  not  a  quantitative  result,  but  simply  indicates  that 
the  excitonic  mechanism  described  above  can  provide  the  pairing  necessary  for 
superconductivity  within  a  physically  reasonable  range  of  parameters. 

Summary  of  the  Excitonic  Mechanism 

The  specific  structural  and  electronic  configuration  required  for  this  exci¬ 
tonic  mechanism  of  high-temperature  superconductivity  has  been  observed  in 
all  currently-known  high-T).  materials,  including  the  hole-doped  cuprates  such 
as  La^goSro.isCuO-i,  YBazCugOy,  B^CaS^C^Og,  and  Tl2Ca2Ba2Cu30io;  the 
electron-doped  cuprates  such  as  Nd2_xCexCu04;  and  the  copper-free  system 
Bai_xKxBi03.  In  every  case,  metallic  layers  are  found  in  close  proximity  to 
insulating  layers,  so  that  excitons  within  the  insulating  region  can  interact  with 
carriers  in  the  metallic  region.  Excitons  in  the  predicted  energy  range  have 
apparently  been  observed  in  the  BaO  layers  of  YBa2Cu307_y .2" 

For  reasonable  estimates  of  the  physical  parameters,  this  mechanism  pro¬ 
vides  a  pairing  interaction  sufficient  to  give  high-Tc  superconductivity.  No 
specific  magnetic  properties  are  required  for  this  mechanism.  We  predict 
that  high-temperature  superconductivity  will  not  be  seen  in  materials  such  as 
Cai_xKxCu02  or  Cai_xYxCu02,  since  these  structures,  although  containing 
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Cu02  planes,  have  no  insulating  metal-oxide  layers  to  support  the  excitons.82 
Quantitative  predictions  within  this  theory  will  require  a  comprehensive  treat¬ 
ment  of  the  superconductivity,  including  the  effects  of  anisotropy  and  large  ex- 
citational  energies  associated  with  the  excitons. 
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CHAPTER  VIII 

CONCLUSION 

The  characterization  of  high-temperature  superconductors  is  a  significant 
experimental  and  theoretical  challenge.  The  present  research  has  concentrated 
on  a  theoretical  model  for  the  electronic  properties  of  these  materials.  We  have 
also  proposed  an  excitonic  mechanism  that  is  applicable  to  all  currently-known 
high- temperature  superconductors. 

In  Chapter  III  vve  developed  the  semiempirical  tight-binding  model  which 
was  used  to  calculate  the  electronic  energy  bands,  the  local  and  total  densities  of 
states,  and  the  atomic  valences  for  each  of  the  materials.  The  results  of  the  cal¬ 
culations  were  presented  in  Chapter  IV  for  the  copper-oxide  superconductors 
La185Sro.i5Cu04,  YBa2Cu307,  Bi2Sr2Cu06,  Bi2CaSr2Cu2Os,  Tl2Ba2Cu06, 
Tl2CaBa2Cu208,  Tl2Ca2  Ba2Cu3Oio,  and  TlCasBa^'u.jOn,  and  for  the  bis¬ 
muth  oxides  BaPbo.75Bio.25O3  and  Ba0.6Ko.4Bi03.  We  find  that  a  single  tight- 
binding  model,  with  fully  transferable  parameters,  provides  a  good  description 
of  the  electronic  structure  for  each  of  these  high-temperature  superconductors. 

All  the  copper-oxide  superconductors  have  two-dimensional  conduction 
bands  arising  from  Cu(d)-0(p)  antibonding  states.  The  metal-oxide  lay¬ 
ers  sandwiching  the  Cu02  planes  (such  as  LaO,  BaO,  and  SrO)  all  display 
ionic  and  insulating  behavior.  The  copper-free  bismuth-oxide  superconductors 
BaPbo.75 Bi0  25O3  and  Ba0, 6K04BiO3  have  valence  bands  of  antibonding  Bi(s)- 
0 (p)  states,  with  cubic  or  almost-cubic  symmetry.  These  materials  also  have 
ionic  and  insulating  layers  of  BaO  in  their  structure.  Although  all  of  these 
systems  have  very  different  crystal  and  magnetic  structures,  the  common  fea- 
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tures  in  the  electronic  structures  provide  support  for  a  single  mechanism  of 
high- temperature  superconductivity. 

Chapter  V  reported  the  study  of  atomic  substitutions  in  YBa2Cu307, 
B^CaS^C^Og,  and  Tl2CaBa2Cu20g.  There  is  reasonable  agreement  with 
expected  chemical  trends  and  with  experimentally-observed  changes  in  the 
superconducting  properties.  The  results  for  substitution  of  Al,  Fe,  Co,  Ni,  and 
Zn  for  Cu  in  YBa2Cu3C>7  show  that  substitutions  on  the  Cu(2)  site  disrupt 
the  Cu02  plane  conduction  bands  more  than  do  substitutions  on  the  Cu(l) 
site.  Results  for  the  replacement  of  Ba  by  Sr  in  YBa2Cu3C>7  support  a  purely 
structural  effect  for  the  suppression  of  Tc.  We  find  that  La  substitution  for  Ba 
leads  to  destruction  of  hole  charge  carriers  by  the  donated  electrons,  resulting 
in  direct  suppression  of  superconductivity.  Substitution  of  F  or  N  on  the  oxygen 
sites  shows  that  the  electronic  structure  of  YBa2Cu307  may  be  slightly  affected 
by  such  replacements,  but  no  obvious  changes  in  the  superconducting  properties 
are  correlated  with  these  results. 

Substitution  of  Pb  for  Bi  in  Bi.l/CaS^C^Og  gives  an  increase  in  the  num¬ 
ber  of  hole  carriers  within  the  Cu02  planes,  assuming  the  same  idealized  struc¬ 
ture  for  BiPbCaSr2Cu20g  and  B^CaS^C^Og.  The  observed  enhancement  of 
Tc  with  Pb  substitution  in  the  bismuth  cuprates  has  to  date  been  correlated 
with  the  stabilization  of  the  higher- Tc  phase  Bi2Ca2Sr2Cu30io,  making  a  direct 
comparison  with  these  results  somewhat  difficult.  The  strongest  effect  of  doping 
with  Hg  or  Pb  in  Tl.\/CaBa2Cu20g  is  the  shift  in  the  Fermi  energy,  with  Hg 
creating  holes  and  Pb  destroying  them. 

Chapter  VI  presented  the  results  for  the  electronic  effects  of  oxygen  vacan¬ 
cies  in  Lai.g5Sro.i5Cu04 _y  and  YBa2Cu3C>7_y,  and  of  lanthanum  vacancies  in 
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La2-xCu04.  The  modification  of  the  density  of  states  and  the  shift  in  the  Fermi 
energy  were  calculated  for  each  case.  Oxygen  vacancies  act  as  electron  donors  in 
both  systems,  raising  the  Fermi  energy  and  decreasing  p{Ep)  for  y  <  0.5.  The 
concentration  of  holes  is  thus  decreased  for  the  oxygen-defect  systems,  consistent 
with  the  observations  that  oxygen  vacancies  degrade  the  superconductivity  and 
metallic  conductivity  in  these  materials.  Lanthanum  vacancies  in  La2-xCu04 
donate  holes,  lowering  the  Fermi  energy  and  increasing  p{Ep)  for  x  <  0.3.  Small 
concentrations  of  La  vacancies  may  thus  have  the  same  effect  electronically  as 
divalent  metal  dopants  or  excess  oxygen,  leading  to  enhancement  of  supercon¬ 
ductivity  in  nominally  stoichiometric  La2Cu04. 

Chapter  VII  presented  a  specific  mechanism  for  high-temperature  super¬ 
conductivity  which  requires  tightly-bound  excitons  in  the  insulating  metal-oxide 
layers  adjacent  to  the  superconducting  planes.  Excitons  with  energies  in  the  pre¬ 
dicted  range  appear  to  have  been  observed  in  the  BaO  layers  of  YBa2Cu307-y. 
No  specific  magnetic  structure  for  the  metallic  region  is  needed  for  this  mecha¬ 
nism,  making  it  equally  applicable  to  the  copper  oxides  and  the  bismuth  oxides. 
The  required  structural  and  electronic  configuration  has  been  observed  in  all 
currently-known  high-Tc  materials,  including  the  hole-doped  cuprates  such  as 
La!  85Sr0.i5CuO4,  YBa2Cu307,  B^CaS^C^Og,  and  Tl2Ca2Ba2Cu30!o;  the 
electron-doped  cuprates  such  as  Nd2-xCexCu04,-  and  the  copper-free  system 
Ba!_xKxBi03.  In  every  case,  metallic  layers  are  found  adjacent  to  insulating 
layers,  so  that  carriers  in  the  metallic  region  can  experience  a  net  attractive 
interaction  mediated  by  the  excitons.  A  rough  estimate  with  reasonable  val¬ 
ues  of  the  physical  parameters  gives  Tc  200  K.  This  mechanism  predicts 
that  high-temperature  superconductivity  will  not  be  seen  in  systems  such  as 


Cai_xKxCu02  or  Cai_xYxCu02,  since  these  layered  C11O2  structures  do  not 
have  the  insulating  metal-oxide  layers  necessary  to  support  the  excitons. 
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